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QUANTITATIVE  RAY  METHODS  FOR  SCATTERING  OF  SOUND  BY 

SPHERICAL  SHELLS 

Abstract 

by  Steven  Gregory  Kargl,  Ph.  D. 

Washington  State  University 
August  1990 

Chair:  Philip  L.  Marston 

The  application  of  ray  methods  to  the  scattering  of  high-frequency  plane  waves 
from  evacuated  elastic  spherical  shells  is  investigated.  The  investigation  of  ray  methods  for 
spherical  shells  is  a  precursor  to  the  application  of  such  methods  to  shells  having  more 
complicated  shapes.  The  scattered  pressure  in  the  faifield  of  the  shell  is 
Psc  =  Pi(P/^f‘)fi^)cxp(.ikr)  where  p,-  is  the  plane  wave  amplitude.  The  outer  radius  of  the 
shell  is  a,  />  =  0  -  b  is  the  shell's  thickness,  and  r  is  the  distance  to  an  observation  point. 

Ray  models  are  developed  to  synthesize  the  form  function /(0,  ka)  where  k  is  the 
wavenumber  of  the  incident  wave  and  0  is  the  scattering  angle.  The  forward  scattering 
amplitude, /(9  =  0),  is  related  to  the  extinction  cross  section,  o^,  by  the  optical  theorem.  If 
the  absorption  by  the  scatterer  is  negligible,  then  <7e  is  equal  to  the  total  scattering  cross 
section  0/.  A  ray  synthesis  partitions  f(d  =  0)  into  a  component  for  ordinary  forward 
diffraction  about  the  shell,//:©,  and  contributions  from  surface  guided  elastic  waves.  For 
high-frequency  scattering,  the  relevant  surface  guided  elastic  waves  are  leaky  Lamb  waves. 

A  similar  ray  synthesis  of  the  backscattering  amplitude y(0  =  n)  contains  a  specular 
reflecdon  component, /^^(S  =  7t),  and  leaky  Lamb  wave  contributions.  A  generalization  of 
the  geometrical  theory  of  diffraction  is  employed  to  synthesize //(9  =  0,  ka)  and 
//(0  =  7t,  ka)  for  the  Uh  leaky  Lamb  wave  contribution.  The  syntheses  for  forward  and 
backwards  scattering  correctly  describe  the  leaky  Lamb  wave  contributions  and  are 
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expressible  in  a  Fabi>'-Perot  resonator  form.  While  the  ray  description  uf  backscattering 
ordinarily  accurately  reproduces  exact  computations  and  experiments  with  tone  ba^st, 
certain  anomalies  are  discussed,  A  ray  synthesis  of/^,  demonstrates  a  significant 
longitudinal  resonance  effect  when  kih  =  tin,  n  =  1,  2, vyhere  ki^  =  oaici  is  the 
longitudinal  wavenumber  within  the  shell.  The  analysis  of  fsp  is  for  an  elastic  mattrial 
with  vanishing  shear  velocity.  The  relevant  range  of  is  7  <  to  <  100.  The  rhell  is 
surrounded  by  water  and  is  composed  of  440c  stainless  steel  with  i  iner-to-outer  radii  ratio 
b/a  =  0.838. 
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Chapter  1 
Introduction 

1.1  Review  of  the  exact  partial  wave  series  for  the  form  function 

The  application  of  acoustic  ray  methods  to  scattering  problems  can  provide  a  simple 
understanding  of  the  physical  interaction  of  an  incident  pressure  wave  and  a  scatterer. 

When  a  scatterer  has  an  arbitrary  shape,  direct  solution  of  the  appropriate  boundary  value 
problem  is  not  usually  possible.  While  numerical  solutions  can  be  performed  for  a  scatterer 
of  known  composition  and  shape,  physical  insight  into  the  underlying  interaction  between 
an  incident  wave  and  a  scatterer  is  often  obfuscated.  Ray  methods  may  be  directly  applied 
to  describe  the  scattering  from  an  object  of  arbitrary  shape  and  the  physical  mechanisms 
involved  in  the  scattering  process  may  easily  be  identified  and  understood.  When 
developing  new  ray  techniques,  it  is  advantageous  to  test  these  techniques  by  comparison 
to  a  canonical  problem  where  an  exact  solution  is  known.  The  comparison  of  a  ray  model 
with  an  exact  solution  should  allow  one  to  identify  regions  where  the  simple  ray  model  is 
useful.  This  dissertation  investigates  novel  ray  techniques,  based  on  an  elastic 
generalization  of  the-  geometriGsl  theory  of  diffraction,^’^  which  describe  the  forward  and 
backwards  scattering  from  evacuated  elastic  spherical  shells.  It  is  anticipated  that  these  new 
ray  techniques  may  be  easily  generalized  to  elastic  objects  with  sm'xith  convex  shapes. 

Before  developing  quantitative  ray  techniques  for  the  scattering  of  a  plane  wave 
from  a  submerged  evacuated  elastic  spherical  shell,  it  is  appropriate  to  review  the  exact 
solution  to  this  canonical  problem.  A  plane  wave  in  the  surrounding  water  propagates  in 
the  positive  z-direction  and  is  incident  on  an  elastic  spherical  shell  situated  at  the  origin  of  a 
coordinate  system.  The  physical  parameters  that  describe  the  shell  are  the  inner-to-outer 
radii  ratio  bla,  the  longitudinal  sound  speed  cl.  the  shear  sound  speed  Ci,  and  the  density 
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Pc.  The  parameters  for  water  are  the  speed  of  sound  c  and  density  p.  After  scattering,  the 
total  pressure  in  the  water  isp=  pinc  +  Psc  where  pinc  =  ptexpiikz)  is  the  incident  plane 
wave  and  p^c  is  tlie  scattered  pressured.  The  pressure  amplitude  of  pi^c  is  Pi  and  k  =  2ic/k 
is  the  wavenumber  of  pinc  in  water  where  X  is  the  wavelength.  The  harmonic  time 
dependence,  exp(-ia)t),  has  been  suppressed  and  co  is  the  angular  frequency.  Using 
elasticity  theory,  Goodman  and  Stem  derived  an  exact  expression  for  the  steady-state 
pressure  as  a  partial  wave  serics.^'"^  Their  analysis  considered  an  elastic  spherical  shell 
with  a  PiUid  interior  which  gave  art  expression  similar  to  (see  Eqs.  (4a)  and  (6)  of  Ref.  3) 

P  "  Pi  X  h^*\kr)JP„(cos(e)) ,  (1) 

where  and  are  spherical  Bessel  and  Kankel  functions  of  the  fu  st  kind  and  is  a 
Legendre  polynomial.  The  functions  Bn  and  Dn  of  the  dimensionless  size  parameter  ka  are 
6x6  determinants  for  a  shell  with  a  fluid  interior  and  5x5  determinants  for  an  evacuated 
shell.  Tlie  elements  of  these  determinants  are  contained  in  Refs.  3  and  5.  The  first  term  in 
the  square  brackets  of  Eq.  (1)  is  the  partial  wave  expansion  of  pinc  while  the  second  term  is 
Psc-  In  the  farfield  of  the  elastic  shell,  Psc  has  the  spherically  diverging  formf^ 

Psc  =  Pi(a/2r)f(0)e‘“  ,  (2) 

where  the  distance  from  the  center  of  the  shell  to  a  distant  observation  point  is  r  and /(0)  is 
the  scattering  form  function  or  scattering  amplitude.  [In  general,/ is  a  function  of  the 
spherical  coordinates  0  and  (]),  but  the  azimuthal  symmetry  of  the  spherical  shell  implies  that 
/is  independent  of  ([).]  Inspection  of  Eq.  (2)  reveals  that  any  structure  evident  in  I  Psc  I  must 
be  a  result  of  structure  in  I  /(0)  I.  An  exact  partial  wave  series  representation  of /(0)  in  the 
farfield  can  be  obtained  from  Eq.  (1)  by  introducing  the  asymptotic  form 


hn  ->  i~^^^^^exp{ikr)lkr,  ifer  »  1.  Comparing  the  asymptotic  form  of  the  second  term  in 
Eq.  (1)  with  Eq.  (2)  gives 
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f(e)  = 


_2_ 

ika 


X<2n  +  1) 


n  =  0 


B„(ka) 

D„(ka) 


Pn(cos(e))  . 


(3) 


It  is  well  known  that  a  partial  wave  series  tends  to  be  slowly  converging  for  large  ka  (see 
Secs.  2.5  and  4.3  and  Appendix  B  of  Chapter  3).  Also,  from  inspection  of  Eq.  (3),  it  is 
evident  that  the  underlying  physical  interaction  of  pi^c  with  the  scatterer  is  obfuscated.  That 
is,  the  complicated  partial  wave  series  representation  of /is  not  amenable  to  a  simple 
physical  interpretation  of  the  scattering  process. 

In  1964,  Hickling  calculated  the  form  function  for  backscattering,/(0  =  it),  from 
iron  and  aluminum  shells.^  The  range  of  frequencies  investigated  corresponds  to  0  <  to  < 
30  and  the  thickness  of  these  shells  varied  from  0.05  ^  hla  <  0.8.  These  form  function 
calculations  contain  complicated  structure  tliat  was  attributed  to  the  elastic  response  of  the 
shell.  Using  Fourier  transform  techniques,  Hickling  then  consmicted  the  backscattered 
pressure  for  short  sine  wave  tone  bursts  incident  on  the  shells.  [Equations  (1)  and  (3) 
describe  the  steady-state  pressure.  As  discussed  by  Hickling,  these  equations  may  be  used 

to  triv#*  a  of  thf*  cr*attF.rino  of  tont»  Hiirctc  6t  nPh#*  oomniitf^H  Kar*lrcr»Qtt#»T#»H 

echoes  contained  a  contribution  due  to  a  specular  reflection.  The  specular  reflection  is 
essentially  the  reflection  of  the  incident  burst  from  the  outer  surface  of  the  shell.  Also, 
evident  in  the  computed  backscattered  pressures  are  echoes  that  are  delayed  in  time  with 
respect  to  the  specular  reflection.  Hickling  concluded  that  these  observed  echoes  were  due 
to  a  flexural  type  wave  within  the  shell  that  propagated  around  the  shell  re-radiating  energy 
into  the  water.  Later,  Diercks  and  Hickling  gave  experimental  confirmation  of  the 
computed  echoes  for  backscattcring  from  air-  or  water-filled  aluminum  shells.”^  Their 
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experiments  involved  several  shells  of  various  thicknesses  and  tone  bursts  of  varying 
duration  with  a  frequency  corresponding  to  ^  =  20. 

Subsequently,  for  high-frequency  .scattering  the  waves  that  circumnavigate  within 
an  elastic  shell  have  been  classified  as  leaky  Lamb  waves.^>^*^  [Franz  waves  (creeping 
waves)  may  be  guided  by  the  shell's  surface.  Franz  waves  circumnavigate  the  shell,  but 
for  these  types  of  waves  the  strain  is  primarily  in  the  surrounding  water  and  the  waves  are 
usually  heavily  radiation  damped.]  Leaky  Lamb  waves  may  be  thought  of  as  a 
generalization  of  Lamb  waves  in  a  flat  elastic  plate  in  vacuumi^'^  ^  to  the  case  of  a  fluid- 
loaded  curved  plate.  That  is,  the  particle  displacements  for  leaky  Lamb  waves  in  a  fluid- 
loaded  shell  of  thickness  k  =  a-hdst  analogous  to  the  displacements  associated  with  tlie 
equivalent  Lamb  waves  in  a  flat  plate  of  the  same  thickness.  Figure  1  illustrates  the  surface 
displacements  of  the  lowest  antisymmetric  and  symmetric  Lamb  waves  in  a  plate  in 
vacuum.  While  the  amplitudes  have  been  exaggerated,  the  values  of  kih  =  kaiclci){\  -  bla) 
used  for  Fig.  1  are  indicative  of  those  encountered  in  the  experiments  of  Ref.  5  and  those 
reviewed  below. 

1.2  Back.scattering  of  short  tone  bursts  from  a  shell 

This  section  briefly  reviews  a  ray  model  for  backscattering  from  an  elastic  sphere  or 
shell  and  supporting  experiments  from  Ref.  5.  The  scattering  process  is  illustrated  in  Fig. 
2.  Ray  diagrams  similar  to  Fig.  2  are  described  in  detail  in  Refs.  1,  2,  and  5.  The  center 
of  a  sphere  or  shell  of  outer  radius  a  is  located  at  the  origin  of  a  coordinate  system.  [The 
inner  radius  of  the  shell  is  b  and  for  the  solid  sphere  b  =  0.]  The  total  scattered  pressure  is 
a  superposition  of  a  specular  reflection  from  the  region  near  point  C  and  contributions  due 
to  the  circumnavigation  of  the  Uh  class  of  surface  guided  elastic  waves  around  the  shell. 

For  high-frequency  scattering  from  elastic  spheres,  Rayleigh  and  Whispering  Gallery 


Fig.  1  Surface  displacements  of  the  oq  antisymmetric  (a)  and  symmetric  (b)  Lamb  wave 
on  a  flat  elastic  plate  of  thickness  h  in  vacuum.  The  dashed  lines  ind’cate  the  equilibrium 
position  for  the  surfaces  of  the  flat  plate. 
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waves^'2  are  the  relevant  surface  guided  elastic  waves  Aihile  leaky  Lamb  waves^-^*^  are  the 
appropriate  elastic  waves  for  a  shell.  The  incident  plane  wave  excites  a  surface  guided 
elastic  wave  near  point  B  in  Fig.  2  that  propagates  about  the  shell  radiating  energy  back  into 
the  water.  In  the  vicinity  of  point  B',  the  radiation  is  directed  in  the  backscattered  direction. 
The  points  B  and  5'  are  determined  from  the  phase  velocity  trace  matching  condition 
sin{^l)  =  dci  where  6/  is  the  local  angle  of  incidence  and  c/  is  the  surface  guided  elastic 
wave  phase  velocity.  ^2.5  'I'he  contribution  to  the  form  function  due  to  the  backscatteiing 
of  a  surface  guided  elastic  wave  can  be  synthesized  by  considering  the  propagation  path 
length  relative  to  an  exit  plane  perpendicular  to  the  z-axis  through  C  in  Fig.  2  and  it  can  be 
expressed  in  a  Fabry-Perot  resonator  form L2, 12 


-2(jt-0,)B,  in,  -2«P,  i2rtkac/c, 

f^  =  [-G^e  ^  ^^^e  ’‘]/[l+e  c  ‘], 


(4) 


where  T)/  is  a  propagation  phase  delay.  A  complex  coupling  coefficient  Gi  describes  the 
efficiency  of  the  coupling  between  the  acoustic  wavefield  in  water  and  the  surface  guided 
elastic  wave  and  p;  is  the  radiation  damping  parameter  (in  Npfrad)  for  the  surface  guided 
elastic  wave.  A  virtually  exact  expression  for  G/  is  available  for  elastic  spheres,  but  its 
dependence  on  the  physically  relevant  parameters  c’Jc  and  P/  is  non-trivial.^  Marston  has 
developed  an  approximation  for  G/  that  appears  to  give  excellent  results  for  elastic  spheres 
as  well  as  shells  (see  Appendix  I).L13 

In  backscattering  and  near  backscattering  experiments  described  in  Ref.  5,  a  3  or  4 
cycle  sine-wave  tone  burst  was  incident  on  a  sphere  or  shell  with  frequency  co/2tt.  Figure 
3a  contains  hydrophone  output  voltage  versus  time  records  for  the  backscattering  of  a  burst 
from  a  solid  tungsten  carbide  sphere;  while  Figs.  3b  -  d  are  for  the  backscattering  from  a 
hollow  440c  stainless  steel  shell.  The  radius  of  the  tungsten  carbide  sphere  is  ^wc  =  12.7 
mm  and  the  outer  radius  of  the  440c  stainless  steel  shell  is  Oss  =  19.05  mm  while  the  shell's 
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Fig.  2  Ray  diagram  of  a  surface  guided  elastic  v/ave  contribution  to  backscattcring  from  a 
sphere  ip  =  0)  or  shell  of  inner-to-outer  radii  ratio  bla.  The  incident  plane  wave  launches  a 
surface  guided  elastic  wave  on  shell  at  B.  The  surface  wave  repeatedly  circumnavigates  the 
sphere  or  shell  continually  shedding  energy  back  into  the  water.  At  B',  the  radiated  energy 
is  directed  in  the  backscattered  direction  and  appears  to  originate  from  a  virtual  ring-like 
source  generated  by  rotating  the  point  F/  about  the  CC  axis.  The  ring-like  source 
describes  a  toroidal  wave  front  associated  with  the  axial  focusing  or  glory  scatteritig  of  the 
surface  wave.  There  are  also  contributions  to  backscattering  associated  witli  the  specular 
reflection  and  repeated  internal  reflections  in  the  vicinity  of  C. 
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Fig.  3  (Upper)  Measured  echoes  for  oackscattering  a  3  or  4-cyle  sinowave  tone  burst 
from  a  solid  tungsten  carbide  sphere  (a)  and  an  evacuated  spherical  shell  composed  of  440c 
stainless  steel  (b  -  d).  The  specular  reflection  is  labeled  S.  Rayleigh  wave  and  lowest 
antisymmetric  and  symmetric  leaky  Lamb  waves  echoes  are  denoted  by  R,  oq,  and  .vq, 
respectively.  The  tail  immediately  following  the  specular  reflection  in  (b)  and  (c)  is 
partially  due  to  a  longitudinal  resonance.  (Lower)  Schematic  of  the  experimental  apparatus 
employed  for  backscattering  and  near  backscattering  experiments. 


thickness  ish  =  ass  -b  =  3.09  mm.  The  size  paranneter,  ka  =  (Mi/c,  for  each  trace  is  as 
follows:  (a)  43.2,  (b)  64.7,  (c)  68.8,  and  (d)  36.4.  In  Fig.  3,  the  specular  reflection  from 
tivc  sphere  or  shell  is  labeled  S.  The  Rayleigh  wave  contributions  radiated  from  the 
tungsten  carbide  sphere  are  labeled  by  I  =  R  while  the  lowest  antisymmetric  and  symmetric 
leak)'  Lamb  wave  contributions  radiated  from  the  shell  are  designated  /  =  ao  and  .s  Oi 
respectively.  The  identification  of  a  specific  surface  guided  elastic  wave  with  an  echo  in 
Fig.  3  is  in  accordance  with  their  arrival  time  relative  to  the  specular  reflection.  Detailed 
analyses  of  backscattering  records  from  a  tungsten  carbide  sphere  similar  to  Fig.  3a  are 
given  in  Refs.  12  and  14.  Altiiough,  Figs.  3a  and  3d  are  similar  to  the  scattered  pressure 
computed  by  Hickling  for  an  incident  tone  burst,^d5  the  records  in  Figs.  3b  and  3c  contain 
a  feature  not  evident  in  Hie kling's  computations  for  shells.  Inspection  of  Figs.  3b  and  3c 
shows  a  long  decaying  pulse  train  immediately  following  the  specular  reflection.  This 
decaying  echo  appears  to  be  partially  associated  with  the  "ringing  down"  of  a  bulk, 
longitudinal  resonance  that  is  investigated  in  detail  in  Chapter  4. 

Distinct  echo  contributions  are  enumerated  /n  =  0, 1, ...,  by  considering  the  number 
of  complete  circumnavigations  associated  with  each  echo.  The  backscattered  pressure 
amplitude  for  the  mth  distinct  contribution  from  the  Ith  class  of  surface  guided  elastic  wave 
has  the  predicted  form 

-2(7C-0,)p, -27trnp,  ^  ^ 

IP^I  =  IP|IA^a/2r,  A^  =  IG,le  ‘ i.yu)  -  (5a, b) 

where  the  functions  Joiu)  and  Ji{u),  u  =  kcr/c  'ji,  are  cylindrical  Bessel  functions  and  yis 
the  local  backscattering  angle,  y  =  7t  --  0.  The  normalization  in  Eq.  (5a)  is  such  that  the 
specular  reflection  from  a  large  immovable  rigid  sphere  is  a  wave  whose  amplitude  is 
I  Pi  \(af2r).  For  exact  backscattering  y=  0  and  I  Joiu)  -  i^pfJiiu)  I  =  1.  Figure  4  contains 
the  normalized  on-axis  amplitudes  of  the  m  =  0  echo  for  the  I  =  ao  (solid  line)  and  ao 
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Fig.  4  Theoretical  normaiized  backscattering  amplitudes  for  oq  (solid  line)  and  .vo  (dashed 
line)  leaky  Lamb  wave  echoes  are  compared  with  experimental  data.  The  normalization  is 
with  respect  to  a  tungsten  carbide  sphere.  The  on-axis  hydrophone  distances  for  each  data 
set  from  the  point  C  in  Fig.  2  is  indicated.  The  data  points  were  obtained  from  voltage 
measurements  on  traces  similar  to  those  in  Fig.  3. 
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(dashed  line)  leaky  Lamb  waves.  Since  absolute  pressure  amplitude  measurements  could 
not  be  achieved  over  the  entire  range  of  frequencies,  both  the  theoretical  and  experimental 
amplitudes  of  the  backscattering  from  the  440c  stainless  steel  shell  were  normalized  by  the 
specular  reflection  from  the  tungsten  carbide  sphere.  The  normalized  measured  Aqi  are 
given  by  (Voi/^sv/cXdwcUsv^c  l/<3ss)  where  Vq/  and  V^wc  are  the  measured  peak-to-peak 
voltages  of  the  (/  =  oq,  sq)  leaky  Lamb  wave  echo  and  specular  reflection  from  tungsten 
carbide,  respectively.  The  magnitude  of  the  specular  reflection  form  function  for  the 
tungsten  carbide  sphere  is  I/5WC I  and  the  ratio  of  radii  appears  because  of  the  difference  in 
the  size  of  the  sphere  and  shell  used.  [For  the  frequency  range  of  interest  I  /^wc  I  =  iPe^L  ~ 
pc)/(PeCL  +  pc)  ~  0,969.]  The  agreement  for  ka  <  44  between  Eq.  (5b)  and  experiment  is 
good,  but  for  thtka>44  the  measured  values  deviate  substantially  from  the  predicted 
amplitudes.  One  plausible  cause  for  this  discrepancy  is  the  possible  excitation  of  the  /  =  fli 
and  leaky  Lamb  waves  for  ka>4\  and  65,  respectively.  These  waves  were  neglected  in 
the  original  analysis  in  Ref.  5.  A  second  probable  cause  of  the  disagreement  between  the 
measured  and  predicted  values  is  that  the  oq  and  50  echoes  overlap  in  time  for  50  <  <  70 

and  the  long  decaying  pulse  train  for  ka  >  65  may  also  interfere  with  the  ao  echo. 

Equation  (5b)  indicates  that  A^i  should  have  a  maximum  for  y  =  0  and  it  decreases 
when  the  observation  point  is  moved  away  from  the  axis  (so  that  y  0).  The  localization 
of  the  pressure  amplitude  near  the  backward  axis  is  a  manifestation  o^  the  acoustical 
glory.54'*  The  m  =  0  amplitude  for  the  /  =  oq  leaky  Lamb  wave  was  studied  to  verify  this 
prediction.  Again,  the  theory  and  measurements  were  normalized  by  the  on-axis  specular 
reflection  from  a  tungsten  carbide  sphere.  From  Eqs.  (2)  and  (5a),  the  normalized  pressure 
is 

Ip  ,(r„,y)l  A  , 

^tnr  0  _ _ ^_a _ ni 

[  1  +  tan^y] 


(6) 
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Fig.  5  Experimental  confirmation  of  the  axial  focusing  predicted  by  Eq.  (6).  The  m  =  0 
amplitude  of  the  ao  leaky  Lamb  wave  echo  was  measured  as  the  hydrophone  was  shifted 
away  from  backscattcring  (y  =  0).  For  the  440c  stainless  steel  shell  studied,  the  frequency 
of  the  incident  4-cycle  sine-wave  tone  burst  conesponds  to  ka  =  24.3 
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The  on-axis  distance  of  the  hydrophone  from  point  C  in  Fig.  2  is  z  and  ro  =  z  +  a.  Figure 
5  compares  the  theoretical  prediction  from  Eq.  (6)  with  experimental  data  for  ka  =  24.3. 
Inspection  of  Eqs.  (5b)  and  (6)  reveals  that  I  />m/(^0»Y)l/'Pswc(^0)l  depends  only  on  7  (note 
ro/r  =  [1  +  and  that  it  contains  no  adjustable  parameters.  The  agreement 

between  theory  and  experiment  is  excellent  for  a  measurement  of  this  type.  A  detailed 
discussion  of  the  theory  and  experimental  methods  used  in  obtaining  Figs.  4  and  5  is 
contained  in  Ref.  5  as  well  as  other  theoretical  and  experimental  results.  The  importance  of 
these  results  is  that  a  qualitative  ray  model  of  the  backscattered  pressure  described  the 
experimental  observations. 

The  remainder  of  this  dissertation  gives  a  more  detailed  analysis  of  the  forward  and 
backwards  scattering  of  a  plane  wave  from  an  evacuated  spherical  shell.  Chapter  2 
develops  quantitative  ray  methods  used  in  synthesizing  the  forward  scattering  amplitude 
and,  ultimately,  the  total  scattering  cross  section  via  the  optical  theorem.  Chapters  3  and  4 
investigate  quantitative  ray  methods  for  the  modeling  of  the  form  function  for 
backscattering.  The  contents  of  chapters  2  and  3  have  been  accepted  for  publication  in  the 
Journal  of  the  Acoustical  Society  of  America.  The  present  form  of  chapter  4  has  been 
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chapters  are  independent  papers,  they  may  be  read  separately.  As  such,  each  chapter 
contains  its  own  abstract,  introduction,  sections,  appendices,  and  references.  The 
numbering  of  equations,  figures,  and  tables  are  unique  to  a  given  chapter.  Mathematical 
symbols  may  differ  slightly  between  chapters,  however,  each  symbol  is  clearly  defined 
within  the  appropriate  chapter.  It  is  suggested  that  chapter  3  be  read  prior  to  chapter  4, 
since  some  results  from  chapter  3  are  used  in  chapter  4.  Two  appendices  follow  the  main 
text.  Each  appendix  contains  material  germane  to  the  ray  synthesis  of  the  forward  and 
backwards  scattering  from  an  evacuated  elastic  spherical  shell.  This  material  could  not  be 
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adequately  presented  within  chapters  2  -  4,  but  it  is  deemed  important  for  completeness  of 
the  analysis.  Appendix  I  discusses  Marston's  approximation  for  the  complex  coupling 
coefficient  and  comparison  to  an  apparent  exact  complex  coupling  coefficient  for  each  leaky 
Lamb  wave  is  presented.  The  dependence  of  the  leaky  Lamb  wave  resonance  spacing  on 
the  phase  and  group  velocities  is  also  discussed  in  Appendix  I,  Finally,  Appendix  II 
contains  tables  of  the  leaky  Lamb  wave  parameters  necessary  for  the  ray  synthesis.  These 
parametei’S  were  generated  by  extending  certain  results  of  the  Sommerfeld-Watson 
transformation  of  the  partial  wave  series  for  backscattering  from  a  solid  elastic  sphere  to  the 
shell. 
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Chapter  2 

Ray  synthesis  of  Lamb  wave  contributions  to  the  total  scattering  cross 
section  for  an  elastic  spherical  shell 


Abstract 

The  optical  theorem  relates  the  extinction  cross  section,  Ceika),  to  the  forward 
scattering  amplitude, /(0  =  0,  ka).  Here  0  denotes  the  scattering  angle,  k  is  the 
wavenumber  of  the  incident  sound,  and  a  is  the  radius  of  the  scatterer.  If  the  absorption  by 
the  scatterer  is  negligible  so  that  the  scatterer  is  elastic,  is  equal  to  the  total  scattering 
cross  section  By  applying  this  theorem  to  the  partial  wave  series  for f(0,ka),  we  can 
obtain  an  expression  for  O/  for  an  elastic  spherical  shell  in  water.  However,  the  series 
representation  of  O/  does  not  facilitate  a  direct  understanding  of  the  rich  structure  caused  by 
the  shell's  elastic  response.  In  paiticular,  the  elastic  response  is  attributable  to  leaky  Lamb 
waves.  We  employ  a  generalization  of  the  geometrical  theory  of  diffraction  [P.  L. 

Mar.ston,  J.  Acoust.  wSoc.  Am.  25-37  (1988)]  to  synthesize /(O,  ka).  This  simple  ray 
acoustic  synthesis  contains  a  component  for  ordinary  diffraction  by  the  shell  and  distinct 
contributions  for  the  individual  Lamb  waves  that  can  be  excited  on  the  shell.  A  comparison 
of  numerical  computations  for  O/  utilizing  the  exact  partial  wave  series  and  the  ray 
synthesis  shows  good  agreement  in  the  description  of  the  resonance  structure.  The  relevant 
range  of  ka  for  this  comparison  is  7  <  ka  <  100.  The  elastic  material  of  the  shell  is  440c 
stainless  steel  and  the  inner-to-outer  radius  ratio  is  b/a  =  0.838.  Dispersion  curves  and 
radiation  damping  for  Lamb  waves  were  calculated  by  Watson  transform  methods.  The 
structure  in  ©/(ka)  due  to  Lamb  waves  may  also  be  depicted  as  a  manifestation  of  forward 
glory  scattering  and  experimental  evidence  for  the  forward  glory  is  noted. 
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2.1  Introduction 

Several  articles  have  concentrated  on  the  surface  elastic  wave  contribu  tions  to 
backscattering  from  solid  elastic  spheres^*  2. 3  and  cylinders'^  as  well  as  spherical^’  ^  and 
circular  cylindrical  shells.^’  ^  The  surface  elastic  waves  on  solid  spheres  and  cylinders  are 
leaky  Rayleigh  and  whispering  gallery  waves;  while  the  surface  elastic  waves  on  spherical 
and  cylindrical  shells  are  a  generalization  of  Lamb  waves.  (Although  the  strains  due  to 
Lamb  waves  are  not  strictly  confined  to  the  surface  of  the  shell,  it  may  be  appropriate  to 
label  l.amb  waves  as  surface  elastic  waves,  since  the  surface  of  the  shell  guides  the  Lamb 
wave  propagation.)  In  particular,  some  of  these  authors  have  considered  the  significance 
of  surface  elastic  wave  resonances  to  the  backscattering  form  function, /(0  =  jc)  defined  in 
Eq.  (1)  below.  Since  the  structure  found  in  the  backscattering  form  function  is  attributed  to 
the  presence  of  surface  elastic  wave  resonances,  then  these  resonances  should  also 
contribute  to  the  scattering  in  the  forward  direction.  In  this  article,  we  are  concerned  with 
the  manifestation  of  such  contributions  to  the  total  scattering  cross  section  for  an  elastic 
spherical  shell. 

When  a  scatterer  is  placed  in  an  incident  acoustic  plane  wave  of  intensity 
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section.^^  The  extinction  cross  section,  denoted  as  Og,  expresses  the  total  power  removed 
(scattered  or  absorbed)  from  the  incident  plane  wave  as  lOg,  From  the  conservation  of 
energy  Og  =  o,  +  Gabs  where  I  Gabs  is  the  power  absorbed  by  the  scatterer  and  Gabs  is  known 
as  the  absorption  cross  section.  The  extinction  cross  section  is  related  to  tixc  forward 
scattering  amplitude  by  the  optical  theorem  reviewed  in  Sec.  2,2  below,  iLl2  (This 
theorem  has  also  been  referred  to  as  the  extinction  theorem.)  Some  early  examples  of  the 
application  of  the  optical  thexirem  to  acoustical  problems  are  given  in  Refs.  13-15.  For  a 
scatterer  with  no  absorption.  Gats  =  0,  so  that  Gg  =  g^  and  the  optical  theorem  may  be  used 
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to  calculate  O/  directly.  (Even  for  nonabsorbing  scatterers,  however,  some  authorsi6-l8 
resorted  to  integration  of  differential  cross  section  over  471  sr  to  obtain  0[.)  The  scatterer 
considered  in  the  present  paper  is  modeled  as  perfectly  elastic. 

In  this  paper,  the  forward  scattering  amplitude  is  obtained  from  the  exact  partial 
wave  series  and  from  a  synthesis  in  terms  of  ray  contributions.  The  synthesis  is  partially 
facilitated  by  the  elastic  generalization  of  the  geometrical  theory  of  diffraction  given  in  Ref. 
9. 

This  paper  is  organized  as  follows.  In  Sec.  2.2,  the  optical  theorem  is  stated  for  the 
case  of  scattering  from  spheres.  Using  the  optical  theorem  and  partial  wave  series 
representation  for  the  forward  scattering  amplitude,  o,  for  an  elastic  spherical  shell  is 
obtained.  The  basis  of  the  ray  synthesis  of  O/  is  outlined  at  the  end  of  Sec.  2.2.  In  Sec. 
2.3,  the  component  of  Ci  due  to  ordinary  forward  diffraction  is  developed  for  the  ray 
model.  Individual  Lamb  wave  contributions  to  the  ray  representation  of  the  forward 
scattering  amplitude  are  determined  in  Sec.  2.4.  Section  2.5  compares  the  results  of  the 
exact  partial-wave  series  calculation  with  the  ray  model.  The  pertinent  conclusions  from 
the  comparison  in  Sec.  2.5  are  contained  in  Sec.  2.6  along  with  a  discussion  of  relevant 
aspects  of  forward  glory  scattering  for  sphere.s.  Section  2.6  also  surnmarizes  an 
experimental  observation  pertinent  to  the  forward  glory  of  shells  and  comments  on  the 
cause  of  certain  frequency  dependent  stmcture  in  a/.  Appendix  A  discusses  the  Lamb 
wave  parameters  used  in  the  synthesis  as  given  by  Watson  transform  methodology. L3 ,5,9 

2.2  The  optical  theorem  for  acoustic  scattering 
A.  Partial  wave  series  analysis 

A  unit  amplitude  plane  wave,  propagating  in  the  positive  z-direction,  in  water  is 
incident  upon  an  elastic  spherical  shell.  The  total  pressure  in  the  water  can  be  expressed  us 
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a  superposition  of  the  incident  plane  wave  and  a  scattered  wave.  In  particular,  in  the 
farfield  of  the  shell,  the  total  pressure  is  ^ 

Pj  =  +  (a/2r)f(0)e‘‘^.  ( i ) 

The  first  term  is  the  incident  plane  wave  with  wavenumber  k  =  co/c.  The  second  term  is  a 
spherically  diverging  scattered  wave.  The  outer  radius  of  the  shell  is  denoted  by  a  and/(6) 
is  the  complex  dimensionless  scattering  amplitude.  The  harmonic  time  dependence,  expi- 
zo)f),  has  been  and,  henceforth,  will  be  suppressed. 

The  optical  theorem  provides  a  fundamental  relationship  between  the  scattering 
amplitude  and  a^.  Several  authors  give  various  expressions  for  the  optical  theorem  (see, 
e.g..  Refs.  11,  12,  19,  and  20).  For  the  amplitude  normalization  used  in  Eq.  (1),  it 
follows  that 

Oe  =  (2rta/k)Im{f(e=0)},  (2) 

where  Im{ }  means  the  imaginary  part  of  the  enclosed  quantity.  The  significance  of  Eq.  (2) 
is  that  the  only  scattering  amplitude  needed  is  the  one  for  forward  scattering.  It  is 
convenient  to  normalize  Og  by  the  geometrical  cross  section  of  the  shell,  -lufi, 

(Cg/7ta^)  =  (2/x)  Im  { f(0  =  0) } .  (3) 

This  normalization  intrcxiuces  the  dimensionless  size  parameter,  x  =  ka,  into  Eq.  (3).  In 
the  discussion  which  follows  we  consider  perfectly  elastic  scatters  so  that  =  a^. 

The  partial  wave  series  (PWS)  representation  of  the  exact  forward  scattering 
amplitL  ^e  for  an  evacuated  elastic  spherical  shell  is 

f(0  =  O)  =  (2/ix)y  (2n  +  l)-^, 


(4) 
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where  n  is  the  partial  wave  index.  The  functions  Bn  and  Dn  are  5x5  detemunants  whose 
elements  are  complicated  combinations  of  spherical  Bessel  functions  of  the  first  and  second 
kind.5  Using  Eqs.  (3)  and  (4),  the  normalized  total  scattering  cross  section  for  the  shell 
becomes 

(a/7ia^)  =  (-4/x^)  ^  (2n  +  1)  Rej  D  (x)  ]  ’  (5) 

where  Re[ }  denotes  the  real  part  of  the  enclosed  quantity.  The  dashed  curve  in  Fig.  1,  that 
contains  the  complicated  structure,  is  the  normalized  O;  for  a  440c  stainless  steel  shell  as 
computed  from  Eq.  (5)  witli  the  material  parameters  given  in  Sec.  2.5.  There  are  two 
imponant  features  in  Fig.  1.  First,  <5t  is  approaching  the  expected  value  of  twice  the 
geometrical  cross  section  of  the  shell  as  x  This  asymptotic  result  for  <5t  is  known  as 
the  extinction  paradox.  Physical  interpretations  of  the  extinction  paradox  have  been 
discussed  by  various  authors.^®"^^  Second,  and  perhaps  more  important,  the  structure 
observexi  in  CTj  for  a  440c  stainless  steel  shell  is  a  manifestation  of  the  elastic  nature  of  the 
scatterer.  That  is,  a  comparison  of  C/  for  an  elastic  shell  and  a  rigid  sphere  of  radius  a  in 
Fig.  1  shows  that  the  elastic  properties  of  the  shell  are  important. 

Tlie  expression  for  the  normalized  cross  section  in  Eq.  (5)  has  properties  which 
would  make  a  geometrical  representation  worthy  of  investigation.  First,  a  simple  physical 
interpretation  of  tlie  structure  in  <5i  is  not  easily  obtained  directly  from  the  PWS 
representation.  Although,  it  is  known  that  the  incident  plane  wave  can  excite  Lamb  wave 
resonances  in  a  shell,  the  quantitative  description  for  this  coupling  is  not  immediately 
apparent  from  the  PWS  expression.  Second,  it  is  well-known,  that  for  large  to,  the  PWS 
is  a  slowly  converging  series;  the  number  of  terms  required  exceeds  to.  See  Eq.  (16) 
below.  We  examine  below  a  simple  expression  for  a<  based  on  a  generalization^  of  the 
geometiical  tlieory  of  diffraction23.  It  should  be  possible  to  modify  this  ray  representation 


Fig.  1  Normalized  total  scattering  cross  sections  for  empty  elastic  spherical  shell  in  water 
and  a  rigid  sphere,  each  of  outer  radius  a.  The  short  dashed  line  is  the  exact  partial  wave 
series  result  for  the  shell,  Eq.  (5),  with  the  stainless  steel  parameters  listed  in  Sec.  2.5  and 
an  inner  to  outer  radius  bja  =  0.838.  The  solid,  monotonic  curve  corresponds  to  the  rigid 
sphere,  where  the  panial  wave  series,  Eq.  (9),  is  calculated.  The  long  dashed  line  is  the 
first  three  terms  of  Beckmann  and  Franz's  approximation,  Eq.  (10),  for  the  normalized 
cross  section  for  the  rigid  sphere 
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of  the  total  scattering  cross  section  for  objects  of  other  shapes  where  the  cross  section  is 
not  easily  obtained  from  a  PWS. 

B.  Ray  synthesis  of  the  total  scattering  total  scattering  cross  section 

It  has  been  proposed,  that  the  form  function  describing  the  steady -state  high- 
frequency  backscattering  amplitude  from  spheres  can  be  approximated  by  three  ordinarily 
distinct  classes  of  contributions.^*^’^  These  contributions  arc;  a  specular  reflection, 
transmitted  bulk  wave  contributions,  and  surface  elastic  wave  contributions.  A 
superposition  of  these  components  gives  the  scattered  pressure  term  of  the  total  pressure  in 
the  backward  diiection.  (If  ka  is  smaller  than  about  15,  a  contribution  to  backscattering 
from  Franz-type  creeping  waves  can  also  become  significant  for  some  materials  or  shell 
thicknesses  but  such  contributions  were  not  separately  analyzed.)  By  superposing  the 
appropriate  contributions,  the  forward  scattering  amplitude  may  be  synthesized  as 
discussed  below. 

To  construct  the  geometrical  synthesis  of  the  forward  scattering  amplitude  for  the 
spherical  shell,  the  specular  reflection  term  for  backscattering  is  replaced  by  a  contribution 
that  describes  the  ordinary  forward  diffraction  about  the  shell  (see  Sec.  2.3  below).  The 
transmitted  bulk  wave  contributions,  that  backscattering  foim  functions  may  contain,  do 
not  occur  in  the  forward  scattering  amplitude  of  a  thin  elastic  air-filled  or  evacuated  shell. 
That  is,  there  is  no  mechanism  for  the  transmission  of  a  bulk  wave  through  the  interior  of 
the  shell.  Since  surface  elastic  waves  leak  energy  continuously  back  into  the  surrounding 
water,  then  these  contributions  must  also  be  included  in  the  forward  scattering  amplitude. 
Hence,  the  forward  scattering  amplitude  has  the  approximate  form 

f(8=0).fp^  +  £f,, 

/ 


(6) 
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where  the  summation  is  over  all  possible  surface  elastic  wave  contributions.  Since  the 
operator,  /m{ ),  is  linear,  then  by  inspection  of  Eqs.  (2)  and  (6)  the  total  scattering  cross 
section  is 


The  forms  of  the  individual  components  in  Eqs.  (6)  and  (7)  are  considered  in  Secs.  2.3  and 
2.4. 


2.3  Forward  diffraction  term  for  the  total  scattering  cross  section 

The  specific  fonn  of  the  contribution  to  due  to  ordinary  forward  diffraction  about 
an  elastic  spherical  shell  is  not  currently  available  from  a  Watson  transformation,^  The 
simplest  approximation  would  be  the  total  scattering  cross  section  of  a  rigid  disk  of  radius 
a.  That  is,  the  total  scattering  cross  section  for  a  rigid  disk  (computed  from  the  optical 
theorem)  can  only  be  attributed  to  ordinary  forward  diffraction.  The  asymptotic  value  (i.e., 
Ica»\)  for  the  normalized  total  scattering  cross  section  of  the  rigid  disk24  is  2.  As  noted 
previously,  (C;  /  na^)  2  as  /kj  oo  is  known  as  the  extinction  paradox.  Although 
forward  diffraction  from  a  rigid  disk  gives  the  expected  asymptotic  value,  this 
approximation  is  inappropriate  for  the  ka  range  of  interest,  7  <  ka  <  lOO. 

A  more  suitable  approximation  for  the  normalized  ordinary  forward  diffraction 
contribution  (Opiy  is  the  total  scattering  cross  section  for  a  fixed  rigid  or  acoustically 
hard  sphere.  To  obtain  (Opiy  na^),  we  employ  the  optical  theorem,  as  given  by  Eq.  (3), 
and  the  forward  scattering  amplitude  for  the  rigid  sphere.  The  forward  scattering  amplitude 
for  the  rigid  sphere  is  given  by  the  PWS^’ 

f^  =  (-2/ix)X(2n-M)-^, 

n=0  h„\x) 


(8) 


where  and  are  spherical  Bessel  functions  of  the  first  and  third  kind,  respectively.  The 
prime  in  Eq,  (8)  denotes  differentiation  with  respeet  to  the  argument  of  the  function.  The 
result  of  inserting  Eq.  (8)  into  Eq.  (3)  is 
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(apj^Tta^)  =  (4/x2)  ^  (2n  +  1)  Rej  ‘  (i)' 


n  =  0 


h„  (X) 


(9) 


The  solid,  monotonicaUy  increasing  curve  in  Fig.  1  is  ofd/JW^  for  a  rigid  sphere  of  radius 
a.  From  Fig.  1  one  observes  that  the  cross  section  for  the  rigid  sphere  asymptotically 
approaches  2.  Although  Eq.  (9)  is  a  suitable  candidate  for  (Opjy  its  dependence  on  x 
is  complicated  and  ,  hence,  not  a  simple  parameterization  of  the  ordinary  forward 
diffraction  contribution  to  O/  of  an  elastic  sphere. 

Beckmann  and  Franz  performed  a  modified  Watson  transformation  on  a  rigid 
sphere  result  similar  to  Eq.  (9).^^’  ^7  The  Watson  transformation  transforms  the  slowly 
convetging  PWS  into  a  more  rapidly  converging  series.  The  result  of  their  calculation  is  a 
series  in  inverse  powers  of  x, 

(apjj/7ta^)  =  2  -  1.7284  -  2.0104  +  0(x'^) .  (10) 


The  first  three  terms  in  Eq.  (10)  are  plotted  as  the  long-dashed  line  in  Fig.  1.  From 
inspection  of  Fig.  1,  we  observe  tliat  the  Beckmann  and  Franz  result  is,  perhaps,  the 
simplest  approximation  of  die  ordinary  forward  diffraction  term  in  Eq.  (7)  throughout  the 
region  of  interest  1  <x<  100.  Thus,  for  the  numerical  computations  given  in  Sec.  2.5,  we 
use  the  fi.rst  three  tenns  in  Eq,  (10)  for  the  contribution  of  ordinary  forward  diffraction  to 
a,. 

It  may  be  argued  that  the  second  term  in  Eq.  (10)  represents  approximately  thex 
dependence  of  contributions  of  grazing  rays  (also  known  as  edge  rays)  to  IrrJfi^  =  0)]  for 
a  rigid  spheie.  In  addition,  there  could  be  distinct  ray  contributions  which  wrap  completely 


around  a  rigid  sphere.  The  waves  associated  with  those  rays  are  creeping  waves  having 
phase  velocities  close  to  the  speed  in  the  outer  media  (water).  The  damping  of  such  waves 
is  sufficiently  large  that  contributions  from  a  complete  circumnavigation  to  forward 
scattering  are  generally  thought  to  be  negligible  for  impenetrable  spheres. 

2.4  Leaky  surface  wave  contributions  to  the  forward  scattering  amplitude 

In  this  section,  we  use  the  methodology  of  geometrical  theory  of  diffraction 
generalized  to  elastic  objects^  to  derive  an  expression  for  an  individual  surface  elastic  v/ave 
contribution,//,  to  the  forward  scattering  amplimdc.  The  acoustic  ray  diagram  in  Fig.  2 
facilitates  an  understanding  of  this  ray  synthesis  of //.  The  ray  diagram  is  similar  to  those 
discussed  in  Refs.  1,  5  and  9.  An  acoustic  plane  wave  propagates  along  tlie  ray  AB  in  the 
positive  z-direction.  At  the  point  B  the  acoustic  wavefield  in  water  couples  to  the  elastic 
shell  and  launches  a  surface  elastic  v'ave.  The  surface  elastic  wave  propagates  from  B  to 
B'  shedding  energy  continuously  back  into  the  water.  At  B’  the  radiated  energy  is  directed 
along  the  forward  scattering  direction  (parallel  to  the  +z-axis).  The  surface  elastic  wave 
repeatedly  circumnavigates  the  shell  radiating  energy  each  time  in  the  forward  direction. 
The  points  B  and  are  determined  by  the  phase  velocity  trace- matching  condition, 

0^  =  arcsin(c/cp,  >  c,  ( 1 1 ) 

where  0/  is  the  local  angle  of  incidence  and  c/c/  is  the  ratio  of  the  speed  of  sound  in  water, 
c,  and  the  phase  velocity  along  the  shell's  outer  surface  ol  the  Lamb  v/ave,  c/.  The  phase 
velocity  is  assumed  to  be  supersonic. 

The  form  of  the  individual  surface  elastic  wave  contribution  to  the  forward 


scattering  amplitude  is 
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Fig.  2  Acoustic  ray  diagram  for  Lamb  wave  contributions  to  the  forward  scattering 
amplitude  for  Lamb  waves  having  c/  >  c.  The  outer  radius  of  the  shell  is  a  and  the  inner 
radius  is  b.  The  incident  acoustic  plane  wave  launches  a  Lamb  wave  in  the  vicinity  of  point 
B.  The  Lamb  wave  propagates  along  the  shell  and  radiates  in  the  forward  direction  at  point 
B'.  The  points  B  and  B'  are  determined  by  Eq.  (11).  The  qualitative  features  of  the  ray 
diagram  are  made  quantitative  through  Eqs.  (12)  -  (14)  to  describe  the  scattering.  The  ray 
synthesis  also  contains  a  contribution  due  to  diffraction  about  the  shell.  That  includes 
contributions  due  to  rays  which  touch  the  shell  near  E  at  the  edge. 
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f  -  -  r,  V  J2mnxc/c^ 

^r  G/®  ^  ^  ^  (12) 

m  =  0 


The  coupling  efficiency  of  the  acoustic  wavefield  in  water  with  the  elastic  shell  at  the  points 
B  and  B'  in  Fig.  2  is  characterized  by  a  complex  coupling  coefficient  Gi.  The  geometrical 
phase  shift, 


-  x(c/cp(TC  -  20p  -  2x  cos(0p 


(13) 


accounts  for  the  difference  in  phase  of  a  surface  elastic  wave  traveling  from  BioB'  along 
the  shell  and  a  plane  wave  in  water  traveling  from  B  to  B'  as  if  the  shell  were  not  present. 
The  factor  exp\r  (ji  -  20/  )P/]  is  associated  with  the  radiation  damping  of  the  surface  elastic 
wave  propagating  along  the  arc  BB'.  The  radiation  damping  parameter  is  denoted  by  P/ 
and  has  the  units  of  np/radian.  The  summation  occurs  in  Eq.  (12)  from  the 
circumnavigations  of  the  surface  elastic  wave.  The  m  =  0  term  is  the  first  partial 
circumnavigation  of  the  shell;  while  m  - 1, 2,...  is  for  each  successive  complete 
circumnavigation.  The  terms  exp(-2mn^i)  and  expipLrmvcc/ci)  are  the  additional  radiation 
damping  and  phase  shift.  Finally,  the  term  expiimn)  arises  from  the  passage  of  the  surface 
elastic  wave  through  the  caustics  at  the  points  C'  and  C". 

It  is  observed  that  the  summation  in  Eq.  (12)  is  a  geometric  series  and  may  be 
expressed  in  a  simple  analytic  form.^  Summing  the  series  in  Eq.  (12),  we  have 
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TTie  form  of  the  denominator  in//  is  that  of  a  Fabry-Perot  resonator.3.9  Equation  (14)  is 
the  desired  ray  synthesis  for //, 

The  coupling  coefficient,  G/,  for  the  elastic  solid  sphere  was  determined  from  the 
Sommerfeld-Watson  transfonnation  on  the  appropriate  PWS.l’  ^  At  the  time  of  this 
writing,  an  expression  for  the  coupling  coefficient  for  the  elastic  spherical  shell  is  not 
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available  directly  from  the  Watson  transformation.  However,  Marston  has  derived 
approximations  for  Gi  for  an  elastic  sphere  and  cylinder.9.29  xhe  approximations  were 
supported  by  various  numerical  results  for  solid  spheres  and  cylindrical  shells.^  The 
relevant  result  for  spheres  may  be  written 

=  87tp^  (c/cp  e'*^'  (15) 

where  <])/  denotes  the  phase  of  G/  and  it  may  be  argued^^  that  ~  0.  For  the  calculation  in 
Sec.  2.5,  <{>/  for  the  Ith  surface  elastic  wave  is  taken  to  be  identically  equal  to  zero.  It  is 
noteworthy  that  the  magnitude  of  G/  depends  on  only  the  damping  parameter  and  the  ratio 
of  the  velocities.  This  result  for  I  G/ !  was  confirmed  in  Ref.  5  for  the  stainless  steel  shell 
considered  here  by  use  of  measured  backscattering  amplitudes.  The  computations 
discussed  in  Sec.  2.5  give  additional  support  for  Eq.  (15)  with  (])/  =  0  as  does  a  synthesis 
of  the  steady  state  backscattering  amplitude  for  the  shell  considered.39 

The  ray  synthesis  of  the  contribution  of  a  particular  surface  elastic  wave  is  given  by 
Eqs.  (11),  (13)  and  (14)  with  the  approximation,  Eq.  (15),  used  for  the  coupling 
coefficient.  Inspection  of  Eqs.  (1 1)  -  (15)  shows  that  a  particular  surface  elastic  wave 
contribution  depends  on  only  the  damping  parameter  and  the  phase  velocity  ratio.  The 
method  used  for  determining  p/  and  c/c/  is  described  in  Appendix  A. 

2.5  Numerical  results  and  discussion 

The  calculations  shown  in  Figs.  3  -  8  are  for  a  440c  stainless  steel  shell  immersed 
in  water.  For  these  computations,  the  material  parameters  for  440c  stainless  steel  were 
selected  because  of  previous  work  with  a  shell  composed  of  this  material.^  The 
longitudinal  velocity,  c^,  for  440c  stainless  steel  is  5.854  mm/\xs  and  the  shear  wave 
velocity,  Cg,  is  3.150  mm/iXs.  The  density  is  ps  =  7.84  g/crn^.  The  shell's  outer  radius  is 


30 


a  =  19.05  mm  and  the  ratio  of  the  inner-to-outer  radius  is  b/a  =  0.838.  The  speed  of  sound 
for  the  water  is  c  =  1.479  mml\is\  while  the  density  is  p  =  1.00  glcm^- 

For  the  PWS  computations,  we  use  Eq.  (5),  where  the  determinants  and  are 
given  in  Ref.  5.  The  PWS  calculation  for  (O/  /luP-)  is  perfonned  over  the  range  0  <  a:  < 

100.  In  Figs.  3-8,  the  dashed  curves  are  the  exact  PWS  results.  When  summing  PWS, 
the  following  maximum  values  for  the  partial  wave  index  were  tested  and  found  to  ensure 
convergence: 

n  =  2  +  [  X  +  4.0  x^^  ],  X  <  8, 

max  .. 

1/3  (1^) 

n  =  3  +  [  X  +  4.05  X  '  ],  x  >  8, 

moA 

where  the  square  brackets  imply  the  integer  part  (i.e.,  rounded  up  or  down)  is  to  be  used. 
These  values  are  similar  to  ones  used  when  evaluating  the  Mie  PWS  for  light  scattering.^  ^ 
That  the  maximum  partial  wave  index  in  the  optical  and  acoustical  cases  should  be  similar 
may  be  argued  from  the  localization  principle.  (The  proper  convergence  of  the  PWS  for 
the  shell  considered  was  verified  by  substantially  increasing  nmax-  Equation  (16)  should  be 
used  with  caution,  however,  for  thin  shells  having  significant  subsonic  wave 
contributions.)  Since  a  Lamb  wave  resonance  can  contain  sharp  peaks,  the  increment  of 
needs  to  be  small.  Hence,  it  is  apparent  that  the  implementation  of  Eq.  (5)  can  be  relatively 
computationally  intensive. 

The  particular  Lamb  waves,  that  can  be  coupled  to  a  shell  with  the  above  material 
properties  and  in  the  range  1  <x<  1(X),  are  the  two  lowest  antisymmetric  or  flexural  modes 
Oq  and  G]  and  the  two  lowest  symmetric  or  dilatational  modes  aq  and  Aj.  The  present  ray 
synthesis  is  restricted  to  j:  >  7  for  two  reason.  First,  the  numerical  procedure  used  to 
obtain  Lamb  wave  damping  parameters  and  phase  velocities  begins  to  breakdown  as  noted 
in  Appendix  A.  That  is,  the  numerical  accuracy  of  the  results  become  suspect.  Second,  the 
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oq  mode  becomes  subsonic  for  x<l.  The  present  ray  synthesis  is  limited  to  supersonic 
surface  elastic  waves  as  evident  from  Eq.  (11)  and  the  related  analysis  given  in  Refs.  1  and 
3. 

The  solid  curves  in  Figs.  3  -  8  are  the  ray  synthesis,  where  {a^Tia^)  is  given  by 
the  first  three  terms  in  Eq.(lO)  (a  monotonic  function  of  jc)  and  the  appropriate  Lamb  wave 
contributions  are  included.  In  Fig.  3,  the  ray  result  has  only  two  terms.  These  terms  are 
the  ordinary  forward  diffraction  component  and  the  lowest  antisymmetric  Lamb  mode 
contribution,  It  is  apparent  that  the  course  structure,  for  1  <x<  25,  is  attributed 

entirely  to  the  ao  fle.'rural  Lamb  wave.  The  sharp  resonances  in  Fig.  3,  that  are  not 
modeled,  are  from  the  contribution.  Figure  4  includes  The  agreement  of  the  ray 
synthe.sis  and  the  exact  PWS  for  (<t,  is  good  considering  the  ray  model  is  based  on 
high-frequency  methexis.  Observe  that  the  above  synthesis  appears  to  be  accurate  down  to 
jc  *  9.  It  is  noteworthy,  however,  that  the  .vq  resonances  near  ka  of  17.60,  21.28,  and 
24.95  may  not  have  been  fully  resolved  due  to  the  narrowness  of  their  line  widths. 

The  region  25  <  a:  <  50  is  plotted  in  Figs.  5  and  6.  The  vertical  axis  has  been 
magnify  to  enhance  the  visibility  of  narrow  resonance  features  and  to  better  illustrate  the 
agreement  between  the  ray  and  PWS  results.  In  Fig.  5,  the  ray  i-odel  contains  the  ordinary 
forward  diffraction  term  and  the  Oq  and  sq  Lamb  wave  contributions.  Again,  the  broad 
structure  is  associated  with  the  «q  Lamb  v/ave  while  the  sharp  resonances  are  attributable  to 
the  Sq  Lanrib  wave.  However,  as  noted  in  Appendix  A,  the  Cj  antisymmetric  mcjde  can  be 
excited  on  the  shell  for  >  40.  In  particular,  the  resonance  features  atx~  43, 44, 46  and 
50  are  not  associated  with  either  the  uq  or  sq  modes.  By  adding  the  aj  contribution  in  this 
region,  we  recover  all  the  resonance  features  of  the  PWS  result  in  the  ray  synthesis  of  the 
cross  section. 
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Fig.  4  A  comparison  of  the  same  ka  range  as  in  Fig.  3,  but  the  CTsq  contribution  has  now 
been  included  in  the  ray  synthesis.  The  broad  and  narrow  structures  appear  to  be  correctly 
synthesized.  The  resolution  of  each  plot  was  enhanced  relative  to  the  plots  in  Figs.  2  and  3 
so  as  to  more  accurately  describe  the  narrow  sq  lamb  wave  resonances. 


Fig.  5  The  clashed  curve  is  the  partial  wave  series  result  for  the  normalized  total 
scattering  cross  section.  The  ray  synthesis  (solid  curve)  contains  the  contributions  Opj), 
Oao  and  <73^.  The  broad  structure  is  associated  with  the  oq  Lamb  wave;  while  the  narrow 
pealcs  and  dips  are  from  the  jq  Lamb  wave.  This  level  of  synthesis  fads  to  reproduce 
features  at  ka  ~  43,  44,  46,  48,  and  50. 
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The  ray  synthesis  for  50  <  x  <  75  is  shown  in  Fig.  7.  The  exact  PWS  computation 
is  the  short  dashed  line.  The  solid  line  is  a  ray  synthesis  of  the  normalized  a/;  where  the 
relevant  contributions  are:  (cpD/Tta^),  ao,  sq,  and  ai.  A  comparison  of  the  ray  and  PWo 
results  demonstrates  that  these  contributions  are  sufficient  to  describe  both  the  broad  and 
narrow  resonances.  In  particular,  the  features  in  the  normalized  O/  forx  <  71  arc  entirely 
due  to  the  ao,  sq,  and  ai  Lamb  waves.  For  jc  >  7 1,  the  acoustic  wavefield  in  water  can 
couple  into  the  Lamb  wave.  The  significance  of  the  Lamb  wave  to  the  ray  synthesis 
is  observed  from  the  long  dashed  line  in  Fig.  7.  By  adding  the  Lamb  wave,  the  ray 
model  is  found  to  give  a  better  approximation  of  the  PWS  calculation  for  50  <  x  <  75. 

The  last  figure.  Fig.  8,  corresponds  to  the  range  15<x<  100.  The  ray  model 
contains  the  (opD/ita^),  ao.  ■^0.  and  .vi  contributions.  An  examination  of  the  individual 
components  in  the  synthesis  allows  one  to  identify  the  particular  surface  wave  responsible 
for  a  specific  structure.  It  is  important  to  notice  the  synthesis  is  essentially  the  same  as  the 
PWS  computation  except  for  a  small  offset.  That  is,  the  ray  synthesis  accurately 
reproduces  the  resonance  features  in  this  region.-  Finally,  the  ray  synthesis  appears  to  be 
approaching  the  extinction-paradox  value  of  2. 

Tne  offset  of  the  ray  result  from  the  PWS  computation,  evident  in  Figs.  5-8, 
does  not  appear  to  be  related  to  the  ray  model  of  the  Lamb  wave  contributions  to  since 
the  shape  of  the  modeled  structure  agrees  with  the  PWS  computation.  The  offset  appears 
instead  to  be  related  to  the  use  of  a  rigid  sphere  model  of  ordinary  forward  diffraction. 
Numerical  tests  show  that  the  approximation  for  the  total  scattering  cross  section  of  a  rigid 
sphere,  Eq.  (10),  converges  to  the  exact  rigid  PWS  result,  Eq.  (9).  For  example,  at;c  = 
75,  Eq.  (10)  yields  a  value  of  1.8965  for  apo/Tca^,  while  the  exact  result  from  Eq.  (9)  is 
1.8971.  The  difference  of  these  two  calculations  is  0.0006;  while  the  offset  atx  =  75  is 
0.0344.  Clearly,  the  use  of  Eq.  (10),  instead  of  Eq.  (9),  for  the  ordinary  forward 
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diffraction  does  not  account  for  the  offset.  Furthermore,  a  comparison  of  a  synthesis  (not 
shown),  employing  Eq.  (9),  with  the  exact  PWS,  Eq.  (5),  shows  a  similar  offset.  Hence, 
the  total  scattering  cross  section  for  the  rigid  sphere  is  not  entirely  appropriate  for  the 
description  of  the  ordinary  forward  diffraction  contribution  to  O/  of  an  elastic  shell.  The 
ordinary  forward  diffraction  component  is  also  not  accurately  modeled  by  the  total 
scattering  cross  section  of  a  soft  sphere.  That  cross  section  for  a  soft  sphere  approaches  2 
from  above  as  x  — >  «>  and  diverges  to  +«>  asx  It  appears  that  the  correct  model 

for  OpD  should  lie  between  the  values  of  die  cross  sections  for  rigid  and  soft  spheres.  This 
conjecture  may  also  be  supported  by  computations^  indicating  that  for  shells,  poles  in  the 
complex  ka  plane  "tentatively  associated  with  the  Franz  wave  are  considerably  removed 
from  their  counterparts  for  a  rigid  sphere,  in  the  low-ka  region." 

It  is  noteworthy  that  the  effect  of  including  the  5i  Lamb  wave  near  the  right-hand 
side  of  Fig.  7  is  to  reduce  an  even  greater  offset  between  curves  which  is  otherwise 
present.  Inspection  of  Fig.  A1  shows  that  in  this  region  |3/  for  /  =  5i  is  relatively  large 
which  suggests  that  rays  having  large  (3/  can  give  rise  to  a  smooth  shift  in  0[. 

'Fhe  structure  manifest  in  Fig.  4  in  the  region  x<7  merits  discussion.  It  is 
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phase  and  magnitude  of/ in  Eq.  (3)  can  vary  rapidly  with  x.  Without  attempting  to  model 
detailed  resonance  manifestations  it  is  possible  to  comment  on  the  nature  of  the  relevant 
modes.  It  is  shown  in  Eq.  (A6)  that  the  frequency  of  the  lowest  purely  radial  or 
"breathing"  model  occurs  at  x  ~  6.3.  Calculations  of  contributions  to /(0  =  7t)  given  by 
Sammelmann  et  al.,^  based  on  a  generalization  of  resonance  scattering  theory,  indicate, 
however,  that  the  breathing  resonance  contributes  weakly  to  scattering  unless  the  shell  is 
much  thinner  dian  the  one  considered  here.  As  evident  from  elementary  considerations  and 
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the  discussion  of  Eqs.  (A4)  -  (A7),  the  breathing  mcde  is  associated  with  die  I  =  so  Lamb 
wave. 

The  structure  in  tlie  region  x^7  is  evidently  associated  with  subsonic  flexural  (/  = 

oo)  Lamb  waves.  Note  that  the  peaks  in  the  plot  of  are  relatively  similai'  in 
appearance  (or  x<l  so  it  may  be  anticipated  that  the  underlying  mechanism  is  similar  for 


each  peak.  Furthemtore  it  is  physically  plausible  that  each  of  these  peaks  should  be 
associated  in  some  way  with  resonance  scattering,  especially  since  we  are  viewing  the  total 
scattering  in  a  region  of  x  where  Cpo/mP-  is  relatively  small.  What  may  be  surprising  is 
that  Gi  (and  hence  G^)  has  peaks  of  over  four  times  the  physical  cross  section,  na^,  of  the 
sphere.  The  maximum  having  the  lo west  x  occurs  at  x:- 2. 5.  The  elastic  mode  of  a 
spherical  shell  having  the  lowest  natural  frequency  is  ordinarily  associated  with  the  n  =  2 
(or  quadrupole)  partial  wave  with  I  =  ao  (or  flexural)  motion  of  the  shell.  (See,  e.g.,  Ref. 
6, 32,  and  33.)  In  conjunction  with  the  aforementioned  peak  magnitude  of  gJts^,  the 
computations  of  Skelton  and  Waterhouse^^  of  the  acoustic  energy  streamlines  in  water  near 
slightly  absorbing  steel  shells  are  noteworthy.  Those  calculations,  which  were  performed 
with  an  incident  wave  at  the  co  of  the  n  =  2  resonance,  suggest  that  the  energy  streamlines 


may  be  deflected  toward  the  SuCll  over  a  cross-sectional  area  rnucli  laigci  than  that  of  the 


shell.  (Computations  of  electromagnetic  energy  flow  for  light  incident  on  particles^  have 
similarly  been  useful  for  understanding  large  optical  Ge.)  'Fhe  other  peaks  in  the  region  x  < 


1  are  evidently  associated  with  resonances  with  «  >  2.  Spherical  shells  typically  have 
several  flexural  resonances  with  n  >  2  at  co  below  the  breathing  resonancc.^^ 


2.6  Forward  glory  scattering,  physical  interpretation,  and  conclusions 

The  ray  synthesis  of  the  total  scattering  cross  section  O/  for  a  44()c  stainless  steel 
shell  demonstrates  the  usefulness  and  power  of  this  method.  The  simple  parameterization 
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of  the  ray  model  correctly  gives  the  surface  elastic  wave  resonance  features  of  the  Lamb 
wave  contributions.  Furthermore,  this  parameterization  of  Ct  reduces  the  problem  to  only 
the  individual  Lamb  wave  damping  parameters,  (5/,  and  the  phase  velocity  ratios,  cict. 

Thus,  a  simple  physical  interpretation  of  the  observed  spiicture  in  Figs.  3  -  8  can  be 
achieved. 

These  computations  confirm  the  validity  of  the  numerical  methods  used  in  obtaining 
p/  and  c/ci  for  the  Ith  type  Lamb  wave.  The  above  calculations  give  further  support  for 
Marston's  approximation,  Eq.  (15),  for  the  coupling  coefficient.^>2^  Also,  Figs.  3  -  8 
demonstrate  the  correctness  of  the  assumed  phase  tp/  of  G/.  That  is,  if  (p/  differs 
significantly  from  zero,  then  the  contribution  of  the  Ith  Lamb  wave  to  O/  would  be  altered 
by  the  optical  theorem,  Eq.  (3).  For  example,  we  have  confirmed  that  setting  (p/  =  7t 
inverts  the  resonance  structure  in  the  synthesis  and  destroys  the  agreement  with  the  partial 
wave  series  result. 

The  ray  shed  in  the  forward  direction  from  point  B'  in  Fig.  2  appears  to  emanate 
from  within  the  sphere.  The  location  of  the  virtual  source  may  be  constructed  by 
considering  the  crossing  of  the  adjacent  dashed  ray  which  differs  infinitesimally  in 
direction.  By  generalization  of  a  discussion  given  in  llie  Appendix  of  Ref.  2  lo  the  prescui 
case  of  forward  scattering  it  may  be  shown  that  the  virtual  focus  F[  where  the  dashed  ray 
crosses  the  forward  ray  is  located  a  distance  equal  to  the  outer  radius  a  of  the  sphere  behind 
tile  vertical  line  through  C".  Since  Fig.  2  may  be  rotated  about  the  C'C"  axis,  the  virtual 
source  is  ring-like  and  the  forward  directed  wavefront  has  locally  the  shape  of  a  circular 
torus.  Such  outgoing  toroidal  wavefronts  from  spheres  were  previously  investigated  for 
acoustical  backscattering^'^.^S  and  give  rise  to  what  has  been  termed  "glory  scattering" 
because  the  enhancement  of  light  scattering  from  cloud  droplets  (known  as  the  optical 
glory)  is  attiibutable  to  toroidal  outgoing  optical  wavefronts. 
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From  the  discussion  above  it  is  evident  that  the  structure  in  Figs.  3-8  attributa,bie 
to  Lamb  wave  contributions  is  also  a  manifestation  of  forward-directed  glory  scattering. 

The  important  point  is  that  the  rays  associated  with  the  Lamb  v/aves  couple  on-to  and  off-of 
the  sphere  with  nonzero  impact  parameters  so  the  associated  outgoing  wavefront  is  toroidal 
and  the  contribution  to  the  form  function /(0  =  0)  is  enhanced  by  axial  focusing.  While  the 
backward  directed  glory  has  long  been  of  interest, Nussenzveig  and  Wiscombe^^ 
relatively  recently  noted  that  the  forward  optical  glory  of  dielectric  spheres  gives  rise  to  an 
oscillation  of  the  normalized  optical  extinction  cross  section  which  is  quasi-periodic  in  the 
optical  ka.  (The  mechanism  for  producing  optical  glory  rays  in  solid  dielectric  spheres  is 
somewhat  different  from  that  considered  here  for  elastic  shells.)  The  scattering  pattern 
associated  with  a  forward  directed  optical  glory  has  been  observed  in  polarized  light 
scattered,  from  bubbles  in  a  viscous  silicone  oil.37  Direct  experimental  evidence  of  forward 
glory  scattering  of  sound  from  shells  due  to  Lamb  waves  has  been  obtained  for  scattering 
of  short  tone  bursts  in  experiments  similar  to  tho.se  de.scrib»xi  in  Ref.  5  but  with  the 
hydrophone  placed  on  the  forward  axis.^^>^^  What  follows  is  a  brief  summary  of  an 
observation  relevant  to  the  present  discussion. 

Q  cHr>\l/C  Q  firn®  oirFr»r»l  n  Vi  r\wy 

^  *45^**^^  ^  LJ  V*  vkAliv  Wi  i»»V  I-lV/iil  M  piUWU.  V/li 

the  forward  axis.  The  sphere  is  440c  stainless  steel  of  radius  a  =  19.05  mm  and  bfa  =■ 
0.838.  The  incident  burst  was  four  cycles  of  a  653  kHz  sine  vave  which  corresponds  to 
ka  ~  53.  The  exjrerimental  technique  is  similar  to  the  one  described  in  Ref.  5  for 


backscattering  from  the  same  shell.  The  large  burst  is  attributable  to  rays  in  the  water 
which  graze  the  sphere.  (That  may  be  demonstrated  by  removing  the  sphere  and  obser/ing 
the  arrival  time  of  a  signal  diiectly  from  the  source.)  The  burst  which  arrives  1 1  fxs  prior  to 
the  large  burst  is  attributable  to  forward  glory  scattering  associated  with  Lamb  waves  on  the 
sphere.  It  arrives  prior'  to  the  large  burst  since  the  (slower)  grazing  ray  has  its  entire  path  in 
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water.  From  inspection  of  Fig.  2,  the  burst  which  travels  along  path  ABB' A '  is  advanced 
in  time  relative  to  a  burst  which  just  grazes  the  sphere  at  point  £  by  an  amount 

At^  =  (2a/c)  cos(0p  -  (a/c^pCtc  -  20^  , 

where  the  calculation  of  the  group  velocity  Cgi  for  the  Ith  Lamb  wave  is  discussed  in  Ref.  5. 
This  burst  is  associated  with  the  m  =  0  term  in  Eq.  (12).  It  is  noteworthy  that  for 
supersonic  Lamb  waves  {ci  >  c)  where  Fig.  2  is  applicable,  that  At/  >  0  when  (as  is  usually 
the  case)  Cgi  >  c.  For  the  sphere  under  consideration,  Eq.  (17)  predicts  At/  =  1 1  |iJ  for 
both  /  =  OQ  and  sq  when  ka  -  53,  in  agreement  with  the  observations.  From  the  arrival  time 
alone  it  is  not  possible  to  discriminate  between  /  =  ao  and  so  Lamb  waves  mechanisms  at 
this  ka.  (At  other  ka,  the  calculated  At/  for  these  /  differ.)  The  late  suucture  of  the  record 
shown  in  Fig.  9  is  at  least  partially  associated  with  Lamb  waves  which  have  traveled 
completely  around  the  shell  and  radiate  in  the  forward  direction.  A  direct  comparison  of  the 
amplitudes  of  the  Lamb  wave  and  edge  diffracted  bursts  in  Fig.  9  is  not  easily  made  with 
theory  because  of  complications  in  calculating  the  edge  diffracted  ray  amplitude  for  an 
observer  not  in  the  far  field.  Nevertheless,  from  the  arrival  time  and  appearance  of  the 
early  burst  in  Fig.  9,  it  can  be  concluded  that  forward  directed  Lamb  wave  contributions  to 
scattering  may  be  readily  observed.  It  was  also  observed^^-^^  that  the  amplitude  of  the 
early  Lamb  wave  burst  decreases  as  the  receiver  is  moved  off  the  axis  in  the  way 
characteristic  of  axially  focused  (or  glory)  scattering,2.5,35 

Consider  again  the  total  cross  section  and  form  functions  characteristic  of  steady 
state  scattering.  Inspection  of  Figs.  7  and  8  reveals  the  presence  of  a  broad  stmeture 
having  a  quasiperiod  Ax  ~  9.  This  structure  is  primarily  due  to  the  variation  in  phase  with 
ka  of  the //  contribution  with  /  =  flo  as  shown  by  the  analysis  which  follows.  From 
inspection  of  Eqs.  (3)  and  (14)  it  may  be  anticipated  that  a,  contains  a  suucture  periodic  in 
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Fig.  9  Time  record  of  on-axis  forward  scattering  of  a  653  kHz  4  cycle  tone  burst  from  a 
stainless  steel  shell.  The  contribution  labeled  FD  is  associated  with  ordinary  forward 
diffraction.  It  is  preceded  by  a  contribution  labeled  LW  attributed  to  Lamb  waves  because 
of  its  arrival  time.  The  observation  experimentally  confirms  the  existence  of  a  forward 
directed  Lamb  wave  contribution.  The  structure  following  the  FD  burst  is  at  least  partially 
attributable  to  repeated  circumnavigation  of  Lamb  waves,  though  there  may  be  some  weak 
scattering  from  a  three-pronged  wire  mount  which  supported  the  sphere. 
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ci.  anges  of  T)/  by  -2k.  (It  is  noteworthy  here  that  at  resonance3>9  expiilKxc/ct)  =  -1  and  the 
denominator  in  Eq.  (14)  has  a  specific  phase.)  From  Eq.  (13),  the  change  in  jc  associated 
with  All/  =  -271  may  be  approximated  as 

Ax^ «  27t[2  cos(8p  -  (c/cp(7c  -  20p]'^ ,  (18) 

provided  that  c/  depends  sufficiently  weakly  on  x  that  dispersion  may  be  neglected.  Of 
particular  interest  is  the  result  that  for  /  =  ao.  cilc  increases  only  from  1.95  to  2.04  as  x 
increases  from  50  to  100.  Taking  ka  ~  75  gives  ci/c  =  2.028  and  Eq.  (18)  predicts  Ax/  ~ 
9.0  in  agreement  with  the  quasiperiod  manifest  in  the  calculations.  The  apparent  relative 
importance  of  the  /  =  oq  contribution  to  the  underlying  structure  may  be  attributable  to  the 
relatively  large  value  of  the  associatexi  I  G/ 1  evident  in  Fig.  13  of  Ref.  5.  (The  I  =  sq 
contribution  may  also  produce  a  superposed  broad  structure  in  Of,  however,  it  is  more 
difficult  to  estimate  the  associated  quasiperiod  because  of  dispersion.)  Mote  that  the 
separation  Ax/  is  not  directly  caused  by  the  separation  of  specific  resonances  though  the 
finer  structure  in  C(  is  attributable  to  resonances.  When  the  surface  elastic  waves  are  only 
weakly  dispersive  the  spacing  between  resonances  may  be  approximated  as  ^sxf^^  «  Cgitc 
(see,  e.g.,  Ref.  7)  and  when  dispersion  is  negligible  this  becomes^  Ax/'"^^  =  c//c.  Hence 
for  /  =  tJQ  and  the  ka  region  considered  above  &xf^^  ~  2  which  is  much  less  than  the 
quasiperiod  of  the  broad  structure. 

Nusscnzveig  and  Wiscombe^^  note  that  the  structure  in  the  optical  O/  of  dielectric 
spheres  (which  is  quasiperiodic  in  ka)  arises  from  an  "interference  between  the  forward 
diffraction  peak  and  forward  glory  contributions  . . .".  Though  the  details  of  the  glory 
contributions  differ  here,  this  interference  condition  also  leads  directly  to  Eq.  (18)  for 
weakly  dispersive  Lamb  waves.  [This  may  have  been  anticipated  by  comparing  the  form 
of  Eq.  (18)  with  the  consequence  of  setting  in  Eq,  (17),  Ar/  to  an  integer  multiple  of  the 
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wave  period  Itc/co.]  Another  example  of  a  bioad  modulation  of  a  cross  section  is  that  of 
backscattering  from  a  tungsten  carbide  sphere  discussed  by  Williams  and  Marston.^  For 
certain  ka  regions,  resonance  reduces  the  total  amplitude,  but  at  other  regions,  it  causes  an 
increase.  In  that  case  the  relevant  interference  was  between  the  specular  contribution  and 
the  contribution  due  to  Rayleigh  waves  on  the  sphere. 

There  has  been  increased  interest  in  ray  treatments  of  wave  propagation  on 
cylindrical  sheUs.^>^>^0,4l  seems  appropriate  therefore  to  comment  on  the  conn,  ction 
between  ray  methods  for  leaky  surface  elastic  waves  on  cylinders  (Sec.  V  of  Ref.  9)  with 
the  present  results  for  the  forward  glory  of  spheres.  The  connection  is  as  outlined  in  Sec. 
VI  of  Ref.  9  where  it  is  shown  that  ray  results  for  cylinders  may  be  adapted  to  construct  the 
farfield  amplitude  due  to  surface  elastic  waves  on  spheres.  The  G/  used  here  in  Eqs.  (12) 
and  (14)  and  elsewhere  in  Refs.  1, 3, 5,  and  9  for  spheres  is  descriptive  of  farfield 
amplitudes.  The  coupling  coefficient  Gfy  for  cylinders  (Ref.  9,  Sec.  V)  is  at  leasi 
approximately  descriptive  of  near  field  as  well  as  far  field  amplitudes.  The  cylinder 
analysis  facilitates  the  approximation  of  the  local  amplitude  of  the  outgoing  toroidal 
wavefront  near  a  sphere  [Eqs.  (54)  and  (57)  of  Ref.  9]  and  the  approximation  of  the 
diffraction  integral  for  glory  scatlering^S  relates  those  local  amplitudes  to  the  far-field  form 
function//.  [See  Eq.  (55)  of  Ref.  9.]  It  is  tlie  intent  of  one  of  us  (P.L.M.)  to  give  a  more 
detailed  description  of  the  connection  of  results  for  spheres  and  cylinders  in  a  subsequent 

publication.'^^ 

It  is  noteworthy  that  since  the  optical  theorem  in  its  most  general  form^l>^2, 19-21 
relates  to  the  forward  scattering  amplitude,  it  is  applicable  to  objects  having  complicated 
shapes.  In  principle,  at  high  frequencies  ray  methods  could  be  developed  to  estimate 
forward  acoustical  amplitudes  and  give  approximations  for  Og  of  complicated  objects. 
Furthermore,  ray  methods  could  be  used  which  include  absorption  by  the  object. 
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Appendix  A.  Lamb  wave  damping  parameters,  phase  velocities,  and  the 
effect  of  the  breathing  mode 

'Fhe  Lamb  waves  considered  in  Ref.  5  are  the  lowest  antisymmetric  or  flexural 
wave,  (iQ,  and  the  lowest  symmetric  wave,  Sq.  The  material  and  geometrical  parameters  for 
the  particular  shell  studied  in  Ref.  5  are  given  here  in  Sec.  2.5  and  the  range  considered 
was  20  <  JL  <  75.  We  extended,  here,  the  range  of  x  =  lca  for  both  Oq  and  5o  to  7  <  < 

100.  Furthermore,  we  find  that  the  a  2  attd  ^2  waves  can  be  excited  on  the  elastic  shell  for  x 
greater  than  41  and  70,  respectively.  The  conoibutions  of  the  ai  and  ^2  waves  to  Gi  are 
discussed  in  Sec.  2.5. 

The  numerical  method,  implemented  in  the  determination  of  p/  and  c//c  is  given  in 
Appendix  A  of  Ref.  5.  'fhe  equuiions  to  be  solved  arc  based  on  applying  Watson 
transform  methodology  given  for  solid  elastic  spheresL^  to  the  case  of  a  spherical  shell.^*^ 
The  relevant  results  are  as  follows: 

Dy^(x)  =  0  ,  ^^2) 

V;  =  ct^  +  ipy ,  (A2) 

(c^  /c)  =  x/(a^  +  1/2) ,  (A3) 

Equation  (Al)  uses  the  determinant  in  the  denominator  of  the  PWS,  Eq.  (4),  for  the  elastic 
shell;  where  the  integer  index,  n,  has  been  replaced  everywhere  by  the  complex  index,  v/. 
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The  complex  index,  V;,  is  obtained  by  solving  Eq.  (Al)  at  fixed  values  of  x.  The  damping 
parameter,  p/  =  /m{v/),  is  immediately  obtained  (as  a  function  of  x)  for  the  Ith  type  Lamb 
wave.  Figure  Al  is  a  plot  of  v/  for  four  Lamb  waves  that  can  be  excited  on  the  shell 
considered.  It  should  be  noted  that  the  radiation  damping  is  not  necessarily  weak.^  In 
solving  Eq.  (Al),  Bessel  functions  of  the  first  and  second  kind  of  complex  order  and  real 
argument  y  must  be  calculated.  A  discussion  of  the  algorithm  employed  for  these  functions 
is  given  in  an  appendix  of  Ref.  43.  When  the  conditions  I  v/ 1  >  3  and  y/1  v/ 1  >  10'^  hold, 
then  the  calculated  values  are  asserted  to  have  errors  of  less  than  0.00001  +  i0.(X)O0L 
From  inspection  of  Fig.  Al,  it  is  found  that  for  the  sq  and  ai  Lamb  waves  the  condition  I  v/ 

I  >  3  is  violated  at  small  ka.  However,  the  computations  presented  in  Sec.  2.5  suggest  that 
the  accumulated  error  from,  these  two  contributions  are  not  significant.  We  did  not  search 
for  roots  corresponding  to  a  generalization^  of  rigid-sphere  Franz  waves  to  the  shell  since 
those  would  be  expected  to  have  large  p/  and  c/  near  c. 

The  normalized  phase  velocity  is  determined  by  using  Eq.  (A3).  The  dispersion 
curves  for  the  various  Lamb  waves  are  displayed  in  Fig.  A2.  There  are  three  features  that 
are  important  to  note.  First,  the  normalized  velocities  for  the  Uj  and  a’j  Lamb  waves  appear 
to  diverge  as  they  approach  their  cut-off  frequencies.  (These  cut-off  frequencies  are 
analogous  to  the  cut-off  frequencies  for  the  propagation  of  particular  modes  in  a 
waveguide.  There  is,  however,  additional  structure  near  the  cut-off  of  the  si  mode,  the 
detail  of  which  is  not  shown  in  Figs.  Al  and  A2.)  Second,  the  c/  for  I  ~  jq  appears  to 
diverge  as  (o  approaches  the  natural  frequency  cub  of  the  lowest  purely  radial  or  "breathing" 
mode  of  the  shell.  'Phis  divergence  is  plausible  since  radial  motion  of  a  shell  can  have  the 
same  phase  at  all  surface  points  when  co  =  ®b-  For  a  thin  spherical  shell  in  a  vacuum,  (Ob 
is  related  to  the  lowest  purely  radial  or  "ring"  frequency  cur  of  an  infinite  cylindrical  shell 
of  the  same  material  via  the  approximation 
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for  the  oq,  ^i.  Lamb  waves  on  a 
ci  are  respectively  the  solid,  short-dashed,  long- 
le  arc  length  along  each  curve  is  monotonic  in  ka 
ndicated. 
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g.  A.  Normalized  phase  velocMes  aloag  ,he  outer  surface  of  the  shell,  for  ao.  so, 
at.  and  s,  Lamb  waves.  The  legend  is  as  in  Fig.  A1  except  that  adjacent  to  each  cttrve  is 
one  with  shoa  dashes  which  was  calculated  front  plate  theoo-  coneoted  approximately  for 
curvature  as  described  in  Ref.  «  and  Fq.  tA8),  Note  that  the  hue  c,  /  =  so,  appears  to 
diverge  as  the  ka  of  the  breathing  mode  is  approached. 
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co3  =  [2(1+V)]^^\, 


(A4) 


where  V  denotes  Poisson's  ratio  which  for  440c  stainless  steel  is  0.296.  Equation  (A4) 
follows  from  Eq.  (7, 1 14)  of  Ref.  32.  For  a  thin  shell^^ 

>  ""pL  =  =  5.31  mnVHs,  (A5) 

where  E  denotes  Young's  modulus,  cpi  is  the  low-frequency  phase  velocity  of 
compressional  waves  in  an  elastic  plate  and  a  =  (a  +  b)/2  is  the  middle  radius  of  the  shell. 
(The  combined  result  for  cog  also  follows  from  Love."^)  If  the  effects  of  fluid  loading  on 
resonance  frequency  are  completely  neglected,  the  estimated  ka  of  the  breathing  resonance 
is 


Xg  [2(1  +  V)]^^^[l  -  (h/2a)]’\cpj^/c)  =  6.3  , 


(A6) 


where  /t  =  a  -  is  the  shell  thickness.  For  a  shell  surrounded  by  water,  fluid  loading 
reduces  jcb  by  a  fraction  of  order 

mg/Psh  =  (p/pg)  (aAi)  x^, 

where  mB  =  pa  is  the  relevant  "model  accession  to  inertia"  ratio. ^2  since  m^p$h~ 

0.018  «  1  for  the  sphere  .considered,  fluid  loading  should  cause  xb  to  be  reduced  only 
slightly  from  the  estimate  based  on  (A4).  The  third  feature  is  a  cut  off  frequency  for  the  si 
wave.  It  is  somewhat  analogous  to  the  cutoff  near  the  breathing  mode  since  it  lies  close  to 
the  condition  ka  «  Kci/c{l  -  b/a)  ~  77  of  the  lowest  thickness  resonance  discussed  in  Sec. 
Ill  of  Ref.  5. 

The  shortest  dashed  curves  in  Fig.  A2  are  the  results  of  introducing  a  curvature 
correction'^3  j^to  Lamb's  equations  (given  in  Ref.  5)  for  a  flat  plate  in  a  vacuum.  The 
thickness  of  the  plate  is  the  thickness  h  =  a(l  -  h/a)  3.1  mm  of  the  shell.  The 
approximate  curvature  correction  gives  t'»:e  phase  velocity  ratio  as'^^ 


52 


(c//c)  =  [(cpp/c][l-(h/2a)]  . 

where  (ci)p  is  the  phase  velocity  of  the  Ith  Lamb  wave  on  the  plate.  The  comparison  in  Fig. 
A2  shows  that  a  useful  initial  estimate  of  the  phase  velocity  ratio  can  still  be  obtained 
through  Eq.  (A8)  and  the  plate  equations.  The  deviation  of  the  5o  result  from  (A8)  with  the 
proper  value  from  (A3)  is  a  manifestation  of  the  hoop  stresses  associated  with  the  breathing 
mode  discussed  above.  The  difference  between  the  curves  for  both  the  ai  and  si  Lamb 
waves  may  be  due  to  the  radiation  loading.  The  curvature  of  the  plot  of  c//c  with  I  =  i’l  at 
the  mode  onset  may  be  associated  with  the  relatively  large  value  of  P/  evident  in  Fig.  Al. 

Root  computations  not  shown  in  Figs.  Al  and  A2  indicate  that  there  is  a  mode  for 
which  the  low-frequency  cut  off  is  at  «  83  for  which  c//c  always  exceeds  8  in  the  region 
between  this  cut  off  andjc  =  100.  We  have  tentatively  identified  this  mode  as  the  S2  mode. 
Computations  indicate  that  it  does  not  significantly  affect  Ot  in  the  region  considered, 
evidently  due  to  the  smallness  of  the  c/q  factor  in  Eq.  (15)  and  the  magnitude  of  p/. 
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Chapter  3 

Longitudinal  resonances  in  the  form  function  for  backscattering  from  a 
spherical  shell:  Fluid  shell  case 


Abstract 

Reverberations  of  longitudinal  waves  in  a  hollow  shell  can  strongly  affect  the 
backscattering  amplitudes  at  frequencies  associated  with  a  thickness  resonance.  The 
phenomena  is  studied  for  the  idealized  case  of  vanishing  shear  stresses  in  the  shell  material 
by  taking  that  material  to  be  an  inviscid  fluid.  The  sound  speed  ci  for  the  fluid  is  taken  to 
be  that  of  longitudinal  waves  in  elastic  materials  of  interest;  the  suirounding  fluid  being 
water.  An  exact  partial-wave  series  gives  the  form  function/for  backscattering  and  plots 
of  1/ 1  as  a  function  of  ka  display  resonance  features  where  a  denotes  the  outer  radius  of  the 
shell.  Tliese  features  are  also  recovered  in  a  direct  geometrical  calculation  of/ which  sums 
the  amplitudes  associated  with  rays  multiply  reflected  within  the  curved  shell.  This 
geometric  synthesis  shows  that  the  effects  of  curvature  are  essential  to  modeling/.  In 
addition  to  numerical  comparisons  with  the  partial  wave  series,  the  geometrical  calculation 
is  tested  by  considering  several  limiting  cases  and  results  anticipated  from  elementary 
consideration  are  recovered. 
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3.1  Introduction 

Various  authors  have  considered  how  resonances  affect  tlic  reflection  of  sound 
from  flat  elastic  plates^  and  fluid  layers.^  Others  have  considered  the  significance  of 
surface  guided  wave  contributions  to  the  backscattering  of  an  incident  plane  wave  from 
elastic  spherical^"^  and  circular  cylindrical  shells.^'*  The  scattering  amplitude  can  be 
represented  in  terms  of  contributions  from  surface  guided  waves  which  circumnavigate  the 
shell,  continually  radiating  energy  into  the  surrounding  water.  This  surface  guided  wave 
representation  can  give  a  simple  quantitative  description  of  resonances  which  has  been 
studied  in  detail  for  solid  spheres.^’l^  When  exploring  such  a  representation  for  elastic 
shells,^!  we  have  seen  a  justification  for  investigating  the  reverberation  of  bulk  waves 
transmitted  through  curved  shells  which  reflect  from  the  inner  surface.  Associated  with 
transmitted  (bulk)  longitudinal  waves  are  thickness  resonances.^ To  isolate  the 
contribution  of  reverberation  of  transmitte<i  longitudinal  waves  from  other  contributions, 
we  restrict  our  attention  to  backscattering  from  an  evacuated,  idealized  spherical  fluid  shell 
in  water.  This  simplification  allows  a  quantitative  ray  representation  of  tlie  scattering 
amplitude  to  be  developed  and  tested  without  considering  contributions  associated  with 
transverse  waves.  The  testing  of  ray  models  for  fluid  objects  can  be  justified  because  of 
insight  gained  into  the  behavior  of  more  complicated  situations. 

The  usual  condition  for  the  existence  of  a  longitudinal  resonance  is  that  the 
thickness  of  an  empty  shell  should  be  an  integral  number  of  kJI  where  A./,  is  the 
longitudinal  wavelength  within  the  fluid.  This  condition  is  applicable  when  the  impedance 
pgCjL  of  the  shell  material  exceeds  pc  of  the  surrounding  fluid.  For  an  incident  acousuo 
plane  wave  with  angular  frequency  the  resonance  condition  can  be  c-'iiiten  as 


k^h  =  nTt. 


(n  ==  1,2,  ...), 
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where  ki  =  (H/ci  is  the  wavenumber  and  ci  is  the  speed  of  sound  in  the  fluid.  3  As  will  be 
shown  below,  the  manifestation  of  thickness  resonances  is  contained  witliin  the  specular 
reflection  form  function  f^p  and  is  intimately  related  to  the  curvatures  of  the  inner  and  outer 
surfaces  of  the  shell.  Examination  of  I  fsp  I  shows  that  the  effects  of  curvature  are  largest 
for  those  k^h  satisfying  Eq.  (1). 

The  organization  of  this  paper  is  as  follows.  A  derivation  of  a  ray  acoustical 
representation  of  the  form  function  is  summarized  in  Sec.  3.2.  This  geometrical  synthesis 
of  fsp  decomposes  into  a  contribution  for  the  reflection  from  a  vacuum-backed  flat  plate  and 
a  curvature  dependent  component, /^c.  In  Sec.  3.2,  the  synthesis  is  compared  with  exact 
computations  employing  the  partial-wave  series  (PWS)  representation  of  the  form  function 
for  backscatter.  Section  3.3  examine  the  significance  of  the  curvature  correction 
With.n  Sec.  3.3,  the  most  significant  contributions  to/cc  obser/ed  to  occur  at  the 
resonance  condition  Eq.  (1).  To  further  support  the  claim  that  the  form  obtained  for  fee  is 
correct,  Sec.  3.4  considers  several  limiting  cases  of  the  fluid  shell  parameters.  Each 
limiting  case  produces  the  expected  result  and  thtJS  supports  the  geometrically  derived  foim 
of  the  curvature  correction.  Section  3.5  contains  concluding  remarks  and  notes  tlie 
importance  of  fee  to  die  scattering  of  an  acoustic  plane  w;  ve  from  an  elastic  spherical  shell. 

The  article  contains  t'nree  apfjendices.  Appendix  A  contains  the  essential  details  of 
the  matrix  method  employed  in  detemiining  the  geometrical  spreading  factor  and  amplitude 
coefficient  introduced  in  Sec.  3.2.  As  noted  in  Appendix  A,  the  matrix  method  is  similar  to 
the  paraxial  matiix  methods  of  geometric  optics.  Appendix  B  i:;  included  for  completeness 
and  contains  information  pertinent  to  the  PWS  representation  of/  for  backscattering  from  a 
fluid  shell.  Finally,  Appendix  C  describes  a  technique  which  rapidly  sums  the  infinite 
series  in/c,;.I‘^  An  extunple  of  the  rapid  summation  technique  is  compared  to  the  exact 
summation  of  the  series  where  the  convergence  of  the  series  is  explicitly  tested. 


3.2  The  specular  reflection  for  a  fluid  shell  from  geometric  ray  methods. 

Consider  an  empty  spherical  shell  of  inviscid  fluid  embedded  in  a  fluid  (water). 

The  outer  radius  of  the  shell  is  denoted  by  a  and  the  inner  radius  is  given  by  b.  ITie  interior 
of  the  fluid  shell  is  taken  to  be  a  vacuum.  That  is,  the  inner  surface  at  6  is  a  pressure 
release  surface.  The  relevant  parameters  for  the  water  are  the  sound  speed,  c,  and  the 
density,  p.  The  sound  speed  and  density  of  the  fluid  are  and  pe,  respectively.  The 
scattered  pressure  in  the  farfield  is 


(2) 


where  the  harmonic  time  dependence  exp(-iox)  has  been  suppressed The  wavenumber 
of  the  incident  plane  wave  in  water  isk-  (o/c  =  27tA  where  0)  is  the  angular  frequency  and 
X,  is  the  wavelength.  The  pressure  amplitude  of  the  incident  plane  wave  is  po  and  the 
distance  from  the  center  of  the  sphere  to  an  observation  point  is  /.  The  scattering  amplitude 
or  form  function,/,  can  be  approximately  partitioned  asf  ~fsp  ^  fsw  where  fsp  is  a 
contribution  from  a  specular  reflection  andfi-^  is  a  term  associated  with  the  possible 
excitation  of  surface  guided  waves.  As  stated  above,  Sec.  3.2  is  concerned  with  a 
geometric  synthesis  of fsp. 

Figure  1  is  a  simplified  ray  diagram  for  the  backscattering  of  an  acoustic  plane  wave; 
from  an  ideal  fluid  spherical  shell.  The  plane  wave  is  assumed  to  be  traveling  in  the 
positive  z-direction  and  the  origin  of  the  coordinate  system  is  located  at  the  center  of  the 
shell  at  point  O.  Tie  z-axis  coincides  with  the  line  SO.  Tie  specular  reflection  appears  to 
come  from  the  point  S.  This  ixiint  may  be  referred  to  as  the  specular  point.15  ^  ray, 
perpendicular  to  the  wavefront,  with  impact  parameter  s  piopagatcs  from  point  A  to  point  B 
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Fig.  1  The  ray  diagram  for  the  geometric  synthesis  of  tite  specular  reflection  form 
function  from  a  vacuum-filled  fluid  shell  The  point  S  is  the  vertex  of  the  refracting  surface 
and  O  is  the  origin  of  a  coordinate  axis  located  at  the  center  of  the  shell.  The  inner  and 
outer  radii  are  denoted  by  b  and  a,  respectively.  A  ray,  infinitesimally  close  to  the  z-axis, 
SO,  is  incident  on  the  fluid  shell  with  impact  parameter  s.  The  ray  ABC  is  the  normal 
specular  ray;  while  rays  ABDEF  and  ABDEGHl  are  the  first  two  internal  reflection 
contributions.  At  the  points  B,  E  and  H  partial  reflection  and/or  transmission  of  a  given  ray 
occurs.  The  internal  rays  are  totally  reflected  at  all  points  on  the  inner  suiface  r  =  b. 
Intersections  of  the  projection  (dashed  lines)  of  the  outgoing  rays  (flC,  FE  and  HI)  and  the 
2-axis  define  locations  of  virtual  point  sources  which  detemiine  local  curvature  of  the 
outgoing  wavefront  associated  v.'ith  a  given  ray.  The  virtual  point  sources  are  denoted 
(n  =  0,  1, ...)  and  the  local  backscattering  angles  are  Y„.  Superposition  of  the  acoustic 
wavefields  from  the  virtual  sources  in  the  limit  0  gives  an  expression  for  the  specular 
reflection  contribution  to  the  form  function  for  backscatter  as  described  in  Sec.  3.2  and 
Appendix  A. 
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on  the  surface  of  the  shell.  Rotation  of  the  ray  diagram  about  the  line  SO  generates  all  rays 
having  an  infinitesimal  impzct  parameter  s.  At  B  the  ray  is  partially  reflected  back  into  the 
water  and  partially  transmirted  into  the  fluid  shell.  ITie  intersection  of  the  reflected  ray  BC 
(dashed  line)  and  the  z-axis  at  Vq  defines  the  location  of  a  virtual  source.  For  backwards 
and  near  backwards  directions,  spherical  aberration  may  be  neglected  and  this  virtual 
source  may  be  modeled  as  a  point  source.  With  the  appropriate  amplitude  and  phase  shift, 
this  virtual  point  source  desciibes  the  propagation  of  the  local  outgoing  wavefront 
associated  with  ray  BC.  The  transmitted  ray  at  B  is  refracted  and  propagates  along  the  path 
BD.  After  total  internal  reflection  at  D,  the  ray  propagates  to  E  and  is,  again,  partially 
reflected  and  tritnsmitted.  The  transmitted  ray  EF,  after  refraction,  can  be  projected  back  to 
a  second  virtual  source  located  at  V/.  'Fhe  partial  reflection  and/or  transmission  and 
refraction  of  subsequent  rays  at  the  outer  surface  generates  an  infinite  set  of  virtual  point 
sources  located  at  points  V„.  Hence,  each  virtual  point  source  describes  the  local  outgoing 
wavefront  associated  with  a  specific  ray. 

For  each  point  define  the  local  backscattering  angle,  y^,,  as  the  angle  between 
the  outgoing  ray  and  the  z-axis  (sec  Fjg.  l).  As  the  impact  parameter  approaches  zero,  the 

local  harlfSrattprincr  analp  V_  a.«r»  oripc  in  wrn  I-Ipnr-p  ac  c  O  ihp  nnntrihiitinn  r\f  thp 

- - ^  I  *'W  wxrs-w*  .a.  vb  •  v  f  Vj  V*  Vra*  V*  bAaW 

specular  reflection  to  the  backscattered  farfield  pressure  can  be  expressed  as  a  superposition 
of  the  acoustic  fields  from  the  vinual  point  sources.  In  particular,  the  farfield  specular 
pressure  is 


Psp 


/  oo 

n  =  1 


HnCnC 


(3) 


where  h~  {a-b)\s  the  thickness  of  the  shed  and  x-kai'i'd.  dimensionless  size  parameter. 
The  2x  phase  shift  in  Fq.  (3)  accounts  for  the  path  length  difference  between  a  ray 
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propagating  in  water  from  point  S  in  Fig.  1  to  O  and  back  lo  S  and  a  ray  which  is  reflected 
in  the  backward  direction  at  S.  The  quantity  74  is  a  geometric  spreading  factor  associated 
with  the  nth  virtual  source.  An  expression  for  is  derived  in  Appendix  A  using  a 
geometric  method  analogous  to  the  matrix  methods  of  geometric  optics.  'Fhe  factor  Cn 
accounts  for  the  partial  reflection  and/or  transmission  of  the  ray  at  each  surface  (see 
Appendix  A).  Tlie  exponential  factor  within  the  summation  is  the  additional  phase  delay  of 
a  wave  which  revert*erates  n  times  within  the  shell  before  being  transmitted  in  the 
backscattered  direction.  Th  longitudinal  wavenumber  is  =  co/c/^  =2kIXl  where  is 
the  wavelength  of  the  acoustic  field  in  the  fluid.  Finally,  ris  the  reflection  coefficient  for  an 
acoustic  plane  wave  at  nomial  incidence  on  a  flat  liquid-liquid  interface.  The  reflection 
coefficient  is 


Pe^L  ~  pc 
Pc^L  +  pc  ’ 


(4) 


for  the  above  material  parameters. 

Using  Eqs.  (A7)  and  (A8)  for  7/„  and  Cn,  the  contribution  of  the  specular  reflection 
to  the  backscattered  farfield  pressure  is 
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where  a  =  2kih  =  2x{clci){i  -  bid).  Equation  (A6)  defines  li  which  depends  only  on  the 


relative  index  of  refraction  Mi  =  dci  and  the  ratio  of  radii  bja.  Comparison  of  Eq.  (5)  with 
Eq.  (2)  reveals  that  the  specular  reflection  contribution  to  the  backscattered  fomi  function  is 
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Although  Eq.  (6)  is  a  geometrical  synthesis  of  the  specular  reflection  fomi  function,  it  can 
be  expressed  in  a  manner  more  easily  understood.  Simple  algebraic  manipulations  allow 
one  to  rewrite  Eq.  (6).  First,  the  factor  [1+  nB]'^  is  expressed  as  1  -  nB[\+  nB]'^.  With 
this  substitution,  the  infinite  series  in  Eq.  (6)  splits  into  two  series.  The  summand  of  the 
first  of  the  resulting  series  is  (re'o^)"  where  I  re‘“  1  <  1.  Comparing  this  series  with  the 
geometric  series  demonstrates  that  the  first  series  is  summable  and  is  equivalent  to  re‘“[l-  r 
e'^]'T  Substitution  into  Eq.  (6)  leads  to 

f,p=«5'‘^’‘+fcc  m 


where 
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Equation  (7)  is  the  desired  expression  for  the  contribution  of  the  specular  refletnion  to  the 
backscattered  form  function.  Define.  ^  =  ^Rsxp(,-i2x);  the  coefficient  i^is  the  complex 
reflection  coefficient  for  a  vacuum-backed  flat  plate  of  thickness  h  at  normal  incidence.^ 
Tlie  phase  shift  exp{-i2x)  is  due  to  the  phase  reference  at  the  shell's  center  used  in  Eq.  (2). 
Hence  fp  describes  the  reflection  from  a  spherical  shell  where  the  complex  reflectivity  is 
modeled  as  that  for  a  flat  plate.  Consequently,  Eq.  (8b)  accounts  for  the  correction  to  f^p 
due  to  the  curvature  of  the  shell  within  the  geometric  approximation.  At  present,  an 
analytical  expression  for  fee  is  not  available,  however,  a  rapid  summation  teclmique  can  be 
employed  to  perform  the  indicated  summation  (see  Appendix  C). 

Figures  2  and  3  are  comparisons  of  the  exact  PWS  calculation  of  I/I  with  the 
geometrical  synthesis  as  given  by  Eqs.  (7)  and  (8).  Appendbt  B  contains  the  PWS 
representation  for  /and  some  relevant  remarks  concerning  its  computation.  Tlicse 


Fig.  2  The  solid  curve  is  the  exact  PWS  compulation  of  the  form  function  for 
backscattering  for  a  4%  thick  aluminum  fluid  shell  with  the  parameters  listed  in  Table  I. 
The  geometric  synthesis  of  the  specular  reflection  contribution  is  the  dashed  line  and  is 
calculated  via  Eqs.  (7)  and  (8).  The  broad  minimum  at.x  =  340  is  the  lowest  longitudinal 
resonance.  Without  the  curvature  correction  (as  discussed  in  Sec.  3.3),  the  synthesis 
would  not  model  the  longitudinal  resonance  since  the  synthesis  simplifies  to  If^p  1=1- 
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Fig.  3  The  exact  PWS  and  geometric  synthesis  are  the  solid  and  dashed  lines, 
re",  lectively.  These  computations  are  for  a  16.2%  thick  440c  stainless  steel  fluid  shell  (see 
Table  I).  The  broad  dip  at  x  =  76.8  corresponds  to  the  %st  longitudinal  resonance.  The 
oscillatory  behavior  in  the  ranges  0  <  jt  <  40  and  11  <x<  100  is  presumably  attributable  to 
waves  in  the  fluid  shell  shedding  energy  into  the  water  while  circumnavigating  the  shell,  as 
such  contribu^’ons  were  not  included  in  the  synthesis. 
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computations  are  intended  to  give  insight  into  the  more  complicated  problem  of  resonances 
in  backscattering  from  elastic  shells,  so  the  choice  of  the  fluid  parameters  corresponds  to 
elastic  material  parameters  instead  of  parameters  for  actual  liquids.  Tlie  parameters  for  the 
fluid  shells  and  surrounding  water,  used  in  the  various  calcinations,  are  listed  in  Table  I. 
Figure  2  corresponds  to  an  aluminum  fluid  shell  with  bla  =  0.96  where  the  solid  line  is  the 
PWS  result  and  the  dashed  line  is  the  geometrical  synthesis.  The  resonance  condition,  Eq. 
(1),  predicts  the  first  thickness  resonance  occurs  at  Jt  =  340.  Clearly,  the  broad  minimum 
in  Fig.  2  corresponds  to  this  longitudinal  resonance.  The  synthesis  of  \fsp  I  and  1/ 1  via  the 
PWS  are  in  e:  cellent  agreement  for  »  0. 

A  440c  stainless  steel  fluid  shell  with  bla  -  0.838  has  its  first  longitudinal 
resonance  at.t «  76.8  (see  Fig.  3).  'Fhe  material  pai'ameters  of  440c  stainless  steel  and  this 
particular  value  of  the  radii  ratio  correspond  to  the  properties  of  a  spherical  shell  studied 
extensively  in  experiments  with  tone  bursts.^  As  in  Fig.  2,  the  solid  and  dashed  lines  are 
for  the  PWS  and  geometrical  synthesis,  respectively.  Again,  the  broad  dip  at  a:  «  76.8  in 
Fig.  3  is  a  manifestation  of  the  longitudinal  resonance  and  the  synthesis  is  in  good 
agreement  with  the  PWS  calculation.  The  oscillatory  behavior  observed  in  Fig.  3  is 
attributable  to  the  symmetric  (or  vi)  siiiface  guided  wave  which  contributes  significantly  to 
the  backscattering  for  x  >  76.8.  That  is,  the  contribution  of fsw  to  I/I  is  no  longer 
negligible  for  the  440c  stainless  steel  fluid  shell.  The  shell  acts  like  a  curved  fluid 
waveguide  which  leaks  energy  back  into  the  surrounding  water. 

For  fluid  shells  where  <  pc,  the  resonance  condition  in  Eq.  (1)  does  not 
apply.  The  proper  expression  for  predicting  a  thickness  resonance  is  kih  ~  (n  +1/2)k,  (n  = 
0,  1, ...).  Since  internal  rays  reflected  at  r  =  a  do  not  sustain  a  K  phase  shift,  then  the 
resonance  requirement  of  the  internal  rays  adding  in  phase  leads  to  the  additional  7i/2. 
Although  the  resonance  condition  is  different  for  p,;C/^  <  pc,  Eqs.  (7)  and  (8)  still  describe 


TABLE  I.  Material  parameters  used  in  the  calculations. 


Pe 

CL 

c 

P 

Fluid/Water 

(g/cm^) 

(km/s) 

(km/s) 

(g/cm3) 

aluminum/water 

2.70 

6.42 

1.4825 

1.00 

440c  stainless  steel/water 

7.84 

5.854 

1.479 

1.00 
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the  specular  reflection  contribution  to  the  form  function  for  backscattering.  Figure  4  is  a 
comparison  of  the  PWS  result  and  the  synthesis  using  Eqs.  (7)  and  (8).  The  fluid  shell  and 
water  parameters  are  the  same  as  the  440c  stainless  steel  parameters  with  b/a  ~  0.838 
except  Pe  =  0. 1  g/cni^.  The  modified  condition  predicts  the  first  thickness  resonance  axx~ 
38.4,  i.e.  half  the  resonance  value  for  the  high  impedance  case  of  76.8.  Inspection  of  Fig. 
4  demonstrates  tliat  the  geometrical  model  agrees  with  the  exact  PWS  result. 

3.3  The  importance  of  curvature  to  |  fsp  | 

An  examination  of  the  curvature  correction/cc  is  merited.  Tlie  significance  of  the 
shell's  curvature  to  the  manifestation  of  thickness  resonances  can  be  elucidated  by 
examining  its  dependence  on  kih.  First,  it  is  noted  that  the  flat  plate  contribution.^, 
alone,  does  not  reproduce  the  longitudinal  resonance  structure  observed  in  \f\  since 
1.^  I  =  1 .  Tlie  approximation  I/I  » \fp  I  does  not  exhibit  stnicture  at  values  of  x  which 
satisfy  Eq.  (1).  To  obtain  the  minima,  the  curvature  correction /tt  needs  to  be  included 
such  that  1/1*1,^  +fcc  I-  "rhe  resonance  condition  indicates  that  I  fee  I  should  have  its  first 
maximum  at  kih  =  k  when  peCL  >  pc.  For  a  specific  fluid,  the  curvature  correction 
depends  on  the  two  parameters  x  and  b/a.  For  the  discussion  below,  the  material 
parameters  are  those  of  440c  stainless  steel. 

In  Fig.  5, 1  fee  I  is  plotted  as  a  function  of  k/^h  =  x/ifi^(l  -  b/a)-,  where  x  =  20  or  60  is 
held  fixed  while  0  <  b/a  <  1.  For  fixed  x,  o/a  =  1  corresponds  to  kih  ~  0  and  as  b/a  — >  0, 
kih  xMi.  The  peaks  in  Fig.  5  are  at  the  expected  value  kih  =  ic  and  tliese  peaks  tend  to 
be  narrow  resonances  in  the  variable  ki.h-  The  interpretation  is,  tliat  for  fixed  x,  shells  v/ith 
h  in  the  vicinity  of  Xi/2  can  support  the  first  thickness  resonance.  Computations  for  larger 
values  of  kih  exhibit  similar  sharp  peaks  at  the  higher  order  longitudinal  resonances. 
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Fig,4  The  parameters  here  are  like  Fig.  3  except  that  pe  =  0.1  gtcm^.  The  condition  in 
Eq.  (1;  no  longer  applies  for  this  case  since  pgC/,  <  pc.  The  correct  resonance  condition 
becomes  k^h  ~  {n  +  l/2)r.,  (n  =  0, 1 ,...)  and  predicts  a  longitudinal  resonance  at  x  ^  38,4. 
Comparison  of  the  PWS  (dashed)  and  the  synthesis  (solid)  shows  a  minimum  atx  -  38.4. 
Tlie  material  ptirarneters  for  these  computations  are  given  in  Sec,  3.2. 


Fig.  5  The  magnitude  of  the  curvature  correction  fcci^iA)  for  fixed  values  of  x  and 
variable  b/a.  The  solid  curve  is  for  x-~20  and  the  dashed  line  is  for  x  =  60.  The  shatp 
peaks  occur  at  the  expected  value  kih  =  7t  where  the  first  longitudinal  resonance  is 
supported  within  the  fluid  shell.  The  material  properties  of  the  shell  coirespond  to  440ct 


stainless  steel. 
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Thus,  a  thickness  resonance  can  be  excited  within  a  fluid  shell  at  fixed  x  pro  vided  /z  is  in 
the  vicinity  of  an  integral  number  of  X//2  . 

Tlie  resonance  condition  still  holds  if  Ufa  is  now  held  constant  and  tlic  frequency  of 
the  incident  plane  wave  is  vailed.  Figure  6  displays  I  fee  I  where  bla  =  0.838  (solid  line)  or 
0,920  (dashed  line)  and  x  ~  ka  --  coa/c  is  varied.  The  linear  relation.ship  between  z  and  kih 
implies  that  as  x  increases  kik  increases.  Although,  the  sharp  peaks  in  Fig.  <5  occur  at  kj^h 
-■  n,  Ui:e  interpretation  of  the  physical  nature  of  the  peaks  is  more  subtle.  Inspection  of  Fig. 
6  reveals  that  the  width  of  the  peaks  are  comparable;  however,  the  x  values  at  which  the 
two  peaks  occur  are  different.  These  x  are  appioximatcly  76.8  and  155  for  bla  ~  0.838  and 
0.920,  respectively.  From  the  resonance  condition,  it  is  inferred  that  the  ;  ange  of  x,  for 
which  a  longitudinal  resonance  is  excited,  is  wider  for  a  thin  shed  than  a  thick  shell  of  the 
same  fluid.  Furihennore,  comparison  of  Figs.  2  and  3  illitstrates  tlic  width  in  x  is  strongly 
affected  by  the  parameter  ^^(1  -  bla). 

3.4  Limititiig  cases  for  scattering  of  an  acoustic  plane  wave  from  an  ideal, 
fluid  shell 

A,  Scattering  from  a  perfectly  soft  bubble 

For  the  limiting  case  b  a,  the  fluid  shell  becomes  a  perfeciiy  soft  bubble  and  Eq, 
(7)  should  describe  the  linear  backscattcring  of  an  acoustic  plane  wave  from  a  bubble  of 
radius  a  within  the  geometric  limit.  For  b  a,  inspection  of  Eq.  (A6)  reveals  the  constant 
B  —^0  which  leads  to  the  result  that  the  curvature  concction  component  in  Eq.  (7)  also 
vanishes.  In  this  limit,  it  is  stiaightlbrward  lo  show/y^  =  fp  -  -exp(42x).  In  Sec.  3.2,  ifc 
was  found  to  be  the  complex  reflection  coefficient  for  a  vaeiniin- backed  flat  plate.  Since 
the  thickness  of  the  plate  is  approaching  zero,  then  the  plate  becomes  a  water-vacuum 
interface.  The  reflection  coefficieni  for  a  nonnal  incidence  plane  wave  in  water  on  a 
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pressure  release  surface  is  -1 .  Hence,  Eq.  (8a)  reduces  to  the  expected  value  and  the 
backscattered  farfield  pressure  due  to  the  specular  reflection  becomes 


D  =-n 

Psp  Po  2r 


(9) 


The  upper  two  curves  in  Fig.  7  are  1/ 1  and  \fsp  I  where  the  dashed  line  is  the  PWS  result 
within  the  limit  b-^a  and  the  horizontal  line  is  the  magnitude  of  the  form  function 
obtained  from  the  acoustic  ray  method.'®  Both  I/I  and  I  fsp  I  are  scaled  by  bla  =  0.838  in 
Fig  7.  Tlie  reason  for  this  scaling  will  be  discussed  in  Sec.  3.4B.  For  x  >  20,  these 
curves  become  indistinguishable;  while  for  values  of  jf  <  20,  the  PWS  deviates 
substantially  from  the  synthesis.  This  result  may  be  anticipated  since  geometrical  acoustics 
is  implicitly  a  high-frequency  approximation. 

n.  Scafrcriiig  front  ati  iinpcilaiice-inalched  fluid  spherical  shell 

For  the  scattering  tforn  an  iiiipedance-matcli'xl  Ihiid,  llit-  rcrieciion  coefficient 
defined  in  Eq.  (4)  is  identically  zero.  That  is,  =  pc  for  an  impedance-matched  fluid 
and  r~  0.  Equation  (7)  reduces  to 
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where  B  is  defined  by  Eq.  (A6).  Figure  7  contains  computations  of  I /I  and  l/jp  I  for  a 
fluid  with  the  sound  speed  of  440c  stainless  steel  and  pe  was  determined  from  Pe  =  pM^. 
Tlie  lower  two  curves  in  Fig.  7  are  for  this  limiting  case.  The  dashed  line  is  the  exact  PWS 
result  and  the  horizontal  line  corresponds  to  Eq.  (!9).  Although,  the  two  ctuves  ai’e  vastly 
different  forx  <  40,  for  x  ~  50  the  difference  between  the  ray  tracing  s>’nthesis  and  the 
PWS  is  less  than  7.2%.  It  is  anticipated  that  the  agreement  wall  improve  as.x  becomes 
large. 
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Fig.  7  The  upper  two  curves  are  for  the  limiting  case  of  tlie  scattering  from  a  perfectly 
soft  'nubble.  'Oie  horizontal  line  is  \fsp  I  and  the  dashed  line  is  1/1  calculated  via  the  PWS 
result.  Both  curves  have  been  scaled  by  0.838  (see  Sec.  3.4B).  The  lower  set  of  curves  is 
the  limiting  case  of  an  irnpedance-rnatched  fluid  (i.e.,  pgCf,  -  pc).  Agai.n,  the  horizontal 
line  is  I  Jsp  1  and  the  dashed  line  is  the  PWS  for  I  ^  I.  The  two  sets  of  curves  demonstrate 
that  the  synthesis  of fsp  correctly  appioximates  the  high-frequency  scattering  for  these 
limiting  cases. 
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Within  the  impedance-matched  fluid  case,  there  is  the  special  case  Pe  ~  p  and  ci  = 
c.  That  is,  the  fluid  shell  is  essentially  a  perfectly  soft  bubble  of  radius  b.  The 
denominator  of  Eq.  (10)  becomes  (1  +  B)  =  alb  and  the  phase  shift  a  =  2kh.  Equation  (10) 
in  this  spex;ial  case  reduces  to  fsp  =  -(Jbla)e'^^  which  gives 


Psp- 


n  b  --i?.kb  ikj 

■P0  2r®  ®  > 


(11) 


for  the  contribution  of  the  specular  reflection  to  the  backscattered,  farfield  pressure. 
Comparison  of  Eqs.  (9)  and  (11)  gives  the  motivation  of  the  scaling  of  the  upper  two 
ciii-ves  in  Fig.  7  by  hla.  That  is,  Eq.  (1 1)  is  the  expected  backscattered  pressure  from  a 
bubble  of  radius  b  in  the  geometric  limit. 

C.  Scattering  from  a  fluid  sphere 

The  final  limiting  case  is  the  backscattering  from  an  ideal  fluid  sphere.^3d8-20  po^ 
a  fluid  sphere  the  inner  surface  has  vanished  (b  ->  0)  and  the  derivation  of  the  spreading 
factor  Hn  in  Appendix  A*  is  no  longer  applicable.  That  is,  the  rays  are  reflected  from  the 
opposite  side  of  tire  sphere  and  hence  from  a  surface  with  a  radius  of  curvature  of  -a.  To 
properly  determine  the  limiting  fonn  of  rhe  specular  reflection  form  function,  the  refractive 
powers  given  in  Eq.  (A3)  need  to  be  redefined  by  allowing  b  -a.  Alternately,  the 
appropriate  form  function  can  be  obtained  by  taking  the  limit  h  — >  in  Eq.  (7),  The  result 
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where  a  ~  2kih  =  Ak^a.  The  magnitude  of  the  contribution  of  the  nth  intemal  ray  to/\p  is 
!(1  -  r^)r2«-lyVf;,/  (2n  -Mul.  CompaiE  :on  of  Eq.  (34)  Ref.  13  with  the  magnitude  of 
tlie  nth  interna!  ray  here  slinws  that  for  the  fluid  sphere  limit  Eq.  (7)  reduces  to  tht 
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appropriate  form.  It  is  noted  that  the  rays  here  are  enumerated  consecutive  while  the 
numbering  scheme  in  Ref.  13  is  based  on  the  number  of  internal  chords  associated  with 
each  ray.  Finally,  it  can  be  shown  that  the  results  here  for  «  =  1  are  equivalent  to  those 
given  in  Appendix  C  of  Ref.  1 3  for  the  first  axial  ray. 

3.5  Conclusion 

When  the  incident  wave  is  at  a  frequency  of  a  resonance,  the  terms  of  the  series  in 
Eq.  (5)  are  mutually  in-phase:  there  is  constructive  interference  between  all  of  the  internally 
reflected  backward  directed  rays.  At  resonance,  the  externally  reflected  ray  (which  reflects 
from  S  in  Fig.  1)  interferes  either  destructively  or  constructively  with  die  summed 
contributions  of  the  internally  reflected  rays,  depending  on  the  sign  of  r.  For  the  usual  case 
where  r  >  0,  Eq.  (1)  gives  the  resonance  condition.  Unlike  Ike  series  for  sequential  surface 
wave  contributions,  which  may  be  expressed  in  a  Fabry-Perot  form,<^-'^d0  the  series  in  Eq. 
(5)  does  not  genet  ally  take  on  the  form  of  a  geometric  series.  If  the  radius  of  curvature  of 
the  shell  is  duly  negligible  in  comparison  to  the  thickness,  B  =  cihicb  0,  and  the  series 
becomes  a  geometric  series.  That  series,  when  summed,  gives  the  first  term  of  Eq.  (7) 
since  the  curvature  correction  miisr  v.mish  in  that  limit  The  nene.ral  form  of  the  cprie,':  in 

■  - - - - . - - - - -  - - ...i  . 

Eq.  (5)  is  not  that  of  a  geometric  series  since  the  spieading  factor  Hn  and  the  vfrtual  source 
location  depend  on  n. 

'fhe  curvature  of  the  shell  is  essential  to  the  manifestation  of  longitudinal  resonances 
in  the  modulus  of  the  backscattering  fonn  function  of  fluid  spherical  shells.  The  specular 
reflection  form  function  decompos.'.s  into  a  contribution  attributable  to  the  reflection  from  a 
flat  vacuum-backed  fluid  plate  and  a  curvature  conrection  component:  f^p  ~fp+  fee  "The 
dependence  of  fee  on  ki,h  manifests  the  resonance  condition  given  in  Eq.  (1).  Finally, 
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without  the  inclusion  of fee  in  the  synthesis  of  f^p,  the  approximation/ does  not 
exhibit  the  longitudinal  resonances  in  1/1. 

The  verification  of  the  geometrical  synthesis  oifsp  follows  from  comparison  of 
\fsp  I  with  the  exact  calculation  of  the  PWS  representation  of  1/  I.  The  excellent  agreement 
between  1/  I  and  I  fsp  I  demonstrates  the  correctness  of/vp  especially  for  j:  »  0. 
Furthermore,  the  limiting  cases  of  a  perfectly  soft  bubble,  fluid  sphere,  and  impedance- 
matched  fluid  shell  supports  the  claim  that  Eqs.  (7)  and  (8)  correctly  describe  the  high- 
frequency  specular  scattering  of  sound  from  an  ideal,  fluid  shell. 

A  synthesis  of  the  backscattering  of  an  elastic  stainless  steel  spherical  shell, 
summarized  in  Ref.  1 1 ,  neglected  the  curvature  correction /cc  since  the  expression  iotfee 
was  unavailable  at  the  time  of  that  work.  The  synthesis  used  the  method  described  in  Ref. 
6  for  the  conuibutions  from  leaky  Lamb  waves.  It  has  been  subsequently  demonstrated 
that  inclusion  of fee  improves  the  synthesis  and  it  is  our  intent  to  discuss  that  in  a 
subsequent  publication  since  the  synthesis  of  the  Lamb  wave  contributions  is  beyond  the 
scope  of  the  present  paper. 

Though  the  present  analysis  has  been  restricted  to  backscattering,  the  ray  diagram, 
Fig.  1,  gives  insight  into  how  such  rays  contnbute  to  the  angular  scattering  pattern  for  near 
backwards  directions.  The  virtual  sources  Vo»  —  li^^  along  the  backward  axis. 

Consider  first  the  case  of  only  two  such  sources.  The  interference  pattern  near  the  axis  is 
the  well  known  ring  or  bull's-eye  pattern  exhibited  in  an  optical  Michelson  interferometer 
or  Newton's  ring  experiment,^ The  center  of  the  ring  pattern  is  the  backwards 
direction.  This  pattern  was  previously  exhibited  in  a  ray  synthesis  of  backscattering  from 
fluid  spheres  in  cases  where  only  axial  ray  contributions  were  included.*^  (See  the  curves 
labeled  I /o  ■^■f2  I  in  Fig.  5  -  7  of  Ref.  I  'i.)  The  centr  al  scattering  amplitude  is  maximized 
when  the  outgoing  wavefronts  (which  are  locally  spherical  but  of  different  radii)  have  a 
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tangential  point  of  contact  on  the  backward  axis.  It  may  be  anticipated  that  tlie  situation 
depicted  in  Fig.  1  yields  ring  patterns,  but  with  a  richer  structure  due  to  the  additional 
virtual  sources  V2,  V3, ....  The  method  of  analysis  of  the  present  paper  and  these 
comments  on  the  scattering  pattern  should  also  be  applicable  to  understanding  the 
consequences  of  multiple  reflections  in  spherical  domes  which  enclose  sonar  transducers. 

While  the  analysis  given  here  has  been  limited  to  spherical  shells,  Fig.  1  also 
applies  to  a  right  circular  cylindrical  shell.  Extension  of  the  analysis  to  that  case  should 
give  the  curvature  correction  to  the  specular  contribution  previously  modeled  by  Borovikov 
and  Veksler.^ 

This  research  was  supported  by  the  Office  of  Naval  Research, 


Appendix  A.  Determination  of  Hn  and  Cr  from  geometrical  ray  methods 
The  factor  Hn  is  derived  by  applying  the  matrix  method  of  geometric  optics  to  the 
acoustic  ray  representation  of  the  wavefield  in  water  and  fluid  shell.  Each  factor  is 
associated  with  a  virtual  point  source  as  discussed  in  Sec.  3.2,  (see  Fig.l).  In  particular, 
the  ray  associated  with  the  local  outgoing  wavefront  produced  by  the  poim  source  Vn  will 
be  called  the  nth  ray.  The  local  backscattering  angle  y„  is  the  angle  between  the  nth  ray  and 
the  2-axis.  The  square  of  the  spreading  factor  in  terms  of  the  impact  parameter  s  and 
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where  [do  /  d^i]n  is  the  spreading  factor  part  of  a  differential  scattering  cross  section  for  the 
nth  ray.  This  is  a  specialization  to  y  0  of  the  more  general  result  (do/dO.)  -  s{sinyy 
Mds/d^  derivable  from  flux  conservation.^'!  The  purpose  of  this  appendix  is  to  outline  the 
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metliod  used  to  relate  the  impact  parameter  to  the  local  backscattering  ang’c  and,  ultimately, 
to  obtain  a  simple  analytic  expression  for  //«. 

The  matrix  methods  of  geometric  optics  are  employed  to  related  s  and 
Since  the  concern  here  is  backscattering,  then  the  paraxial  form  of  the  matrix  representation 
can  be  used.  Tlie  conventions  as  outlined  in  Refs.  21  and  22  and  summarized  below  are 
adopted  for  the  following  discussion.  The  relevant  matrices  are: 
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’  21 

where  X  is  the  matrix  defining  the  initial  angle  of  inclination  and  displacement  from  the  z- 
axis  of  the  incident  ray  at  point  B  in  Fig.  1.  The  final  angle  of  inclination  and 
displacement  x,i  of  the  nth  ray  are  represented  by  X„.  In  the  paraxial  approximation,  is 
the  ray  displacement  from  the  z-axis  in  the  vertical  plane  through  vertex  S.  The  matrix 
applies  to  the  point  where  the  ray  exits  the  shell.  The  refraction  and  translation  matrices  are 
Ri  and  T2i.  The  subscript  i  =  1,2  corresponds  to  the  surface  r  =  a,b.  For  paraxial  rays 
T21  =  Ti2  and  h  =  a{l  -  bldllMi  where  =  dci  is  the  relative  index  of  refraction. 

Noting  that  a  prime  on  R/  indicates  a  reflection  at  the  surface  i,  the  relevant  constants  are: 


1  -M, 


(A3) 


where  Tq  represents  the  special  case  of  the  zeroth  ray.  The  3^  are  referred  to  as  the 
refractive  powers  of  the  ith  surface.  The  trajectory  of  the  nth  ray  is  described  by  the  matrix 
equations. 
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where  J  =  and  M  =  Rj^R'^R^^ .  The  exponent  in  Eq.  (A4)  implies  a 


repeated  product,  i.e.,  JMJM  for  ra  =  3.  Performing  the  matrix  multiplications  leads  to  the 
result. 
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(n  >  0),  (A5) 


where 

B-  1-b/a  _  h 
M^(b/a)  M^b  * 


(A6) 


Equation  (A5)  contains  the  desired  relationship  between  s  and  y^.  The  analytic  fonn  of  the 
spreading  factor,  determined  by  using  Eqs.  (Al)  and  (A5),  gives 


u  -  (g/2) 
n  1  +  nB’ 


(n  =  0,  1,  ...).  (A7) 


From  inspection  of  Fig.  1,  in  the  paraxial  approximation  the  distance  of  the  virtual  soiu'ce 
Vn  from  Oisa-  (x«/y«).  From  Eqs.  (A5)  and  (A7),  we  find  that  a  -  (Xn/yn)  =  Hence, 
Vn  O  as  n  oo. 

The  amplitude  factor  C,.  in  Eq.  (3)  is  found  by  considering  the  partial  reflection 
and/or  transmission  of  a  ray  at  each  interface.  From  Fig.  1,  it  is  observed  that  the  zeroth 
ray  is  partially  reflected  at  r  =  a  without  further  interaction  with  an  interface.  The 
appropriate  constant  is  Cq  =  4  where  ris  the  reflection  coefficient  defined  in  Eq.  (4).  The 
ray  ABDEF  has  the  amplitude  factor  C;  =  tr2t’.  The  transmission  coefficients  t  and  t'  are 
for  rays  which  propagate  from  water  into  the  fluid  and  from  the  fluid  into  water, 
respectively.  The  transmission  coefficients  are  related  to  the  reflection  coefficient  in  Eq.  (4) 
by  t  =  1  -  rand  f  =  1  +  r.  The  reflection  coefficient  at  the  inner  surface  at  r  =  b  is  r2  =  -1. 
Using  these  relationships  gives  Cj  =  -(I  -  r^).  The  amplitude  constant  for  the  ray 
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ABDEGHI  is  C2  =  tr2  rr2t'  where  r  =  -r  is  the  reflection  coefficient  for  the  ray  incident  in 
the  fluid  on  the  fluid- water  interface.  Tlie  constant  C2  in  terms  of  r is  C2  =-r(l  -  r^).  In 
general,  the  amplitude  constant  for  the  nth  ray  is 

C„  =  -r"-^l-  r2)  (A8) 


Appendix  B.  The  partial  wave  series  representation  of  / 

The  PWS  representation  of  the  backscattered  form  function  is  included  for 
completeness.  This  PWS  is  for  the  backscattering  of  an  acoustic  plane  wave  from  an  ideal 
fluid  spherical  shell.  The  interior  of  the  shell  is  assumed  to  be  a  vacuum.  The  analysis  is 
similar  to  those  in  Refs.  3  and  19  and  gives/as 
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where  and  are  3  x  3  determinants  obtained  from  the  boundary  conditions  at  r  =  a  and 
b.  In  Appendix  B,  n  denotes  the  partial- wave  index.  The  elements  of  the  determinants  are: 


bjj  =  p  (c/c^)jn(x), 

dji  =  -p  (c/Cj^)  hn  \x). 

b2i=jn(x). 

d2^  =  -h^^\x). 

b3,=d3,=0. 

^12  ~  ^12  ~ 

b|3  =  tij3  =  njj(x^), 

^23  ^23 

*^33  ""  ^^33  "n(yL^- 

where  the  parameters  are  defined  as  follows:  p  =  p/Pc»  x  =  ka,  xl  =  xc/cl  and  yi  -  xibla. 
The  spherical  Bessel  and  Neumann  functions  and  Hankel  function  of  the  iii  ;t  kind  of 
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integer  order  are  denoted  hy  in,  nn,  and  hn,  respectively.  The  prime  on  a  function  indicates 
differentiation  with  rs;spect  to  the  argument. 

The  numerical  implementation  of  the  PWS  requires  that  the  infinite  series  must  be 
truncated.  A  discussion  of  a  choice  of  the  maximum  partial-wave  index  for  the  truncated 
series  is  given  in  Ref.  25.  The  result  is 


”max  =  ^  +  f  ^ 

'^max  =  5  +  [  X  +  4.05x^^'  j,  x  >  8, 


(B2) 


where  [  ]  implies  truncation  to  the  nearest  integer.  The  criterion  of  Eq.  (B2)  was  used  for 
the  PWS  compulations  shown  in  Figs.  2-4.  The  number  of  terms  exceeds  x  and  the  PWS 
converges  slowly  until  n>x. 


Appendix  C.  Rapid  summation  technique  for  fee 

Since  an  analytic  expression  is  unavailable  at  this  time  fox  fee,  then  the  infinite  series 
in  Eij.  (8b)  must  be  evaluated.  For  the  material  parameters  given  in  Table  I,  the  infinite 
series  convcigcs  fairly  fast,  but  for  other  values  the  series  may  not  converge  rapidly.  Here, 
we  investigate  a  method  that  may  increase  the  convergence  of  fee  =  where 


(l-r-^)y  nBrV”“ 

r  1  +  nB 

n  =  1 


(Cl) 


or  at  least  give  a  good  approximation  for  the  infinite  sum.  First,  it  is  noted  that  there  exists 
an  integer  n  =  N  such  that  Afi  >  1  for  any  B.  Now,  the  summation  in  Eq.  (Cl)  is  split  into 
tv/o  summations. 


n=:N+l 


(C2) 


where  S  denotes  the  sum  in  Eq.  (Cl)  and  G(N)  is  the  finite  sum  of  the  first  N  terms. 
Shifting  the  summation  index  in  the  second  term  of  Eq.  (C2)  and  dividing  the  numerator 
and  denominator  by  (n  +  N)B,  this  sum  is  rewritten  as 
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Since  [(n  +  <  1,  then  the  denominator  may  be  expanded  by  the  binomial  expansion. 

Neglecting  terms  of  order  [(n  +  ,  the  summation  on  the  right  hand  side  of  Eq.  (C3) 

is  approximated  by  two  infinite  series.  The  first  is  a  geometric  series  and  can  bt^  summed 
analytically;  while  the  second  is  the  first  correction  to  the  rapid  summation  technique. 
Equation  (C2)  becomes 


=  G(N)  + 


1  -  re‘“ 


-  r 


N JNa  r”e’”” 

(n  +  N)B  ‘ 


(C4) 


Finally,  the  rapid  summation  technique  consists  of  approximating  S  by  the  first  two  terms 
in  Eq,  (C4).  An  estimate  of  the  error  in  implementing  the  rapid  summation  teclinique  can 
be  ascertained  by  summing  the  last  term  in  Eq.  (C4). 

Numerical  tests  of  the  above  technique  indicate  that  suitable  converges  occurs  for 
the  choice  of  A  =  }  +  10  where  { }  implies  the  integer  part  of  B' V  (If  B-^  -  (B-l }  > 

0.5,  then  { }  is  increased  by  one.)  For  example,  for  a  fluid  shell  with  the  parameters  of 
440c  stainless  steel  and  bla  =  0.838,  the  constant  B  =  0.765  which  gives  {B'l }  =  1. 

Hence,  the  rapid  summation  technique  involves  the  calculation  of  i2  terms;  while  directly 
summing  5  requires  between  215  to  220  terms  for  convergence  to  six  decimal  places  of 
precision.  Figure  Cl  is  a  comparison  of  the  results  of  the  rapid  summation  technique  and 
tlie  truncated  infinite  series.  Agreement  between  the  two  curves  is  good  and  it  is  anticipated 
that  a  small  increase  in  N  will  give  even  better  agreement. 
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Chapter  4 

Ray  synthesis  of  the  form  function  for  backscattering  from  an  elastic 
spherical  shell:  Leaky  Lamb  waves  and  longitudinal  resonances. 

Abstract 

An  acoustic  ray  analysis  is  employed  in  synthesizing  the  form  function  for 
backscattering, /(0  =  n,  ka),  from  a  fluid-loaded  evacuated  elastic  spherical  shell  where  k  is 
the  wavenumber  of  the  incident  plane  wave  and  a  is  the  outer  radius  of  the  shell.  The 
synthesis  contains  a  component  associated  with  a  specular  reflection, and  contributions 
from  leaky  Lamb  waves.  The  contribution//  of  the  Ith  leaky  Lamb  wave  is  expressible  in  a 
Fabrj'-Perot  resonator  form  [P.  L.  Marston,  J.  Acoust.  Soc.  Am.  £2,  25  -  37,  (1988)]. 

The  present  synthesis  differs  from  previous  results  by  including  the  effects  of  longitudinal 
resonances  on  fsp-  A  novel  ray  synthesis  of fsp  indicates  a  significant  resonance  effect  near 
the  condition  kih  =  /jjt:  («  =  1,  2, ...).  The  thickness  of  the  shell  is  h  and  ki  =  co/c/  is  the 
longitudinal  wavenumber  where  ci  is  the  longitudinal  speed  of  sound  in  the  elastic 
material.  The  ray  synthesis  demonstrates  that  the  curvature  of  the  shell  is  essential  to  the 
modeling  of  longitudinal  resonances.  A  comparison  of  the  ray  synthesis  for  f{ka)  with  the 
exact  partial-wave  series  representation  for  a  440c  stainless  steel  shell  displays  the 
usefulness  of  the  ray  synthesis.  Although  acoustic  ray  modeling  is  generally  a  high- 
frequency  technique,  the  ray  synthesis  of f{ka)  for  a  440c  stainless  steel  shell  appears  to  be 
applicable  for  ka  as  small  as  7.  Certain  anomalies  in  the  synthesis  are  investigated  to  better 
understand  the  limitations  of  the  present  ray  model. 
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4.1  Introduction 

The  application  of  novel  ray  techniques  to  canonical  scattering  problems  allows  one 
to  test  the  integrity  and  efficiency  of  ray  representations.^'^  Since  ray  representations  of  the 
interaction  of  sound  with  a  scatterer  give  a  simple  picture  of  the  scattering  process,  they 
merit  investigation.  Ray  representations  have  the  potential  to  be  generalized  to  non- 
canonical  problems  not  easily  treated  by  other  methods.  Furthemiore,  the  need  for 
generalized  ray  theories  has  been  recognized  in  stmctural  acoustics.^  In  developing  such 
ray  techniques,  comparison  of  an  available  exact  solution  and  a  ray  model  for  a  canonical 
problem  gives  physical  insight  into  regions  where  the  techniques  are  applicable.  This 
article  investigates  a  novel  ray  synthesis  of  the  canonical  problem  of  high-frequency 
backscattering  from  a  fluid-loaded  evacuated  elastic  spherical  shell.  It  is  anticipated  that  the 
ray  synthesis  developed  here  may  be  generalized  to  scattering  from  other  smooth  convex 
elastic  objects.  Hence,  botli  the  success  and  failure  of  the  present  level  of  ray  synthesis 
may  be  instructive  for  more  advanced  ray  models. 

The  ray  synthesis  is  based  on  the  application  of  the  principles  of  the  generalized 
geometric  theory  of  diffraction,  It  has  been  proposed  and  demonstrated  that  for 

sufficiently  high  frequencies,  the  steady-state  form  function /for  backscattering  from 
elastic  spheres  can  be  partitioned  into  three  distinct  components.  These  contributions  are  an 
ordinary  specular  reflection,  transmitted  bulk  waves,  and  surface  guided  elastic  waves.  In 
the  development  of  a  ray  synthesis  of /for  an  elastic  shell,  the  ordinary  specular  reflection 
and  transmitted  bulk  wave  contributions  are  grouped  to  form  a  generalized  specular 
reflection  component  denoted  by  f^p-  The  present  ray  synthesis  oifsp  differs  from  a 
previous  my  rnodel^'*  ^  by  introducing  a  curvature  coirection,/cc.  Inclusion  of fcx  in  the 
synthesis  shows  that/^p  can  be  significantly  affected  near  longitudinal  resonances  within 
the  shell.  The  conditions  for  excitation  of  longitudinal  resonances  are 


(la) 

(lb) 
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=  nTT,  (n  =  1,  2,  „.), 

k^h  =  (n  +  1/2)7c,  (n  =  0, 1,  -.), 


where  h  =  (a-  b)is  the  thickness  of  the  shell.  Equation  (la)  corresponds  to  the  usual  case 
(the  one  considered  here)  where  the  acoustic  impedance  of  the  elastic  material  is 
greater  than  the  impedance  of  the  surrounding  water  pc.  When  p^c^  is  less  than  pc  then 
Eq.  (lb)  is  the  appropriate  condition.  These  resonance  conditions  are  obtained  from  the 
requirement  of  constructive  interference  of  consecutive  internally  reflected  rays  upon 
transmission  back  into  the  water.  Finally,  for  high-frequency  scattering  from  the  elastic 
shell  considered  here,  the  relevant  surface  guided  elastic  waves  are  leaky  Lamb  waves. 

The  steady-state  scattered  pressure  in  the  farfield  from  an  evacuated  elastic  spherical 
shell  has  the  fonn^^ 


n  =  n  ■Sf. 

Ps  Pq  2r  ’ 


(2) 


where  pQ  is  the  pressure  amplitude  of  the  incident  acoustic  plane  wave.  The  wavenumber  k 
of  the  incident  plane  wave  is  defined  by  ^  =  co/c  where  to  is  the  angular  frequency  and  c  is 
the  speed  of  sound  in  water.  The  harmonic  time  dependence  exp{-mt)  has  been 
suppressed.  The  outer  radius  of  the  spherical  shell  is  a  and  r  is  the  distance  from  the  center 
of  the  shell  to  some  distant  observation  point.  The  complex  scattering  amplitude  or  form 
function  in  the  backscattered  direction  has  the  exact  partial-wave  series  representation^^ 


f  =  ^£(-l)"(2n.-l) 


Bn(^) 


(3) 


where  x  =  ka-  Inafk  and  X  is  the  wavelength  of  the  incident  plane  wave.  The  functions 
B^ix)  and  Dn{x)  are  5  x  5  determinants  obtained  by  satisfying  the  appropriate  boundary 
conditions. ’2  The  elements  of  these  determinants,  which  are  complicated  expressions  of 
spherical  Bessel  functions  and  spherical  Hankel  functions  of  the  first  kind,  are  listed  in 


95 


Ref.  12.  The  material  parameters  for  the  elastic  shell  are  the  longitudinal  sound  speed  ci, 
the  shear  or  transverse  sound  speed  Cs  and  the  density  pg  while  the  density  of  water  is 
denoted  by  p. 

The  organization  of  this  article  is  as  follows.  In  Sec.  4.2.A,  the  ray  synthesis  of  an 
individual  leaky  Lamb  wave  contribution,//,  to  the  form  function  for  backscattering  is 
developed.  The  expression  for//  is  cast  in  a  fomi  analogous  to  a  Fabry-Perot  resonator. 
Section  4.2.B  summarizes  a  ray  synthesis  of/jp  and  demonstrates  that  fsp  can  be  separated 
into  a  term  associated  with  reflection  from  a  vacuum-backed  f  at  elastic  plate  and  a 
curvature  dependent  correction.  Using  parameters  corresponding  to  a  440c  stainless  steel 
shell  studied  in  Ref.  12,  the  results  of  the  ray  syntnesis  and  exact  partial- wave  scries  are 
compared  in  Secs.  4.3  and  4.4.  Since  a  longitudinal  resonance  does  not  significantly  affect 
the  ranges  0  <  jc  <  60  and  80  <  x  <  100,  the  computations  in  Sec.  4.3  consider  these 
regions.  The  presence  of  a  longitudinal  resonance  in  the  range  60  <  x  <  80  complicates  the 
analysis  and  this  region  is  investigated  in  more  detail  in  Sec.4.4.  A  general  discussion  and 
concluding  remarks  are  contained  in  section  4.5.  Appendix  A  contains  a  discussion  of  the 
radiation  damping  parameter  and  normalized  phase  velocity  for  each  leaky  Lamb  wave 
contribution.  These  leaky  Lamb  wave  parameters  are  essential  to  the  ray  synthesis 
calculations  and  are  included  for  completeness.  The  first  antisymmetric  and  syrmnetric 
leaky  Lamb  wave  contributions  are  isolated  in  Appendix  B.  The  Fabry-Perot  form  of//  is 
investigated  and  a  discussion  of  off-resonance  contributions  from  a  Lamb  wave  to  the  form 
function  for  backscattering  is  presented.  Finally,  Appendix  C  studies  a.  peak  near  a:  =  71  in 
the  ray  synthesis  of/ which  is  not  contained  in  the  exact  partial-wave  series  result.  The 
discussion  in  Appendix  C  is  based  in  part  on  a  localization  principle  analysis  of  ir  dividual 
partial  waves  similar  to  that  used  by  van  de  Hulst  for  light  scattering. 
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4.2.  Ray  synthesis  of  the  form  function  for  backscattering 
A.  Leaky  Lamb  wave  contributions 

The  methodology  of  the  generalized  geometrical  theory  of  diffraction  is  applied 
to  the  backscattering  of  a  plane  wave  from  a  hollow  elastic  spherical  shell.  For  the 
canonical  problem  of  high-frequency  scattering  from  elastic  shells,  leaky  Lamb  waves  are 
the  relevant  surface  guided  elastic  waves  which  produce  the  resonance  structure  observed 
in  the  form  function  for  backscattering.  (In  Sec.  4.3  below,  it  is  demonstrated  that  the  ray 
synthesis  works  well  down  xox~l  for  the  particular  shell  considered.  For  smaller  x  the 
contributions  of  subsonic^^  Lamb  and  Franz-type  waves  may  be  significant.)  The  ray 
diagram  in  Fig.  1  provides  a  fundamental  illustration  of  the  ray  synthesis  of //.  This 
diagram  is  similar  to  those  discussed  in  Refs.  1, 6, 12,  and  15.  The  diagram  represents  an 
elastic  spherical  shell  of  outer  radius  a  and  inner  radius  b  situated  such  that  the  shell's 
center  C  coincides  with  the  origin  of  the  coordinate  axis.  A  plane  wave  propagates  in  the 
positive  z-direction  along  the  ray  AB  in  Fig.  1.  At  points  the  acoustic  wavefield  in  water 
couples  to  the  shell  and  launches  a  leaky  Lamb  wave.  The  Lamb  wave  is  guided  by  the 
shell  along  the  arc  BB'  continuously  re-radiating  energy  into  the  water.  At  B',  the  radiated 
energy  is  backscattered  along  the  ray  B’A\  (Physically,  B  and  B’  actually  locate  surface 
regions  where  the  interaction  takes  place.  The  widtlt  of  the  regions  correspond  to  the  size 
of  Fresnel  zones.  The  leaky  Lamb  wave  repeatedly  circumnavigates  tlie  shell  radiating 

energy  in  the  backscattering  direction  with  each  circumnavigation.  The  points  B  and  are 
obtained  from  the  phase  velocity  trace-matching  condition, 

0^  =  arcsin(c/cp,  c^  S:  c,  ^4^ 

where  6/  is  the  local  angle  of  incidence  and  c/  is  the  phase  velocity  of  the  leaky  Lamb  wave, 
■^fhe  phase  velocity  is  taken  to  be  the  phase  velocity  of  the  leaky  Lamb  v.'ave  along  the  outer 
surface  of  the  shelP^  and  is  assumed  to  be  supersonic. 
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Fig.  1  A  ray  diagram  for  a  contribution  from  a  leaky  Lamb  wave  to  the  form  function 
for  backscattering  from  ?»n  elastic  spherical  shell  of  outer  radius  a  and  inner  radius  b.  A 
plane  wave  represented  by  the  ray  AB  couples  to  shell  at  B  launching  a  leaky  Lamb  wave. 
The  Lamb  wave  is  guided  by  the  shell  through  a  polar  caustic  at  C"  to  B’.  At  B'  the 
radiated  energy  is  shed  in  the  backscattered  direction  along  the  ray  B'A'.  'fhe  points  B  and 
B'  are  determined  from  Eq.  (4)  where  0/  is  the  local  angle  of  incidence.  The  Lamb  wave 
continues  to  circumnavigate  the  shell  shedding  some  energy  at  fl'  with  each  passage.  The 
point  C'is  the  specular  point  and  the  location  of  a  second  polar  caustic.  For  backwards 
and  near  backwards  scattering  the  rays  fi'A'  and  DFi  appear  to  diverge  from  a  virtual 
source  located  at  point  F/.  When  the  diagram  is  rotated  about  C'C,  the  point  F/  traces  out  a 
virtual  ring-like  source.  The  vertical  line  through  C  represents  an  exit  plane  for  an 
outgoing  toroidal  wavefront  associated  with  the  ring-like  source. 
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The  form  of  an  individual  leaky  wave  contribution  to  the  backscattering  amplitude  is 
obtained  by  summing  the  backscattering  amplitudes  for  each  circumnavigation,  1*7 


f^  =  -G^ e'^^ "  ■  J  ^i2nmxc/c 


m  =  0 


The  geometrical  phase  shift  r[i  =  [2x(c/c/)(rc  -  6/)  -  2xcos(Qi)  -  k/2]  accounts  for  the  phase- 
difference  between  a  Lamb  wave  propagating  along  the  arc  BB'  and  a  wave  in  water 
propagating  from  C'  to  C  and  back  to  C'  as  if  the  shell  were  not  present.  The  7t/2  term  in 
r\l  is  associated  with  the  propagation  of  the  leaky  Lamb  wave  through  the  polar  caustic  at 
C".  The  additional  phase  shift,  IkkcIci,  is  attributable  to  repeated  circumnavigation  of  the 
leaky  Lamb  wave  about  the  shell.  The  tc  phase  shift  results  from  the  polar  caustics  at  C 
and  C".  The  radiation  damping  of  the  leaky  Lamb  wave  is  characterized  by  the  radiation 
damping  parameter  p/  with  units  of  Npirad.  The  factor  exp[-2(jc  •  0/  )P/]  is  the  radiation 
damping  of  the  partial  circumnavigation  of  the  Lamb  wave  along  BB';  while  expi-2K^i) 
accounts  for  the  additional  radiation  damping  due  to  repeated  circumnavigation.  TTie 
efficiency  of  the  coupling  between  the  wavefield  in  water  and  the  leaky  Lamb  wave  at  the 
points  B  and  B'  is  characterized  by  a  complex  coupling  coefficient  G/.  Inspection  of  the 
summation  in  Eq.  (5)  shows  that  it  is  a  geometric  series  and  Eq.  (5)  reduces  to 


7 


ill,  -2(7t-e,)B,  i2nxc/c, 

e  'e  ^  ^]/[  1  +  e  U  ^J. 


(6) 


The  form  of  the  denominator  in  Eq.  (6)  is  that  of  a  Fabry-Perot  resonator,  Appendix  B 
contains  a  discussion  of  1// 1  for  the  lowest  antisymmetric  and  symmetric  leaky  Lamb 
waves.  Finally,  Eq.  (6)  is  the  ray  synthesis  of  the  contribution  of  an  individual  leaky  Lamb 
wave  to  the  form  function  for  backscattering. 

Williams  and  Marston^  applied  the  Sommerfeld-Watson  transformation  to  the 
partial-wave  series  representation  of / for  backscattering  from  a  solid  elastic  sphere.  From 
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their  analysis  a  complicated  expression  for  the  complex  coupling  coefficient  Gi  was 
obtained  (see  Eq.  (28)  of  Ref.  6).  Subsequently,  Marston^  developed  an  approximation 
for  Gi  which  gives  explicit  dependence  of  I  G/ 1  on  the  leaky  Lamb  wave  parameters  p/  and 
cilc.  The  approximation  can  be  expressed  as 

G^  =  87tp^(c/cpe‘‘'’^  (7) 

where  0/  denotes  the  phase  of  G/  and  it  may  be  argued  that  0/  =  0  for  elastic  spheres  or 
shells.lLlS  Pqi-  the  computations  of  the  ray  synthesis  of  I/I  in  Secs.  4.3  and  4.4, 0/  for 
the  Ith  leaky  Lamb  wave  is  taken  to  be  identically  equal  to  zero.  Equation  (7)  was 
employed  in  the  modeling  of  backwards  and  near  backwards  scattering  of  short  tone  bursts 
from  a  440c  stainless  steel  shell  in  experiments  discussed  in  Ref.  12.  Further  confirmation 
of  the  applicability  of  Eq.  (7)  to  the  backscattering  from  elastic  spherical  shells  is 
demonstrated  below  in  Secs.  4.3  and  4.4.  An  alternative  derivation  of  Eq,  (7)  with  0/  =  0 
has  been  given  for  a  thin  fluid-loaded  spherical  shell.^ 

Consider  again  the  ray  diagram  in  Fig.  1.  For  an  observation  point  on  the  negative 
z-axis  a  finite  distance  from  C ,  the  ray  that  reaches  the  observation  point  is  DF/  not  ray 
B'A'.  It  can  be  argued  that  the  appropriate  modifications  to  the  ray  synthesis  of fi  are 
negligible  provided  the  observation  point  is  sufficiently  distant.^  The  intersection  of  rays 
B'A'  and  DFi  at  the  point  F/  locates  a  virtual  source.  The  point  F/  is  situated  in  the  vertical 
plane  perpendicular  to  the  z-axis  through  the  point  From  the  symmetry  of  the  shell,  it 
is  apparent  that  the  virtual  source  is  ring-like  and  has  a  radius  bi  =  a  siniQi).  Hence,  the 
backwards  and  near  backwards  contribution  of  a  leaky  Lamb  wave  may  be  modeled  by  a 
virtual  ring-like  source.  The  local  outgoing  wavefront  in  the  vertical  plane  through  point  C 
from  a  ring-like  source  has  a  toroidtil  shape.  A  toroidal  wavefront  is  associated  with  the 
phenomena  of  axial  focusing  or  glory  scattering.  ^^^9-22  Axial  focusing  of  the 
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backscattering  from  solid  elastic  spheres^^*^®  and  elastic  spherical  shells^^  has  been 
observed. 

B.  Specular  reflection  contribution 

The  specular  reflection  contribution in  the  form  function  for  backscattering 
from  an  evacuated  elastic  spherical  shell  can  be  approximated  by  ray  methods.  Figure  2 
shows  a  simplified  ray  diagram  which  facilitates  the  ray  synthesis.  The  diagram  is  similar 
to  the  one  discussed  in  Ref.  23.  The  orientation  of  the  shell  with  respect  to  the  incident 
plane  wave  is  the  same  as  in  Fig.  1  where  the  points  S  and  0  in  Fig.  2  correspond  to  C 
and  C  in  Fig.  1,  respectively.  The  specular  reflection  appears  to  originate  at  S  and  is  often 
referred  to  as  the  specular  point.24  The  incident  rays  which  significantly  affect fsp  lie  close 
to  the  z-axis.  Consider  the  incident  ray  AB  in  Fig.  2  with  infinitesimal  impact  parameter  s. 
Rotation  of  the  figure  about  the  line  SO  generates  all  the  incident  rays  with  impact 
parameter  s.  At  point  B  the  ray  is  partially  reflected  back  into  the  water  and  partially 
transmitted  into  the  elastic  shell.  The  intersection  of  the  reflected  ray  BC  with  the  axis  SO 
at  Fq  defines  the  location  of  a  virtual  source.  For  backwards  and  near  backwards 
scattering,  spherical  aberration  may  be  neglected  and  the  virtual  source  can  be  described  as 
a  point  source.  Since  the  impact  parameter  s  is  assumed  small,  then  the  angle  between  the 
incident  ray  and  the  outward  unit  normal  at  point  B  is  infinitesimal.  Hence,the  ray  AB  is 
approximately  at  normal  incidence  with  the  shell  and  the  energy  transmitted  into  the  shell 
will  produce  a  (bulk)  longitudinal  wave.  To  simplify  the  ray  synthesis,  the  weak  mode 
conversion  between  a  longitudinal  wave  and  a  transverse  (shear)  wave  upon  reflection  at  a 
surface  is  neglected  since  the  angle  of  incidence  is  infinitesimal.  The  transmitted  ray  at  B  is 
refracted  and  propagates  along  BD.  At  D  (again  neglecting  mode  conversion),  total  internal 
reflection  occurs  and  the  longitudinal  wave  propagates  along  DE.  ITie  ray  is  again  partially 
reflected  and  transmitted  at  E.  Projection  of  the  transmitted  ray  after  refraction  at  E  back  to 
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Fig.  2  The  ray  diagram  for  the  specular  reflection  from  an  evacuated  elastic  spherical 
shell  of  outer  radius  a  and  inner  radius  b.  The  point  S  is  the  vertex  of  the  refracting  surface 
and  O  is  the  origin  of  a  coordinate  axis  located  at  the  center  of  the  shell.  The  ray  AB 
infinitesimally  close  to  the  z-axis  is  incident  on  the  shell  with  impact  parameter  s.  The  ray 
ABC  is  the  ordinary  specular  ray;  while  ABDEF  and  ABDEGHI  are  the  first  two  internal 
specular  reflection  contributions.  Air  =  b,  the  rays  arc  totally  reflected  and  the  rays  are 
partially  reflected  and/or  transmitted  at  the  water-shell  interface.  Intersection  of  the 
projection  (dashed  lines)  of  the  outgoing  rays  and  the  z-axis  define  locations  of  virtual  point 
sources,  which  de.scribe  local  curvature  of  the  wavefront  associated  with  each  ray.  The 
specular  reflection  contribution  to  the  fonn  function  for  backscattering  is  determined  from  a 
superposition  of  the  wavefields  from  the  virtual  sources. 
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the  z-axis  gives  the  location  of  a  second  virtual  point  source  Vi.  By  continuing  to  trace  ihe 
ray  trajectories  of  the  higher  internal  reflections  an  infinite  set  of  virtual  point  source 
locations  Vn  on  the  z-axis  can  be  generated. 

Each  virtual  point  source  Vn  describes  the  propagation  of  a  local  outgoing 
wavefront  associated  with  the  nth  ray.  In  the  limit  s  — >  0,  the  specular  reflection 
contribution  to  the  farfield  pressure  is  a  superposition  of  the  acoustic  wavefields  from  the 
virtual  point  sources. 
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(8) 


where  ki  =  (nfci.  The  2jc  phase  shift  accounts  for  the  path  length  difference  between  a  ray 

propagating  in  water  from  S  to  O  and  back  to  S  and  a  ray  which  is  backscattered  at  S.  The 

additional  phase  delay  within  the  summation  comes  from  a  longitudinal  wave  which 

reverberates  n  times  within  the  shell  before  being  transmitted  in  the  backscattered  direction. 

The  amplitude  of  the  nth  virtual  point  source  is  HnCn  where  Hn  represents  a  geometric 

divergence  factor  and  Cn  accounts  for  the  cumulative  effects  of  partial  reflection  and/or 

transmission.  By  applying  an  acoustic  ray  analysis  of  longitudinal  resonances  within  a 

fluid  shell  the  foUowine  exnressions  for  H«  and  C»  were  obtained23 
^  *  0*  •• 


and 


(a/2) 

1  +nB’ 


“  r, 

=  -r"-Vl-r2), 


(n  >  0), 


(9) 


(10) 


where  B  =  {cihlcb)  and  r  =  [pgC'L  -  pc]/[pecx  +  pc].  The  constant  B  depends  on  only  the 
normalized  longitudinal  sound  speed  of  the  elastic  material  and  tlie  radii  of  curvature  of  the 
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inner  and  outer  surfaces.  The  reflection  coefficient  for  the  scattering  of  an  acoustic  plane 
wave  at  normal  incidence  from  an  interface  between  water  and  an  elastic  half-space  is  r. 

After  substitution  of  Eqs.  (9)  and  (10)  into  Eq.  (8),  comparison  with  Eq.  (2)  gives 
an  expression  for  fsp  which  depends  on  the  geometric  and  material  parameters  of  the  shell, 


^sp 


(1  -  rh 


i2nk,  h 


E-*' ..  L 
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n  =  1 


.-i2x 


(11) 


Although  Eq.  (1 1 )  is  a  ray  synthesis  of  the  contribution  of  the  specular  reflection  to  the 
form  function  for  backscattering,  it  can  be  manipulated  such  that  each  resulting  term  has  a 
simple  physical  interpretation.  With  the  replacemeni  of  [l+«S]‘l  by  1  -  n5[l+nS]'l,  the 
summation  in  Eq.  (11)  can  be  split  into  two  summations.  The  first  summation,  via  a 
geometric  series,  has  the  analytic  form  Ttxp(J2kLh)[\-Ttxp(,i2ki,h)Y^ .  Hence, becomes 

^sp  =  ^cc’  (12) 

where 
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Define,^  =  ‘Iiexp{-i2x)  where  :^is  the  complex  reflection  coefficienr  associated  with  the 
scattering  of  a  normal  incidence  plane  wave  from  a  vacuum-backed  flat  elastic  plate  of 
thickness  )j.25,26  Hence, describes  the  reflection  from  an  elastic  spherical  shell  where 
the  complex  reflectivity  is  mcdeicd  as  that  for  a  flat  elastic  plate.  Consequently,  the 
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influence  of  the  curvature  of  the  shell  on  fsp  within  the  paraxial  approximation  must  be 
contained  in  Eq.  (14).  At  present  a  simple  analytic  expression  for  fee  is  unavailable,  but  a 
rapid  summation  teclinique  has  been  described  and  implemented  in  Appendix  C  of  Ref.  23. 

Figure  3  demonstrates  the  importance  of  fee  to  the  ray  synthesis  of  \fsp\.  The 
material  and  geometric  parameters  used  in  the  computation  correspond  to  the  parameters  of 
a  440c  stainless  steel  shell  studied  in  Ref.  12.  The  values  of  these  parameters  are  as 
follow:  Cl  =  5.854  km/s,  Cs  =  3. 150  km/s,  pg  =  7.84  g/cm^  and  bla  =  0.838  for  the  elastic 
shell;  c  -  1.479  km/s  and  p  =  1.00  gterrfi  for  water.  The  resonance  condition  given  in  Eq. 
(la)  predicts  that  the  a  longitudinal  resonance  should  occur  at  =  76.8.  Clearly,  the 
minimum  observed  in  Fig.  3  at  x  77  is  a  manifestation  of  this  longimdinal  resonance.  It 
is  noted  that  the  approximation  f^p  ~fpis  not  sufficient  to  produce  the  observed  minimum. 
It  is  relatively  simple  to  demonstrate  that  1)^  I  =  1  for  all  x.  The  importance  of  fee  m  the 
synthesis  of  i/I  for  a  440c  stainless  steel  shell  will  be  examined  in  more  detailed  below  in 
Sec.  4.4.  Hence,  the  curvature  of  the  shell  appears  to  be  intimately  related  to  the  presence 
of  a  longitudinal  resonance  in  the  specular  contribution  to  the  form  function  for 
backscattering. 

4.3.  Comparison  of  the  ray  synthesis  and  exact  calculation  of  |  / 1  for  jc 
outside  the  region  of  a  longitudinal  resonance 

To  facilitate  a  discussion  of  a  comparison  between  the  ray  synthesis  and  the  exact 
partial-wave  series  computation  of  I/I,  the  range  0  <  x  <  100  is  split  into  five  intervals.  By 
displaying  these  smaller  regions,  the  finer  details  of  the  high  Q  resonance  structure  caused 
by  some  of  the  leaky  Lamb  wave  contributions  are  more  easily  resolved.  The  range  60  <  x 
<  80  is  examined  in  more  detail  in  Sec.  4.4  because  the  presence  of  a  longitudinal 


Fig.  3  The  modulus  of  specular  reflection  fsp  from  an  elastic  spherical  shell  as  defined 
by  Eqs.  (12)  -  (14).  The  broad  minimum  at  jt  =  77  is  a  manifestation  of  a  longitudinal 
resonance.  The  resonance  condition  in  Eq.  (la)  predicts  a  longitudinal  resonance  occurs  at 
^LR  ~  76.8.  As  discussed  in  4.2.B,  if  the  curvature  correction  component/cc  were  omitted 
in  Eq.  (12),  then  l/^p  I  =■-  I/p  I  =  1  for  all  x. 
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resonance  obfuscates  the  cause  of  certain  resonance  features.  In  the  calculations  below, 
unless  otherwise  stated,  the  ray  synthesis  is 

fray  ^sp  5^  f/ =  fp  +  fee  +  5],  V  (15) 

/  / 

where  the  summation  in  Eq.  (15)  is  over  the  possible  leaky  Lamb  wave  contributions.  The 
material  of  the  elastic  shell  is  440c  stairiless  steel  with  the  parameters  listed  above  in  Sec. 
4.2.B. 

The  computation  of  the  partial-wave  series  for  the  form  function  for  backscattering 
requires  the  truncation  of  the  infinite  summation  in  Eq.  (3).  The  minimum  number  of  terms 
(or  maximum  partial-wave  index)  nmax  retained  for  sufficient  convergence  of^pws  exceeds 
X.  The  following  criterion  was  tested  and  determined  to  ensure  adequate  convergence  of 
the  series  for  the  shell  considered: 

Unjax  =  2  +  [  X -t- 4.0x^^^  ],  x<8, 

where  the  square  brackets  imply  truncation  to  the  nearest  integer.  For  very  thin  shells 
caution  in  the  use  of  Eq.  (16)  for  determining  rimax  naay  be  required  since  subsonic  guided 
wave  contributions  can  be  significant.  From  the  localization  principle  subsonic  guided 
waves  are  associated  with  partial  waves  having  n>x.  Discussions  concerning  the  choice 
of  Eq.  (16)  are  contained  in  Refs.  15  and  23.  Secondly,  to  ensure  adequate  resolution  of 
the  high  Q  resonance  structure,  the  increment  Ax  needs  to  be  small.  For  the  computations 
of  lypWS  I  shown  below  Ax  =  0.01  except  in  the  vicinity  of  the  very  narrow  resonance 
structure  caused  by  the  jq  leaky  Lamb  wave  (see  Appendix  B).  Near  these  resonances,  it 
was  necessary  to  take  steps  in  Ax  as  small  as  0.00001.  Hence,  it  is  evident  from  Eq.  (16) 
and  the  small  size  of  Ax  that  the  partial-v/ave  representation  of/in  high-frequency 
backscattering  calculations  can  be  numerically  intensive. 
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The  rclev£int  leaky  Lamb  waves  are  analogous  to  Lamb  waves  that  can  be  excited  on 
a  flat  elastic  plate  of  thickness  h  in  vacuum  27.28  Using  the  tenninology  associated  with  the 
flat  plate  results,  the  leaky  Lamb  waves  excited  on  the  shell  for  0  <  x  <  100  are  designated 
ao  and  ai  for  antisymmetric  or  flexural  Lamb  waves  and  5o.  ■si.  and  S2  for  symmetric  or 
dilatational  Lamb  waves.  As  the  frequency  increases,  higher  Lamb  wave  resonances  can 
couple  with  the  acoustic  wavefield  in  water  and  these  higher  modes  can  easily  be  included 
in  the  synthesis.  The  computation  of  Eq.  (15)  requires  the  parameters  p/(x)  and  c/(x)/c  for 
the  Ith  leaky  Lamb  wave.  These  parameters  can  be  determined  by  extending  certain  results 
of  the  Sommerfeld- Watson  transfomiation  analysis  of  the  backscattering  from  a  solid 
elastic  sphere.^'7  The  pertinent  results  from  this  analysis  are  surrimarized  in  Appendix  A 
and  p/(x)  and  c/(jc)/c  for  the  Ith  leaky  Lamb  wave  are  examined.  The  values  of  p/(x)  and 
ci0c)/c  depicted  in  Figs.  Al  and  A2  are  employed  in  the  present  calculations.  Finally,  Eqs. 
(6)  and  (12)  -  (14)  are  the  expressions  necessary  for  the  calculation  of  I/ray  I. 

Figure  4  displays  the  comparison  of  I  /pws  •  I  /^y  I  where  the  dashed  curve 
corresponds  to  the  exact  partial-wave  representation  and  the  solid  line  is  the  ray  synthesis 
result.  In  the  interval  0  <x  <  20,  Fig.  4a,  the  shell  can  suppon  both  the  ao  and  leaky 
Lamb  waves.  The  ray  synthesis,  Eq.  (15),  is  truncated  at  x  =  7  since  the  phase  velocity  for 
the  ao  leaky  Lamb  wave  becomes  subsonic.  From  inspection  of  Eq.  (4),  the  present  ray 
synthesis  is  implicitly  a  supersonic  theoiy.  Comparing  Figs.  4a,  Bl,  and  B2  demonstrates 
that  the  broad  resonance  structure  is  attributable  to  the  ao  mode;  while  the  Lamb  wave  is 
resjjonsible  for  the  sharp  (high  Q)  resonance  features  atx  =  10, 14,  and  18.  It  is 
noteworthy  that  the  high-frequency  ray  synthesis  gives  fairly  good  agreement  even  though 
X  is  not  very  large.  To  accurately  model  the  exact  l./pws  I  for  x  <  7,  both  the  subsonic 
portion  of  the  ao  Lamb  wave  and  Franz-type  (creeping)  waves  would  need  to  be 
considered.  Figure  4b  shows  the  range  20  <  x  <  40  where  the  relevant  leaky  Lamb  wave 
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contributions  are  oq  and  sq*  Again,  the  broad  resonance  structure  is  due  to  oq  and  the  sharp 
features  are  associated  with  the  Lamb  wave.  The  agreement  between  the  synthesis  and 
exact  solution  is  excellent.  Finally,  Fig.  4c  compares  I  jlpws  I  and  I /ray  I  where  the  next 
possible  leaky  Lamb  wave  can  contribute  to  the  backscattering.  This  new  contribution  is 
the  ai  leaky  Lamb  wave  which  appears  to  have  a  cutoff  frequency  near  ;t  =  41.  Inclusion 
of  the  OQ.  ■^Ov  and  a\  leaky  Lamb  waves  not  only  reproduces  all  the  resonance  features,  but 
the  agreement  between  the  ray  synthesis  and  the  exact  partial-wave  series  representation  is 
excellent 

The  range  80  <  jc  <  100  is  shown  in  Fig.5.  The  relevant  contributions  to  the  ray 
synthesis  contained  in  Fig.  5a  are  oq,  ai,  5o.  and  as  well  as  the  specular  reflection. 
Although  the  ray  synthesis  seems  to  reproduce  the  resonance  features  in  I  ^ws  some 
differences  are  evident.  Forx  >  83  the  ray  synthesis  can  be  improved  by  acluding  the  S2 
leaky  Lamb  wave  which  appears  to  have  a  cutoff  frequency  near  a:  =  83.  By  including  the 
S2  leaky  Lamb  wave  contribution  in  the  synthesis  of  Fig.  5b,  the  exact  result  is  modeled 
more  accurately  by  i /ray  I  particularly  for  jc  >  88.  The  difference  between  I/ray  I  and  l/>ws  I 
for  .r  <  85  might  be  associated  with  a  longitudinal  resonance  discussed  in  Sec.  4.4.  The 
broad  minimum  in  the  soecular  reflection  contribution  shown  in  Fie.  3  caused  bv  a 

A  w'  y 

longitudinal  resonance  dXxui  ~  76.8  may  significantly  affect  the  form  function  for  a:  <  85. 
Inspection  of  Fig.  5  denionstrates,  however,  that  the  geometric  ray  synthesis  contains  the 
resonance  structure  caused  by  various  leaky  Lamb  wave  contributions. 

4.4.  Form  function  in  the  vicinity  of  a  longitudinal  resonance. 

The  resonance  stmcture  in  I />ws  I  for  the  interval  60  <  a  <  80  is  perhaps  the 
most  difficult  region  to  interpret  physically.  Within  this  range,  a  longitudinal  resonance  can 
be 
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Fig.  4  Trie  dashed  line  is  the  exact  partiai-wave  representation,  i^ws  i.  and  the  solid 
line  is  the  ray  synthesis,  I  /ray  I  for  the  stainless  steel  shell  considered  with  b/a  =  0.838. 

The  pertinent  leaky  Lamb  waves  included  in  each  synthesis  are:  (a)  flo  and  jq!  (W  ^  and 
5();  (c)  aQ,  so,  and  a\.  In  (a)  the  ray  synthesis  is  truncated  below  x  =  7  since  ao  becomes 
subsonic  and  the  ray  synthesis  is  implicitly  a  supersonic  theory.  The  agreement  between 
the  high-freriuency  ray  synthesis  and  the  exact  results  indicate  that  the  ray  synthesis  may  be 


useful  for  non-canonical  scattering  problems. 
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Fig.  5  The  exact  I  /pws  I  I  /ray  I  are  the  dashed  and  solid  cuives,  respectively.  In 
(a)  the  relevant  leaky  Lamb  wave  contributions  arc;  ag,  aj,  5o>  and  while  in  (b)  the 
leaky  Lamb  waves  are:  aQ,a\,  sg.  i.  and  S2-  The  specular  reflection  employed  in  I  /ray  1  is 
given  by  Eqs.  vl2)  -  (14).  The  importance  of  these  figures  is  that  the  ray  synthesis 
accurately  models  the  resonance  saaicture  in  the  exact  result  and  the  inclusion  of  the  si 
leaky  Lamb  wave  in  (b)  improves  the  synthesis. 
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supported  by  the  elastic  shell.  From  Eq.  (la),  the  longitudinal  resonance  is  predicted  to 
occur  atx/j?  ~  76.8.  It  is  observed  in  Fig.  3  that  the  longitudinal  resonance  can  affect  a 
relatively  broad  range  of  x  in  the  vicinity  of  x/jr.  Also,  within  this  region  of  x  the  440c 
stainless  steel  shell  can  support  the  oq,  ai,  sq,  and  si  leaky  Lamb  waves.  Of  these  Lamb 
waves,  the  physical  nature  of  the  Lamb  wave  contribution  to  I /ray  I  becomes  difficult  to 
understand.  The  Lamb  wave  on  a  440c  stainless  steel  flat  plate  of  thickness  h  in  vacuum 
has  a  cutoff  frequency  corresponding  to^^  x  «  76.8,  but  the  dispersion  curve  in  Fig.  A2  for 
,si  does  not  demonstrate  an  abrupt  cutoff.  (The  dispersion  curve  near  a  cutoff  should  be 
particularly  sensitive  to  the  fluid  loading.)  Furthermore,  the  group  velocity  for  the  leaky 
Lamb  wave  becomes  negative  for  jc  <  71  (see  Appendix  A).  As  discussed  below,  the 
presence  of  a  longitudinal  resonance  and  the  behavior  of  the  4’i  leaky  Lamb  wave  for  x  <  7 1 
obfuscates  a  physical  interpretation  of  the  scattering  process. 

Figure  6a  compares  I  ^pws  I  and  I  /ray  I  where  the  curvature  correction /cc  is  omitted 
in  Eq.  (15).  The  ray  synthesis  works  well  forx  <  70,  but  for  x  >  70  the  synthesis  fails  to 
model  the  exact  result.  In  particular,  in  the  vicinity  of  x  =  76.8  the  ray  synnesis  gives  a 
completely  erroneous  representation  of  the  form  function  for  backscatiering.  However,  it 
is  evident  from  Fig.  6a  that  the  positions  of  the  sharp  resonance  features  caused  by  the 
presence  of  leaky  Lamb  waves  seem  to  be  correctly  predicted.  Figure  6b  compaies  I /ray  I 
and  I /pws  I.  hut  fee  is  now  included  in  the  synthesis.  It  i;.  immediately  evident  that  the 
inclusion  of  the  curvature  correction  improves  the  synthesis.  Although  I /ray  I  does  not  fully 
replicate  I /pws  Fig-  hb  indicates  that  the  contribution  from  a  longitudinal  resonance  (as 
presently  modeled)  is  an  important  contribution  in  the  exact  result.  Furthermore,  the  ray 
method  employed  in  deteimining  an  expression  for fsp  necessarily  includes  the  influence  of 
the  shell's  curvature  explicitly  in  fee- 

An  examination  of  the  absolute  error  may  give  some  insight  into  the  cause  of  the 


Fig.  6  The  relevant  leaky  Lamb  v/ave  contributions  are  cq,  ai,  sq,  and  J i.  In  (a)  the 
curvature  correction/cc  is  omitted  in  Eq.  (15).  The  agreement  between  i/ray  I  (solid  line) 
and  l/pws  I  (dashed  line)  is  good  for  z  <  70,  but  I/ray  I  fails  to  model  I /pws  I  forx  >  70. 
Inclusion  of  fee  in  I/ray  I  in  (b)  replicates  l^ws  I  for  jc  <  70  anda:  >  74,  but  the  agreement 
above  jc  =  74  is  marginal.  At  present,  the  cause  of  the  anomaly  in  I  /ray  I  near  a:  =  71  is  not 
fully  understood  (see  Secs.  4.4  and  4,5  and  Appendix  C). 
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discrepancy  between  I  ^pws  I  and  I /ray  I  nearjc  =  71.  A  measure  of  the  absolute  error  in  the 
ray  synthesis  is  given  by  Ej  =  I  ^ws  -/ray  I.  j  =  1, 2.  Figure  7  contains  three  curves.  The 
solid  lint  is  £i  where  the  curvature  correction  in  Eq.  (15)  has  been  omitted.  The  short- 
dashed  line  is  £2  when/cc  is  included  iwfray-  The  long-dashed  line  corresponds  to  \fcc  I 
where  Eq.  (14)  is  employed.  The  reduction  in  Ej  when/cc  is  included  is  given  roughly  by  I 
fee  I.  The  cause  of  the  peak  near  j:  =  76.8  in  £1  is  due  to  the  omission  of  fee  ^^^fray  and  the 
peak  near  jc  =  71  must  not  be  related  exclusively  to  a  longitudinal  resonance. 

The  identification  of  specific  resonance  features  associated  with  the  Ith  leaky  Lamb 
wave  can  be  achieved  by  comparing  1/ 1  in  Figs.  B1  and  B2  with  l^pws  I  In  Fig.  6a.  The 
narrow  resonance  structure  (at  x  =  62.5,  65.0,  67.5,  69.8,  72.0,  74.2,  76.3,  and  78.7)  is 
a  result  of  the  ui  leaky  Lamb  wave  contribution.  The  broad  structure  is  a  combination  of 
the  ao  and  jq  leaky  Lamb  waves.  Since  the  widths  of  the  oq  and  so  resonances  are  fairly 
broad,  it  is  difficult  to  unambiguously  identify  a  feature  in  I  /pws  I  wi^  either  Lamb  wave. 
However,  it  is  apparent  from  the  magnitude  of  each  contribution  that  neither  the  oq  nor  the 
50  leaky  Lamb  wave  is  responsible  for  the  significant  difference  between  l/iay  I  and  I  ypws  I 
near  j:  =  71.  The  final  contribution  to  the  ray  synthesis  is  the  5i  leaky  Lamb  wave.  The 
Fabry-Perot  representation  as  developed  in  Sec.  4.2. A  accounts  for  the  partial 
circumnavigation  and  all  subsequent  circumnavigations  of  a  Lamb  wave  about  the  shell. 
Furthermore,  inspection  of  Eq.  (6)  and  Fig.  A1  demonstrates  that  the  large  radiation 
damping  of  the  si  Lamb  wave  effectively  makes  its  contribution  to  l/^y  I  negligible  for  Jc  < 
72.  It  is  noteworthy  that  our  ray  model^^  of  the  snucture  in  the  forward  scattering 
amplitude  (and  the  total  scattering  cross  via  the  optical  theorem)  for  the  same  shell  has  no 
anomalies  near  either  jt  =  71  or  76.8.  Hence,  the  deficiency  in  the  present  model  does  not 
affect  the  forward  scattering. 


Fig.  7  The  absolute  error  between  the  ray  synthesis  and  the  exact  partial-wave  series  is 
given  by  Ej  =  I  Tpws  '/ray  •>)  =  1.2.  The  solid  line  is  E\  where  fee  in  Eq,  (15)  has  been 
omitted.  The  short-dashed  line,  £2.  includes  the  effects  of  fee  on  /ray.  For  comparison,  the 
long-dashed  line  is  I  fee  I.  The  inclusion  of  the  curvature  correction  in/ray  accounts  for  the 
presence  of  a  longitudinal  resonance  at  xm  ~  76.8.  Furthermore,  the  longitudinal 
resonance  is  not  responsible  for  the  anomaly  near  a:  =  71. 
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Inspection  of  Figs.  1  and  A2  suggests  a  possible  cause  of  the  anomalous  behavior 
near  x  =  7 1.  As  is  decreased  to  xlr  =  76.8  and  below,  cilc  becomes  large  for  the  ^ i 
Lamb  wave.  As  noted  above,  the  frequency  for  which  x  =  xlr  corresponds  to  the  cutoff 
frequency  of  the  mode  for  die  corresponding  flat  plate  in  vacuum  of  the  same  thickness  k 
=  a-b.  Hence,  the  trace  velocity-matching  angle  0/  becomes  small  and  the  wave  may 
radiate  nearly  backwards  directly.  This  corresponds  to  the  radiation  along  ray  B"A"  (with 
B  and  B'  shifted  closer  to  C’  than  illustrated  in  Fig.  1).  Such  radiation  was  not  included  in 
Eq.  (15)  or  elsewhere  in  the  present  analysis.  It  is  plausible  that  such  rays  can  contribute  to 
backscattcring  even  without  being  backward  directed  if  the  rays  intercept  the  Fresnel 
volume^^  of  a  backwards  directed  ray.  The  essential  concept  is  that  rays  in  effect  occupy  a 
region  of  space  (the  Fresnel  volume)  and  that  sound  radiated  primarily  in  some  direction 
(say  along  B"A")  can  contribute  to  the  scattering  in  other  directions  (e.g.,  the  backwards 
axis).  Further  support  for  this  mechanism  is  evident  from  inspection  of  Fig.  A1  which 
shows  that  P/  for  /  =  ^  i  is  very  large  for  x  in  the  region  of  interest.  The  large  p/  indicates 
the  51  leaky  Lamb  wave  is  strongly  coupled  to  the  acoustic  Field  near  the  shell,  thus 
enhancing  the  radiation  in  the  general  direction  of  the  ray  reflected  with  a  local  angle  of 
incidence  «  0/. 

While  the  aforementioned  mechanism  for  the  anomaly  nearx  «  71  has  not  been 
quantitatively  modeled,  additional  support  can  be  seen  by  inspection  of  the  group  velocity 
plotted  in  Fig.  A3  for  the  /  =  5i  mode  of  this  shell.  This  mode  is  seen  to  have  a  negative 
group  velocity  for  j:  <  71  when  computed  by  the  Watson  transfonn  methodology.  Tliis 
suggests  that  energy  can  be  radiated  back  towaids  the  source  without  circumnavigating  the 
shell  It  can  be  argued  by  inspection  of  Eqs.  (39)  and  (40)  of  Ref.  1  that  the  final 
approximation  for  G/,  Eq.  (7),  may  break  down  if  Cgi  is  negative.  The  group  velocity 
anomaly  is  further  discussed  in  Appendix  A. 
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The  influence  of  a  longitudinal  resonance  on  the  form  function  has  also  been  seen 
for  a  set  of  material  parameters  for  an  aluminum  shell.  Figure  8  shows  l/pws  I  and  \fsp  I 
where  the  parameters  are:  ci  =  6.42  hn/s,  =  3.04  km/s,  pe  =  2.70  g/crn^,  c  =  1.4825 
km/s,  and  p  =  1.00  glcrn^.  The  thickness  of  the  shell  is  (hfa)  =  (1  -  b/a)  =  0.04.  The  solid 
line  in  Fig.  8  is  l^ws  I  and  the  dashed  line  is  I  fsp  I.  From  Eq.  (la),  a  longitudinal 
resonance  is  predicted  to  occur  at  xjji  ~  340.  Clearly,  the  dip  near  x  =  340  corresponds  to 
the  presence  of  the  longitudinal  resonance.  Although,  a  detailed  analysis  of  the  surface 
guided  elastic  wave  contributions  for  the  aluminum  shell  has  not  been  carried  out, 
presumably  the  calculated  resonance  structure  is  a  result  of  leaky  Lamb  waves.  The  large 
peak  near  jc  =  3(X)  may  be  a  manifestation  of  prompt  radiation  along  path  B"A"  as  noted 
above. 


4.5  Discussion  and  conclusion 

The  canonical  problem  of  the  backscattering  of  an  acoustic  plane  wave  from  an 
elastic  spherical  shell  has  been  re-examined  by  novel  ray  techniques.  A  form  function  for 
backscattering  including  leaky  Lamb  waves  and  a  specular  reflection  has  been  developed 

qriH  TfiA  pQf»r\r-P#»rr»t  #»vnr#»ccir%n  in  Pn  ffw  a  If^alrv  T  .amH  wavft  hflH 

previously  been  established^' and  the  subsequent  approximation  of  Gi  by  Eq.  (7)  had 
been  verified  for  backscattering  from  spheres  and  shells.^®  The  calculations  presented  here 
give  further  verification  of  the  applicability  of  Eqs.  (6)  and  (7)  in  describing  leaky  Lamb 
wave  contributions  to  backscattering.  The  ray  synthesis  of  the  specular  reflection  is 
composed  of  a  component  associated  with  reflection  from  a  flat  elastic  plate  of  a  normal 
incidence  plane  wave  and  a  novel  curvature  correction  contribution /cc.  A  comparison  of 
the  partial-wave  representation  and  ray  synthesis  of  the  form  function  in  Sec.  4.3 
demonstrated  regions  of  x  where  the  exact  result  and  synthesis  were  in  excellent  agreement. 


It  is  shown  in  Sec.  4.4  that  inclusion  of/cc  in  fray  accounted  for  the  resonance  structure  in 
the  form  function  caused  by  a  longitudinal  resonance.  In  Sec.  4.4,  a  comparison  of  the 
exact  and  synthesized  results  illustrated  that  the  present  ray  model  breaks  down  in  a  range 
of  X  where  the  5i  leaky  Lamb  wave  was  found  to  have  a  negative  group  velocity'.  Finally, 
although  a  ray  model  is  a  high-frequency  approximation,  the  synthesis  displayed  here  for  a 
440c  stainless  steel  shell  was  applicable  dowt.  tox-l. 

While  the  inclusion  of  fee  in  Eq.  (15)  (and  Fig.  6b)  improves  the  synthesis  in  the 
vicinity  of  x  =  xlr,  it  is  noteworthy  that  the  use  of  fee  in  Eq.  (15)  tends  to  increase  the 
discrepancy  between  I  ^ws  I  and  I  fray !  fora:  <  12.  Figure  9  compares  l^ws  !  and  \fray  I 
where  fee  =  0  in  Eq.  (15).  That  is,  the  specular  reflection  contribution  to/is  modeled  as 
the  specular  reflection  from  a  vacuum-backed  flat  elastic  plate  at  normal  incidence. 
Inspection  of  Figs.  4a  and  9  suggests  that  for  x  <  12,  fray  with  fee  =  0  gives  a  more 
accurate  representation  of  the  exact  I ypws  I-  This  niay  be  anticipated  since  the  assumptions 
used  in  the  modeling  of fsp  as)p  +  fee  breakdown  for  sufficiently  low  frequencies. 

Tliree  improvements  to  the  current  ray  syntiiesis  which  should  extend  the  range 
of  applicability  are  as  follows.  First,  extension  of  the  ray  synthesis  to  smaller  values  of  x 
requires  the  inclusion  of  contributions  from  the  subsonic  portion  of  the  oq  leaky  Lamb 
wave  and  Franz-type  waves.  These  were  not  included  in  the  present  model  since  for 
subsonic  waves  Eq.  (4)  predicts  siniQi )  >  1.  (During  the  course  of  this  research,  a  method 
for  including  such  contributions  for  thin  shells  has  been  proposed.^)  Presumably,  these 
types  of  surface  guided  elastic  waves  contribute  to  the  structure  in  I/I  forx  <  7.  Second, 
improvements  in  the  approximation  for  the  specular  contribution /vp  when  x  is  not  large 
may  be  needed.  Third,  a  thorough  investigation  of  the  si  leaky  Lamb  wave  behavior  in  the 
vicinity  of  the  anomaly  discussed  in  Sec.  4.4  should  be  conducted.  As  suggested  in  Sec. 
4.4,  a  modified  ray  picture  may  be  needed  in  this  region.  The  consequetice  of  a 


Fig.  9  The  dashed  line  is  I  fpws  '  ^nd  the  solid  line  is  I/ray  i  for  the  stainless  steel  shell. 
The  leaky  Lamb  waves  included  in/ray  oq  and  sq.  The  present  figure  differs  from  Fig. 
4a  by  setting/cc  =  0  in  Eq.  (15)  so  that  the  curvature  correction  to  the  reflected  wave 
amplitude  is  neglected.  The  improved  agreement  over  part  of  this  x  region  between  I /pWS 
and  I  /ray '  is  discussed  in  Sec.  4.5. 
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negative  group  velocity  on  the  ^  i  contribution  needs  to  be  more  fully  explored  as  well  as 
any  connection  bem'een  the  first  longitudinal  resonance  and  the  cutoff  of  the  si  Lamb  wave 
for  curved  fluid  loaded  shells. 

Comments  on  the  significance  of  the  fractional  thickness  hja  are  appropriate.  For 
the  present  numerical  investigation  for  the  stainless  steel  shell,  we  took  hja  =  0.162  to 
correspond  to  the  shell  studied  in  the  scattering  experiments  in  Ref,  12  and  subsequent 
theoretical  studies.  15^7,23  Some  of  the  changes  which  may  be  anticipated  for  smaller  hla 
will  now  be  noted.  For  flat  plates  in  a  vacuum,  the  cutoff  frequencies  increase  with 
decreasing  h.  While  the  dispersion  relations  near  cutoffs  may  be  strongly  affected  by  fluid 
loading,  as  a  general  rule  it  may  be  anticipated  that  for  a  given  frequency,  fewer  leaky 
Lamb  modes  will  be  required  as  h  is  decreased.  The  frequency  for  which  ci  of  the  lowest 
flexural  (or  antisymmetric)  Lamb  mode  equals  c  of  tire  surrounding  fluid  is  commonly 
referred  to  as  the  "coincidence  frequency."  (See  also  comments  in  Appendix  A  on  the 
bifurcation  of  the  oq  dispersion  curve.)  While  for  the  shell  considered,  coincidence  occurs 
for  X  «  7;  smaller  values  of  hla  would  raise  the  coincidence  frequency.  Even  away  from 
cutoff  and  coincidence  frequencies,  as  hJa  is  decreased,  the  effect  of  the  fluid  loading  on  the 
phase  velocity  curves,  Fig.  A2,  will  be  more  pronounced.  It  may  be  anticipated  that  ciJc 
(as  given  by  the  Watson  methodology)  will  be  shifted  below  the  plate  in-a-vacuum 
curvature-corrected  values.  Finally,  the  longitudinal  resonance  condition,  Eq.  (la),  may  be 
expressed  as  xir  =  n7c(c//c)[/i/a]'L  Consequently,  for  a  given  n,  xjji  increases  with 
decreasing  hja. 

While  the  emphasis  of  the  present  research  has  been  on  scattering  in  the  exact 
backwards  direction,  the  generalization  of  leaky  Lamb  wave  contributions  to  near 
backwards  (but  off-axis)  directions  follows  from  the  discussions  of  axial  focusing  given  in 
Refs.  1,12,  and  19.  Insight  into  the  high-frequency  near-backwards  scattering  patterns  for 
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shells  which  are  only  slightly  spheroidal  may  be  obtained  by  adapting  the  analysis  given  in 
Ref.  22. 
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Appendix  A,  Leaky  Lamb  wave  parameters  p/  and  c//c 

A  Sommerfeld-Watson  transformation  has  not  been  directly  applied  to  the  exact 
partial-wave  series  for  the  backscattering  from  a  evacuated  elastic  spherical  shell. 

However,  some  results  from  a  Sommerfeld-Watson  transformation  of  the  fonn  function  for 
backscattering  from  an  solid  elastic  sphere  can  be  extended  to  the  case  of  a  shell.  It  can  be 
argued  from  the  geometric  similarity  of  the  two  scatterers  that  certain  results  of  a 
SommerfelC.-Watson  transformation  of / for  the  shell  will  produce  analogous  results.  In 
Sec.  4.2. A.  a  ffeometric  intemretatinn  of  the  backscatterinp  of  surface  eiiided  elastic  waves 
from  a  solid  elastic  sphere  was  employed  in  modeling  the  backscattering  from  a  shell.  In 
this  Appendix,  the  methodology  for  the  determination  of  p/(x)and  ci{x)lc  is  briefly 
reviewed. 

From  Williams  and  Marston's  analysis,7  the  radiation  damping  parameter  and 
surface  guided  elastic  wave  phase  velocity  are  obtained  by  solving  the  following  equations: 


Dv/(x)  =  0, 

(Al) 

V/  =  a/  +  ip/. 

(A2) 

(c//c)  =  x/(a/  +  1/2). 

(A3) 
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Equation  (Al)  is  the  deteiminant  of  the  denominator  in  Eq.  (3)  where  the  integer  index  n 
has  been  replaced  everywhere  by  a  complex  index  v/.  The  complex  index  V/  for  fixed  Jc  is 
obtained  by  numerical  techniques  described  in  Ref.  12.  These  numerical  techniques 
involve  the  evaluation  of  residues  of  certain  integrals.  Once  v/  is  determined,  then  both  (3/ 
and  (c//c)  via  Eq.  (A3)  are  known.  The  radiation  damping  parameters  P/(jc)  for  the  various 
leaky  Lamb  waves  are  shown  in  Fig.  Al.  An  important  feature  in  Fig.  Al  is  that  (3/  for  the 
ao.  ■SQ. ,  and  si  leaky  Lamb  waves  are  relatively  small  throughout  the  range  of  x 
investigated;  while  P/  for  the  leaky  Lamb  wave  appears  to  diverge.  The  truncation  of  the 
vertical  axis  at  0.8  Np/rad  resolves  the  detail  of  some  the  Lamb  wave  damping  parameters, 
but  it  obscures  the  fact  that  the  damping  parameter  becomes  large  forx  <  72  (P;  ~  3.5  at 
j;  =  70  and  P/ «  7.5  at  x  =  65).  The  derivation!  leading  to  Eq,  (7)  made  the  assum.ption 
27cP/  «  1.  It  may  be  argued,  however,  that  Eq.  (7)  is  a  good  approximation  for  G/  even 
when  2Ttp/  violates  this  assumption,!  Inspection  of  Figs.  4-6  supports  the  claim  that  Eq. 
(7)  with  <j)/  =  0  is  a  valid  approximation  for  G/  even  if  27ip/  is  not «  1, 

The  normalized  dispersion  curves  for  the  leaky  Lamb  waves  are  shown  in  Fig. 
A2.  Mode  identification  was  confirmed  by  comparison  with  curvature  corrected  results  for 
a  flat  plate  in  vacuum  following  the  method  of  Ref.  17.  Several  characteristics  should  be 
noted.  First,  the  ao  leaky  Lamb  wave  (solid  line)  becomes  subsonic  for  x  <  7. 
Sammelmanti  et  al.^®  have  reported  a  bifurcation  of  the  ao  leaky  Lamb  wave  dispersion 
cuive  for  the  fluid-loaded  spherical  shell  near  a  value  of  x  associated  with  the  transition 
from  a  supersonic  to  subsonic  phase  velocity  for  the  dispersion  curve  of  the  ao  Lamb  wave 
on  a  spherical  shell  in  vacuum.  Since  the  relevant  contributions  to  the  present  ray  synthesis 
are  implicitly  supersonic,  the  subsonic  branch  was  not  included  in  Fig.  A2.  The  dispersion 
curves  for  ai  (long-dashed  line)  and  S2  (dot-dashed  line)  are  analogous  to  the  flat-plate 
Lamb  wave  dispersion  curves.  Both  the  aj  and  ^^2  modes  appear  to  approach  well  defined 
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Fig.  A 1  The  radiation  damping  parameters  for  the  various  leaky  Lamb  waves  employed 
in  I/ray  i  for  the  stainless  steel  shell  considered.  Each  leaky  Lamb  wave  is  associated  with 
the  following  line:  ao  solid;  sq  short-dashed; 
lashed. 


ai  long-dashed;  i  j  short-long-dashed;  and  i'2 
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Fig.  A2  The  normalized  dispersion  curves  for  the  various  leaky  Lamb  waves  employed 
in  I/ray  I-  Eacli  Lamb  waves  is  associated  with  the  following  line:  oq  solid;  jq  short- 
dashed;  (i\  long-dashed;  short-long-dashed;  and  S2  dot-dashed.  The  dotted  lines  are 
based  on  corresponding  mode  calculations  for  a  flat  plate  in  a  vacuum  with  curvature 
conections  as  given  by  Marston  (Ref.  17).  These  are  included  so  as  to  confirm  the 
classification  of  modes  given. 
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cutoff  frequencies.  As  in  Fig.  Al,  the  vertical  axis  has  been  truncated  and  the  asymptotic 
behavior  of  these  modes  i.s  not  adequately  represented.  For  example:  for  /  =  ai,  c//c  =  27 
at  X  =  42,  and  for  /  =  S2,  ci/c  -  32  at  x  =  85. 

The  phase  velocity  ratio  for  the  sq  Lamb  wave  is  represented  by  the  short-dashed 
line  in  Fig.  A2.  Two  features  are  observed  from  this  dispersion  curve.  First,  the  5o  Lamb 
wave  is  weakly  dispersive  in  the  ranges  10  <x  <  40  andx  >  80.  As  discussed  in 
Appendix  B,  the  spacing  Ax^^^  between  adjacent  resonance  peaks  in  Fig.  B2  is 
approximately  equal  to  the  normali  zed  group  velocity  Cgi/c  for  weakly  dispersive  Lamb 
waves.  Second,  tlie  dispersion  curve  appears  to  diverge  as  the  frequency  approaches  the 
natural  frequency  0)b  of  a  purely  radial  oscillation  or  breathing  mode  of  the  shell.  Such  a 
divergence  is  plausible  since  the  radial  displacement  of  the  shell  at  all  points  on  the  surface 
is  in  phase  when  (O  =  0)£.  In  Appendix  A  of  Ref.  15,  cog  for  a  thin  spherical  shell  in 
vacuum  is  related  to  the  ring  frequency  O)/?  of  an  infinite  cylindrical  shell,  ^1  The  estimated 
value  of  X  corresponding  to  cog  is  xg  =  6.3.  Funhermore,  it  is  argued  in  Ref.  15  that  the 
fluid-loading  of  the  spherical  shell  will  only  slightly  reduce  xg. 

The  normalized  phase  velocity  for  the  Lannb  wave  is  the  final  dispersion  curve  in 
Fig.  A2  to  be  considered  (the  long-short  dashed  line).  Unlike  the  a\  and  $2  leaky  Lamb 
waves,  the  dispersion  curve  for  the  si  leaky  Lamb  wave  on  the  fluid-loaded  shell  does  not 
approach  a  definite  cutoff  frequency.  The  cutoff  frequency  for  the  equivalent  si  Lamb 
wave  on  a  flat  elastic  plate  in  a  vacuum  corresponds  to^^  x  =»  76.8.  Near  x  =  71,  the  si 
leaky  Lamb  wave  attains  a  maximum  value  and  then  decreases  with  decreasing  x.  The 
group  velocity  is  related  to  the  phase  velocity  by^2 


cgl  =  c/{  1  -  [1  -  (c.7/x)(dc//dx)-lj-I). 


(A4) 
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Numerical  differentiation  indicates  that  on  a  fluid-loaded  shell  has  a  negative  group 
velocity  in  the  region  65  <  a  <  71.  The  normalized  group  and  phase  velocities  are 
displayed  in  Fig.  A3.  The  possibility  of  a  wave  propagating  with  a  negative  group  velocity 
was  first  investigated  by  Lamb.^^  Hackman  and  Sammelmann33  have  reported,  that  an 
analysis  of  die  poles  of  the  S-matrix  in  the  complex  ^'-plane  produces  regions  of  negative 
group  velocities  for  some  leaky  Lamb  waves  on  thin  elastic  spherical  shells.  Finally,  the 
dispersion  curve  has  been  truncated  below  a  =  65  because  the  numerical  algorithm 
implemented  in  the  determination  of  c//c  becomes  inaccura  te. 

Without  resorting  to  the  above  numerical  differentiation,  the  sign  of  the  group 
velocity  for  a  leaky  Lamb  wave  on  a  fluid-loaded  shell  can  be  determined  from  the  slope  of 
a/(jc).  If  Eq.  (A3)  is  differentiated  with  respect  to  x,  then  the  result  is 


(  j  \ 

da,  ^  V  dc, 

dx  -  C,  C,  dx  J. 


(A5) 


Equation  (A4)  relating  Cgi  and  c/  can  be  inverted  to  obtain  an  expression  for  the  second  term 
in  Eq.  (A5).  Substituting  this  result  into  Eq.  (A5)  gives  the  simple  result  a}(a:)  =  cjcgi. 


Although  a  plot  of  (xi(x)  is  not  pre.sented  here,  it  has  been  verified  that  a/(x)  has  a  positive 


slope  for  x>  71,  attains  a  local  minimum  near  a:  =  71,  and  has  a  negative  slope  for  x  <  71. 


Appendix  B.  Fabry-Perot  expression  for  // 

In  this  appendix,  the  ao  and  ;7o  leaky  Lamb  wave  contributions  to  the 
backscattering  amplitude  are  isolated  and  examined.  By  considering  an  individual 
contribution,  the  resonant  nature  of  Eq.  (6)  is  demonstrated  and  the  similarity  between  Eq. 
(6)  and  a  Fabry-Perot  resonator  is  illustrated.  Figure  B 1  is  I// 1  where  /  =  no  and  Fig.  B2 
corresponds  to  I// 1  for  the  leaky  Lamb  wave.  The  width  of  any  given  resonance  pesik 
from  the  flo  Lamb  wave  contribution  in  Fig.  B1  is  representative  of  relatively  low  Q 


C//C,  Cq//C 


F'  3.  A3  The  nornialized 


group  velocity  (solid)  for  the  5]  leaky  Lamb 
and  longitudinal  resonance  at^^R  =  76.8. 


73)  are  a  result  of  numerical 
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resonances  while  Fig.  B2  for  sq  contains  both  high  Q  resonances  (x  <  40)  and  low  Q 
resonances  (x  >  40).  Ordinarily,  the  increment  in  x  for  these  calculations  of  I// 1  (I  =  ao,  -so) 
is  Ax  =  0.015875.  In  Fig.  B2,  the  high  Q  resonances  of  the  sq  Lamb  wave  a.tx~  17.6, 
21.3,  25.0  and  28.6  have  been  enhanced  by  sufficiently  decreasing  Ax  to  resolve  the  very 
narrow  Jinewidths.  It  may  be  inferred  from  Fig.  B2  and  Eq.  (16)  that  a  partial-wave  series 
calculation  will  be  numerically  intensive,  since  Ax  must  be  chosen  small  and  rimax  exceeds 


x. 

Consider  the  form  of  the  denominator  in  Eq.  (6).  If  a  Lamb  wave  on  the  shell  is 
non-dispersive,  then  the  spacing  between  two  adjacent  resonances  would  be  Axres  =  ci/c. 
For  weakly  dispersive  Lamb  waves,  it  is  well-known  that  Axres  “  Cgifc.  The  relationship 
between  Axres  and  Cgi  can  also  be  obtained  from  the  '^orm  of  the  denominator  in  Eq.  (6).^4 
Figure  A2  demonstrates  that  the  oq  leaky  Lamb  wave  is  weakly  dispersive  forx  >  50  and 
Cgilc  «  2. 178  (see  Fig.  6  of  Ref.  12).  Inspection  of  Fig.  B 1  indicates  that  the  spacing  AXres 
forx  >  50  is  Axres  ”  2.123  (for  comparison  1.947  <  c//c  <  2.037,  50  < x  <  100).  Again, 
the  dispersion  curve  in  Fig.  A2  indicates  the  5o  leaky  Lamb  wave  is  weakly  dispersive  in 
the  regions  10  <  x  <  40  and  x  >  80.  The  spacing  between  adjacent  peaks  in  the  region  10  < 

V  ^  4  O  ^  M  Lrfk  /T  13  O  ^4^  W  A  A  V  ^  1 A  1  t'1  Ol~l 

^  "TV/  xo  VXViVt  iXUiCVVX  IIVXXI  X  1^.  (.V/  l,;v  vvxtixv  ilUiJlVXiV'U.X  \XXX.XV'XVAXMUU.V/«1 

predicts  Ax^es  “  Cgjc  =  3.616  ( at  x  =  14,  cilc  =  3.978).  A  thorough  investigation  of  the 
dependence  of  AXres  on  c//c  and  Cgilc  for  dispersive  leaky  Lamb  waves  is  beyond  the  scope 
of  the  present  appendix. 

Also,  contained  in  Figs.  B 1  and  B2  are  envelops  of  the  maximum  and  minimum 
values  that  I  // 1  can  attain.  From  Eq.  (6),  I //  Imax  and  I  fi  Imin  occur  when  expiilKxc/ci)  =-l 
and  expiUiixc/ci)  =  1,  respectively.  Inspection  of  Eqs.  (6)  and  (7)  demonstrates  that  I//  \q 
(q  =  min,  max)  is  an  implicit  function  of  x  through  Ute  dependence  of  (i;  and  ci/c  on  x.  The 
importance  of  the  lower  envelop  is  that  it  describes  a  smooth  background  associated 


130 


ka 


Fig.  B1  The  Fabry-Perot  representation  of  the  oo  leaky  Lamb  wave  contribution  to  the 
form  function  for  backscattering  from  an  evacuated  elastic  spherical  shell.  The  long-dashed 
line  is  the  maximum  value  that  I// 1  can  attain  and  the  short-dashed  line  is  a  minimum 
envelop  which  can  be  interpreted  as  a  smooth  background  associated  with  off-resonance 
contributions  from  the  oq  Lamb  wave.  The  widths  of  the  resonances  are  representative  of 
relatively  low  Q  leaky  Lamb  wave  resonances. 
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Fig.  B2  The  Fabry-Perot  representation  of  the  sq  Lamb  wave  contribution  to  the  form 
function  for  backscattering.  The  long-dashed  line  is  the  maximum  value  that  I // 1  can  attain. 
The  short-dashed  line  is  a  minimum  envelop  which  can  be  interpreted  as  a  smooth 
background  associated  with  off-resonance  conoibutions  from  jq-  Th^  widths  of  the 
resonances  for  a:  <  40  are  representative  of  high  Q  leaky  Lamb  wave  resonances. 
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with  the  Ith  leaky  Lamb  wave.  The  resonance  enhancement  rises  above  the  background. 
Unlike  resonance  scattering  theory35-37  v/here  a  suitable  background  contribution  must  be 
chosen  to  synthesis  the  backscattering  amplitude,  the  Fabry-Perot  form  of  Eq.  (6)  correctly 
accounts  for  a  background  contribution  even  for  off-resonance  values  of  x.  That  is,  the 
resonance  peaks  in  I// 1  rise  relative  to  i //  Imin  which  represents  a  off-resonance  smooth 
background.  ^ 


Appendix  C.  Localization  principle  analysis  of  partial  waves  near  a 
longitudinal  resonance 

The  localization  principle  states  that  the  nth  partial-wave  amplitude  for  scattering 
from  a  sphere  may  be  associated  with  the  contribution  to  the  scattering  from  a  ray  having  an 
impact  parameteri"^ 

^  =  a(n  +  l/l)/x.  (Cl) 

For  s>  a,  rays  miss  the  sphere,  and  the  partial  wave  amplitudes  typically  decrease  abruptly 
in  magnitude  as  n  increases  above  x.  Here,  we  analyze  the  elastic  response  of  the  shell  by 
considering  the  contribution  of  individual  partial-wave  amplitudes  in  the  form  function  for 
backscattering.  The  elastic  behavior  of  the  sheU  is  isolated  by  subtracting  a  rigid 
background.  The  relevant  expression  in  the  present  analysis  is 


r 


F(x,n)  =  f„(x)  -  f;,(x)  =  (-l)"(2n+l) 

IX 


Bn(x)  ^  jn(x) 


(1)' 

hn  (X) 


(C2) 


where  the  prime  indicates  differentiation  with  respect  to  the  argument.  The  second  term  in 
Eq.  (C2)  represents  the  subtraction  of  an  acoustically  hard  or  rigid  immovable  sphere 
background. In  particular,  the  cause  of  the  anomalous  peak  in  the  ray  synthesis  near  x  = 
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71  and  a  possible  relationship  between  the  longitudinal  resonance  and  the  ,>  i  leaky  Lamb 
wave  are  considered  below. 

Figure  Cl  shows  I  F(x,n)  I  for  jc  =  71.7  which  corresponds  to  a  value  of  x  at  the 
center  of  the  peak  in  E2  of  Fig.  7.  Figure  Cl  illustrates  the  localization  principle  since  the 
important  contributions  are  the  partial-wave  amplitudes  at  n  =  8  - 11, 15,  27,  and  35.  These 
correspond  to  various  leaky  Lamb  wave  contributions  as  well  now  be  shown.  Inspection 
of  Fig.  1  and  Eq.  (4)  demonstrates  that  the  impact  parameter  is  related  to  the  outer  radius  of 
the  shell  by 

s  =  a  sm(0/)  =  acici.  It  follows  that  a  partial  wave  amplitude  can  be  associated  with  a  leaky 
wave  contribution  by 

n  =  x(c/cp-(l/2).  (C3) 

Table  Cl  gives  the  result  of  applying  Eq.  (C3).  Inspection  of  Fig.  Cl  and  Table  Cl  indicate 
that  the  ao,  sq,  and  a\  leaky  Lamb  wave  contributions  are  isolated  and  associated  with  a 
single  partial  wave  amplitude.  An  interpretation  of  the  n  =  8  - 11  partial- wave  amplitudes  is 
that  a  broad  bundle  of  rays  couples  to  the  leaky  Lamb  wave  on  the  shell  and  and  that 
wave  contributes  to  backscattering.  This  may  be  consistent  with  the  direct  backscattering 
mechanism  proposed  in  Sec.  4.4. 

One  may  anticipate  some  relationship  between  the  longitudinal  resonance  and 
the  Lamb  wave,  since  the  resonance  condition  in  Eq.  (la)  gives  xm  «  76.8  which  is 
equivalent  to  the  cutoff  frequency  of  .vi  on  a  vacuum-backed  flat  elastic  plate.  Information 
concerning  the  longitudinal  resonance  and  the  S]  leaky  Lamb  wave  can  be  obtained  by 
studying  the  contributions  to  the  form  function  due  to  individual  partial  waves.  Figure  C2 
displays  a  sequence  of  partial-wave  amplitudes  as  a  function  of  x  in  the  vicinity  of  xlm- 
The  series  of  large  peaks  starting  near  jc  =  55  and  n  =  16  corresponds  to  the  59  mode  and 
the  set  of  narrow  peaks  is  due  to  the  ai  leaky  Lamb  wave.  Since  an  estimate  of  the 
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Fig.  Cl  Localization  principle  analysis  of  the  elastic  response  of  the  shell  dxx~l  l.1. 
Each  peak  is  associated  with  a  particular  leaky  Lamb  wave  contribution.  A  particular  leaky 
Lamb  wave  can  be  associated  with  a  given  peak  by  applying  Eq.  (C3).  (See  Table  CI\ 
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TABLE  CL  Localization  piiiiciple  identification  of  the  partial  wave  amplitudes  2l\.x  =  111 


plotted  in  Fig.  Cl.  The  normalized  phase  velocity  is  given  by  Eq.  (A3)  for  the  stainless 

steel  shell  and  n  is  given  by  Eq.  (C3). 

/ 

Clic 

n 

oo 

2.023 

34.9 

SQ 

2.620 

26.9 

4.678 

14.8 

6.484 

10.5 

I 
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normalized  group  velocity  is  given  by  Cgilc  ~  Axres/^,  then  the  sq  and  ai  peaks 
demonstrate  the  usual  behavior  of  positive  Cg}/c  for  a  surface  guided  elastic  wave.  The 
remaining  set  of  peaks  is  primarily  associated  with  the  5i  leaky  Lamb  wave.  Two 
important  feature  are  immediately  evident.  First,  the  group  velocity  of  the  5]  Lamb  wave 
can  have  either  positive  or  negative  values  as  demonstrated  above  in  Appendix  A.  This  can 
be  seen  by  considering  the  slope,  &Xresl^,  of  a  smooth  curve  through  the  peaks. 

However,  a  discrepancy  between  the  dispersion  curves  from  the  analysis  in  Appendix  A 
and  the  present  method  for  estimating  cg^c  is  also  evident.  The  result  of  Eqs.  (Al)  -  (A3) 
in  Fig.  A2  produces  a  single- valued  function  of  ci{x)lc.  Equation  (C3)  and  the  partial-wave 
amplitudes  in  Fig.  C2  would  suggest,  however,  that  ciix)lc  is  a  multi-valued  function  of  x 
for  /  =  ^1,  Furthermore,  the  dispersion  curve  for  in  Fig.  A2  is  defined  for  <  71,  but 
Fig.  C2  displays  that  no  corresponding  peaks  in  the  partial  wave  amplitudes  to  account  for 
this  region.  A  second  observation  concerning  these  peaks  is  that  as  n  tends  to  zero,  the 
locus  of  X  corresponding  to  these  peaks  appears  to  approach  xm.  The  weak  peak  in  I 
F(x,0)  I  near  x  =  xm  appears  to  be  a  consequence  of  the  use  of  a  rigid  background  in  Eq. 
(C2)  which  does  not  display  a  longitudinal  resonance.  A  calculation  of  I  F{x,n)  I  (not 
shown)  where  the  rigid  background  has  been  replaced  by  an  evacuated  fluid  shell 
background  demonstrates  that  the  peak  in  I  F(jc,0)  I  nearx  =xlr  is  suppressed.  This 
suggests  that  the  peak  in  I  F(jc,0)  i  displayed  in  Fig.  C2  is  associated  with  the  longitudinal 


resonance. 
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Appendix  I 

A.  Complex  coupling  coefficient:  Approximation  and  an  apparent  exact 
expression 

In  the  development  of  the  ray  methods  presented  in  the  preceding  chapters,  a 
complex  coupling  coefficient  is  defined  and  denoted  as  G/.  This  coefficient  describes  the 
efficiency  of  the  coupling  of  the  incident  acoustic  wavefield  in  water  with  a  surface  guided 
elasdc  wave  excited  on  the  elastic  spherical  object  during  the  scattering  process.  For  the 
case  of  backscattering  from  a  solid  elastic  sphere,  a  Sommerfeld-Watson  transformation  of 
the  exact  partial  wave  series  representation  of  the  fomi  function/ gives  a  virtually  exact 
expressionl*2 


G^(x)  =  87c(c/cp 


1  + 


CX  y, 


iC/P/ 


/ 


2%  (x) 


(1) 


The  radiation  damping  of  the  surface  guided  elastic  v/ave  is  P/  (in  Np/rad)  and  x~kaisa 
dimensionless  size  parameter  where  k  is  the  wavenumber  in  the  surrounding  water  and  a  is 
the  radius  of  the  sphere.  The  normalized  phase  velocity  of  the  surface  guided  elastic  wave 
is  cifc.  The  functions  tP'*'  and  •D~  are  extremely  complicated  expressions  involving 


cGiiibinaticns  of  spheiiCoi  Sessel  and  Neumann  functions  of  complex  Oider  v/  and  leal 


argument.  The  functions  and  (D"  for  a  solid  elastic  sphere  are  given  in  Refs.  1  and  2. 
The  dot  above  indicates  differentiation  with  respect  to  v  and  this  derivative  is  to  be 
evaluated  at  v/  where  v/  satisfies  the  relationships 

iPy^(x)=0,  v^  =  a^  +  ip^,  c^(x)/c  =  x/[a^  + 1/2].  (2a,b,c) 

Equation  (2)  is  a  result  of  the  above  mentioned  Sommerfeld-Watson  transformation.  [The 
V/  that  satisfy  Eq.  (2)  also  satisfy  Dv(x)  =  0  where  Dy  is  the  denominator  of  the  nth  partial- 
wave  amplitude fnix)  -  {-l)''(2n+\)(2/ix)[Bn(x)/Dnix)],  with  the  integer  index  n 
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replaced  by  v.]  The  complicated  forni  of  G/  in  Eq.  (1)  obfuscates  the  dependence  of  G/  on 
the  relevant  surface  guided  elastic  wave  parameters  P/  and  c//c.  A  simple  approximation  for 
Eq.  (1)  would  facilitate  computations  of  forward  or  backwards  scattering  by  the 
quantitative  ray  methods  previously  developed. 

Marston  developed  a  simple  approximation  for  I  G/ 1  by  comparing  the  results  of 
resonance  scattering  theory^.^  and  a  generalization  (based  on  Watson  transform  methods) 
of  the  geometrical  theory  of  diffraction.^  A  short  summary  of  this  analysis  follows  and  it  is 
restricted  to  the  backscattering  direction.  In  resonance  scattering  theory  at  the  location  of  an 
isolated  resonance, /« is  split  into  an  appropriate  background  term  and  a  surface  elastic 
wave  contribution /«/.  The  contribution /„/  is  usually  exptessed  in  a  Breit-Wigner  form 


_  (2n+l)  ,  n.i2^ _ td - 

- T” ^  x-(i/2)r^^ 


(3) 


where  Xni  is  the  location  of  an  isolated  resonance  and  the  width  of  the  resonance  can  be 
related  to  Define  the  scattering  function  Sn{x)  -  {[2B„(jc)/Dn(x)]  -  1 }.  The 
denominator  of  Sn(x)  can  be  written  so  as  to  contain  the  information  about  the  resonance 
location  and  damping;  Xni  is  the  real  root  of  the  real  part  of  that  denominator.  Equation  (3) 


is  a  icsult  of  Epjjlyui^  a  Taylor  Series  expansion  to  the  denominator  of  Sfiix)  about  x  —  Xni 


and  separating /„/  from  a  background  contribution.  For  high-frequency  backscattering 
from  a  solid  elastic  sphere, /n  for  a  rigid  or  immovable  sphere  is  a  suitable  choice  of 
background.  It  can  be  shown  that  i?e{Drt(JC  =  Xn/)}  =  0  where  D„  is  the  denominator  of /«. 
The  subscripts  n  and  /  imply  that  the  isolated  resonance  occurs  in  fn  and  it  is  classified  as 
the  Ith  type  of  surface  elastic  wave  resonance  (Rayleigh,  Lamb,  Franz,  etc.).  From  the  ray 
synthesis,  a  surface  elastic  wave  contribution  to  backscattering  has  the  following  Fabry- 
Perot  form”^ 
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ill,  -2(7t-0,)P,  -2nP,  i2jca, 

f/x)  =  [-G^e  ‘'e  '  ^  ‘  ], 


(4) 


where  Tj/  is  a  propagation  phase  delay.  The  basic  premise^  leading  to  an  approximation  for 
Gi  is  that  if  both  Eq.  (3)  and  Eq.  (4)  describe  the  resonance  contribution  of  a  surface  elastic 
wave,  then  locally  near  the  isolated  resonance//(jc)  ~fni{x).  The  remainder  of  the  analysis 
involves  applying  Taylor  series  expansions  to  |3/(x)  and  a/(x)  about  x  =X;i/.  Here,  Xni  is 
a  real  value  of  x  such  that  ot/(X„/)  =  n  and  the  denominator  in  Eq.  (4)  attains  a  local 
minimum.  With  these  expansions,  Eq.  (4)  can  bt;  cast  into  a  form  similar  to  Eq.  (3) 


f/x) » 


-G.e 


iti, 


i47tp 


nl 


r„,  +  i2S 


nl 


where 


®n/  =  -  Pn  A  “n?  *  I'm  =  “m  1' 


(5) 


The  prime  in  Eq.  (6)  indicates  differentiation  of  a/(x)  and  P/(x)  with  respect  to  x  evaluated 
Sit  x=Xni  and  Prt/  =  P/(X„/), 

Inspection  of  Eqs.  (3)  and  (5)  demonstrates  that  the  location  of  poles  in  the  complex 
.r-plane  for  fni  and /;  coincide  if  Xni  =  Xni  +  5^/  and  Tni  ~  "fnl-  Numerical  computations  for 
tungsten  carbide  and  aluminum  spheres  indicate  I  6^/ !  «  1,  which  suggests  the  following 
conditions  I  bnl  i  «  Jnl  and  (dp;/cU) « (dai/dx).  With  these  conditions  and  pole 
locations,  comparison  of  Eqs.  (3)  and  (5)  reduces  1  G/ 1  at  an  isolated  resonance  to  the 
apprordmate  form 


(n+l/2)87ip 


IG/x=X„,)U  ^ 


nl 


nl 


(7) 


The  form  of  Eq.  (2c)  and  the  identification  of  [(n  +  l/2)PCni\  =  [i^li^nl)  +  1/2)/^/;/] 


suggests  the  following  smooth  continuation  of  Eq.  (7) 
I  G^(x)  I  =  87tP^(x)c/c^(x). 


(8) 
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Equation  (8)  is  a  simple  approximation  for  1  Gi  I  which  explicitiy  exhibits  the  dependence  of 
I  G/ 1  on  the  relevant  surface  guided  elastic  wave  parameters. 

While  this  analysis  gives  a  simple  approximation  for  I  G/ 1,  it  was  initially  deemed 
too  complicated  to  construct  an  approximation  for  the  phase  of  G/  from  the  above  analysis. 
However,  Williams  and  Marston  have  subsequently  demonstrated  that  arg(Gi) »  0  for 
surface  guided  elastic  wave  contributions  to  backscattering  from  solid  or  hollow  elastic 
spheres.8'9  Recently,  Ho  and  Felsen  have  confirmed  Eq.  (8)  with  arg(Gi)  ~  0  for 
backscattering  from  thin  elastic  spherical  shells.  It  is  noteworthy  that  their  analysis  is 
based  on  thin  shell  equations  while  the  above  results  are  from  an  analysis  starting  from  the 
exact  partial  wave  series  for  backscattering  from  elastic  spheres. 

A  Watson  transfonnation  of  the  exact  partial  wave  series  representation  of  the  form 
function  for  backscattering  from  an  elastic  spherical  shell  has  not  been  performed.  From 
the  geometric  similarity  of  the  sphere  and  shell,  it  can  be  argued  that  the  complex  coupling 
coefficient  for  the  shell  will  have  the  same  form  as  Eq.  (1).  The  shell  and  sphere  coupling 
coefficients  should  only  differ  in  the  definitions  of  ©■*■  and  CD"  in  Eq.  (1).  The  appropriate 
expressions  are 


!r{(x)  =  4^\x)  -  Fy(x), 

2\,(x)  =  4^\x)  -  Fy(x), 

(9a) 

(1). ,  Jl)'.  .-(1),  . 

(x)  =  xhy  (x)/h^  (x). 

(2),  ..  \ 
z^j  (x)  =  xhv  (x)/hy  (x). 

(9b) 

9  (1)  (2) 

Fy(x)=-xJpDy(x)/PeD;,' 

(X). 

(9c) 

In  Eq.  (9b),  the  hv  are  spherical  Hankel  functions  of  the  first  and  second  kind  of  complex 
order  and  real  argument  and  the  prime  denotes  differentiation  with  respect  to  the  argument. 
In  Eq.  (9c),  Xi  -xcict  whcie  C;  is  the  trrnsvcrse  velocity  of  the  shell's  material  and 
and  are  the  cofactor  expansions  wdh  icspfx^t  to  the  elements  d\\  and  di]  of  D,i  froro/„ 
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with  n  replaced  by  v.  Foi’  solid  spheres,  and  are  2x2  determinants  while  for 
evacuated  shells  and  would  be  4  x  4  determinants  (see  Refs.l  and  3  for  the 
elements  of  the  determinants).  This  suggests  that  the  approximation  in  Eq.  (8)  can  be 
applied  to  shells  for  certain  classes  of  surface  guided  elastic  waves. 

The  remainder  of  this  section  compares  the  consequences  of  using  Eq.  (8)  with 
<^rg{Gi)  =  0  instead  of  Eq.  (1)  for  an  elastic  spherical  shell.  The  material  of  the  shell  is 
440c  stainless  steel  with  the  following  parameters:  longitudinal  sound  speed  cl  = 
5.854/tmA,  transverse  sound  speed  q  =  2>.\50kTnlSy  density  pe  =  1-^^glcrrfi  and  inner-to- 
outer  radii  ratio  bla  =  0.838.  The  parameters  for  water  are  c  =  \A19kmls  and  p  = 
l.Q^glcnC^.  For  these  parameters  and  the  frequency  range  0  <  a:  <  100,  the  relevant  surface 
guided  elastic  wave  contributions  are  the  oq  and  a\  antisymmetric  (or  flexural)  leaky  Lamb 
waves  and  the  5o>  •^1*  and  ^2  symmetric  (or  dilatational)  leaky  Lamb  waves.  [The  present 
ray  synthesis  is  implicitly  supersonic  through  the  use  of  the  phase  velocity  trace-matching 
condition  5m(0/)  =  c/c/.  Since  the  oq  leaky  Lamb  wave  becomes  subsonic  for  x  <  7,  the 
actual  range  investigated  is  7  <x:  <  100.  Also,  for  jc  <  10  contributions  from  Franz-type 
creeping  waves  rr  contribute  significantly  to  backscattering.]  Finally,  partial  lists  of  the 


relevant  leaky  L<uiib  wave  parametcis  ore  \.;uiitaincd  mi  tables  in  Ajipendix  n.  Nlore  detailed 
discussions  of  these  leaky  Lamb  wave  parameters  are  contained  in  Appendix  A  of  chapters 
2  and  4. 

Comparisons  of  I  G/ 1  as  computed  from  Eqs.  (1)  and  (8)  are  displayed  in  Figs,  la  - 
4a.  These  figures  correspond  to  the  oq,  sq,  ai,  and  S'z  leaky  Lamb  waves,  respectively. 

'fhe  solid  line  represents  the  exact  I  Gi  I  while  the  dashed  line  corresponds  to  the 
approximation.  It  is  immediately  evident  that  the  two  cutves  for  I  G/ 1  are  barely 
distinguishable  for  these  leaky  Lamb  waves.  Hence,  Eq.  (8)  appears  to  be  an  excellent 
approximation  of  I  G/ 1  for  these  leaky  Lamb  waves.  The  phases  of  the  exact  coupling 
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from  Eq.  (8)  is  the  dashed  line.  The  agTcenient  between  Eqs.  (1)  and  (8)  is  excellent.  The 
phase  of  the  exact  Gi  is  shown  in  (b).  The  fine  structure  evident  in  arg(Gi)  is  attributed  to 
the  method  of  calculation  of  the  derivative  in  Eq.(l).  The  matetial  of  the  shell  is  440c 
stainless  steel  with  radii  ratio  b/a  =  0.838. 
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coefficients,  arg{Gt)  =  tan-^[Im[Gi}/Re{Gi]],  for  these  Lamb  waves  are  presented  in 
Figs,  lb  -  4b.  These  figures  indicate  that  the  assumption  arg(Gi)  =  0  (/  =  oq,  sq,  ai,  and  S2) 
is  acceptable  for  the  ray  synthesis  calculations.  The  apparent  deviation  in  Figs,  lb  and  2b 
from  arg{Gi) «  0  (/  =  uq,  jq)  for  Jc  <  20  might  be  anticipated  since  the  Watson 
transformation  is  an  asymptotic  analysis.  The  effects  of  the  fluid  loading  of  the  shell  are  a 
plausible  cause  for  the  structure  evident  in  arg(Gi),  I  =  S2,  since  near  a  cutoff  frequency  the 
fluid  loading  may  significantly  affect  the  behavior  of  the  leaky  Lamb  wave.^^  The  rapid 
variation  in  arg(Gi),  I  =  ai,  is  attributed  to  the  method  used  in  computing  the  derivative  of 
with  respect  to  V.  However,  inspection  of  the  ray  synthesis  of  the  form  function  for 
backscattering  in  chapter  4  and  the  total  scattering  cross  section  in  chapter  2  demonstrates 
the  usefulness  of  Eq.  (8)  with  arg(Gi)  =  0. 

The  leaky  Lamb  wave  contribution  to  backscattering  demonstrates  behavior  that 
merits  a  more  detailed  discussion.  Figure  5  displays  I G/ 1  calculated  using  Eqs.  (i)  «md 
(8).  In  Fig.  5,  the  dashed  line  is  the  approximation  and  the  so’- d  curve  is  the  exact  I G/  i. 
Inspection  of  Fig.  5  indicates  that  Eq.  (8)  is  an  excellent  approximation  for  I  G/ 1  when  x  > 
72,  but  it  appears  to  diverge  from  the  exact  coupling  coefficient  for  x  <  72.  The  resolution 
in  the  region  of  the  divergence  is  enhanced  in  Fig.  5b  to  give  greater  detail  of  the  maximum 
in  the  exact  I  G/ 1.  Furthermore,  the  phase  of  the  exact  T/  in  Fig.  6  supports  the  assumption 
argiQi)  =  0  for  r  >  72,  but  arg{Gi)  appears  to  have  a  ~kI2  phase  shift  in  the  vicinity  of  r  = 
72,  The  apparent  divergence  of  Eq.  (8)  from  the  exact  I  G/ 1  and  the  -7t/2  phase  shift  are 
related  to  the  negative  group  velocity  of  the  si  leaky  Lamb  wave.  In  Appendices  A  and  C 
of  chapter  4,  this  leaky  Lamb  wave  is  shown  to  have  a  negative  group  velocity  for  x  <  72. 
Hence,  the  approximation  in  Eq.  (8)  for  I  G/ 1  with  argiGi)  =  0  appears  to  be  applicable  for 
those  surface  guided  elastic  wave  that  have  positive  group  velocity  and  seems  to 
breakdown  when  the  group  velocity  becomes  negative. 


=  .^1,  from  Eq.  (1)  (solid  Line)  and  the 
Figure  5a  demon  stiaies  the  apparent  divergenc  e 
hile  Fig.  5b  enhances  the  resolution  near  die 


Two  important  numerical  techniques  used  in  the  determination  of  the  exact  G/  are 
noted.  First,  the  derivative,  d!Z>+^/dv,  evaluated  at  the  Ith  leaky  Lamb  wave  pole  v  =  V/  is 
performed  by  a  contour  integration.  The  theory  of  residues  from  complex  analysis  gives 


d©v(x) 


dv 


"^V  =  V 


<(x) 


dv, 


t(v-v/ 


(10) 


where  is  assumed  analytic  everywhere  within  and  on  the  contour  C  and  v/  is  contained 
within  C.  For  the  numerical  algorithm  implemented,  C  is  chosen  so  as  to  isolate  the 
residue  of  the  pole  at  v  =  v/  from  any  other  possible  residue  contribution.  While  Eq.  (10) 
appears  to  be  a  good  method  for  the  evaluation  of  d!Z7^/dv  at  v  =  v/,  there  seems  to  be 
regions  where  numerical  round-off  may  corrupt  the  result  from  Eq.  (10).  To  minimize  this 
type  of  possible  numerical  error,  a  three-point  sliding  average  was  implemented  on !  Gi  I 
and  argiGi).  This  sliding  average  tends  to  smooth  any  fluctuations  in  I  G/ 1  or  argipi) 
introduced  from  implementation  of  Eq.  (10). 

B.  Significance  of  the  group  velocity  and  its  sign 

To  better  understand  the  breakdown  of  Eq.  (8),  it  is  first  necessary  to  review  the 
relationship  between  the  group  and  phase  velocities  and  the  spacing  between  adjacent 
resonances.  The  group  velocity  for  the  Ith  leaky  Lamb  wave  Cgi  is  defined  by  Cgi  =  doi/dki 
where  to  is  the  angular  frequency  and  ki  is  the  wavenumber  of  the  Lamb  wave  propagating 
on  the  shell.  The  phase  velocity  of  the  Ith  leaky  Lamb  wave  is  c/  =  (o/ki  and  the  phase 
velocity  of  the  incident  wave  is  c  =  (X)/k.  The  group  and  phase  velocities  for  the  Ith  leaky 
Lamb  wave  are  related  by 
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It  is  generally  accepted  that  for  a  weakly  dispersive  surface  guided  elastic  wave  the  spacing 
between  adjacent  resonance  peaks  is  AXres  ~  Cgllc-  If  die  leaky  Lamb  wave  propagates 
without  dispersion,  then  Eq.  (11)  reduces  to  the  expected  result  Cgi  =  c/. 

The  simple  relationship  uccres  ~  Cg^c  can  be  established  from  the  Fabry-Perot  form 
of  a  leaky  Lamb  w'ave  contribution  to  backscatiering.  If  x-Xm,  {m  =  1,2),  correspond  to 
different  resonances  in //,  then  ;ci  and  X2  are  adjacent  resonance  if  the  phase  <!>(>:)  = 

27ca/(x)  in  the  denominator  of  Eq.  (4)  changes  by  27i  as  x  -  x\  is  increased  to  X2-  The 
condition  for  adjacent  resonances  becomes  [aKx2)  -  a/(jri)]  =  1  where  &Xres  =X2-  xi.  A 
two  term  Taylor  series  expansion  of  a/(.V2)  about  jc  =  xi  gives  the  expression 
AJtrgj  »  [a}(x  =  Xj)]’^  where  terms  of  order  (Axres)^  higher  can  be  neglected  for 

weakly  dispersive  leaky  Lamb  waves.  The  function  a'^(x)  can  be  evaluated  from  Eqs.  (2c) 
and(ll).  Diflerentiationof  Eq.  (2c)givesa/(x)  5=  (c/c/)ll  -  (x/c^)(dc’^/dx)].  FromEq. 
(11),  the  second  term  in  a'^(x)  becomes  (x/ci){dci/dx)  -  {Cgi  -  c/)/c/.  Combining  these 

results  gives  the  simple  expressions 

a'/x)  =  c/c^,  (12a, b) 

Equation  (12b)  has  been  verified  in  Appe  ndix  A  of  chapter  4  for  regions  of  x  where  the  ao 
and  .yq  leaky  L,amb  waves  demonstrate  weak  dispersion. 

The  relevance  of  Eq.  (12a)  to  tlic  leaky  Lamb  wave  contribution  to  backscattering 

and  the  breakdown  of  Eq.  (8)  for  the  coupling  coefficient  may  be  found  by  examining 
ot/(x).  The  instantaneous  slope  of  ai{x)  is  a'/(x)  which  implies  that  the  group  velocity  of  a 

leaky  Lamb  wave  can  be  directly  obtained  frotn  tiie  Watson  iransfonnation  without  the  use 
of  Eq.  (1 1),  Figures  7a  and  7b  are  a/(x)  and  ()/(x)  dcicimined  by  a  Watson  methodology. 
Figure  7a  shows  that  fx}(x)  is  positive  for  x  >  72,  goes  through  zeio  ncai  x  =  72,  and  is 

negative  for  x  <  72.  As.  x  approaches  x  ~  72  from  above  Cgi/c  divergences  to  +  oo  and 
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divergences  to  ~  oo  as  x  approaches  72  from  below.  Hence,  Fig,  7a  confirms  the  group 
velocity  calculation  displayed  in  Fig.  A.3  of  chapter  4  for  the  i’l  leaky  Lamb  wave.  The 
radiation  damping  parameter  P/,  I  =  s\  is  shown  in  Fig.  7b  for  completeness.  Inspection 
of  Fig.  7b  and  Eq.  (4)  demonstrates  that  for  x  <  72  contribution  from  circumnavigation  of 
the  Lamb  wave  about  shell  is  negligible.  The  locus  of  leaky  Lamb  wave  roots  v^(x), 

I  =  SI  is  displayed  in  Fig.  8. 

Finally,  the  width  of  an  isolated  resonances  is  related  to  in  Eq.  (3).  With  the 
assumed  coincidence  of  the  poles  of /«/(x)  and //(x)  in  the  complex  x~plane  Tni  becomes  Yn/ 
as  noted  below  Eq.  (6).  Inspection  of  Eq.  (6b)  suggests,  however,  that  Ynl  would  go  to 
zero  near  x  =  72  as  a'/(x)  vartishes.  Furthermore,  since  a/(x)  becomes  negative  for  x  < 

72,  then  Yn/  would  be  negative  for  any  resonances  in  this  region.  That  is,  the  width  of 
leaky  Lamb  wave  resonances  would  appear  to  be  negative  for  a  leaky  Lamb  wave  with 
negative  group  velocity  if  the  resonance  is  assumed  isolated  having  the  form  of  Eq.  (3). 
When  Cgfc  is  negative,  however,  the  argument  breaks  down  since  the  resoriance  spacing 
from  Eq.  (12b)  is  predicted  to  be  negative.  It  would  appear  that  the  derivation  that  <  0 
also  breaks  down.  Finally,  when  c^//c  <  0,  it  may  also  be  argued  that  the  combined  result 
of  Eqs.  (4)  and  (8)  with  arg{Gi)  -  0  is  incomplete. 

One  caveat  is  noted.  From  Fig.  (7b),  it  is  observed  that  the  .yi  leaky  Lamb  wave 
has  a  rather  large  radiation  damping  parameter.  In  Ref.  2,  the  assumption  27cP/ «  1  was 
used  in  deriving  Eq.  (8).  Subsequently,  this  condition  was  determined  to  be  too  restrictive 
and  Eq.  (8)  has  been  shov/n  to  give  fairly  good  agreerrient  with  the  exact  I  G/ 1  even  when 
27iP/  «  I  is  not  satisfied.^  (Whether  the  group  velocity  is  positive  or  negative  for  the 
relevant  surface  clastic  waves  in  those  comparison  was  not  pursued.]  Furthermore,  the 
exp\  ~2iv.  -  9/)P/]  factor  in  Eq.  (4)  indicates  that  the  .si  leaky  Lamb  wave  contribution  in 
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the  ray  synthesis  for  the  ray  path  assumed  in  Eq.  (4)  is  negligible  for  x<12  since  P/  is 
large. 

Finally,  the  present  level  of  ray  synthesis  appears  to  be  incomplete.  The  current  ray 
model  includes  contributions  from  a  specular  reflection  and  leaky  Lamb  waves  which 
circumnavigate  the  elastic  spherical  shell.  The  generalized  specular  reflection  contains  the 
ordinary  specular  ray  and  rays  associated  with  repeated  internal  reflections  from  the  inner 
surface  of  the  shell  at  r  =  6  (see  Sec.  2  of  Chapter  4).  The  inclusion  of  these  internal  rays 
seems  to  correctly  describe  the  manifestation  of  longitudinal  resonances  in  the  form 
function  for  backscattering.  The  ray  synthesis  of  the  Ith  leaky  Lamb  wave  contribution 
giving  a  Fabry-Perot  form  accounts  for  the  propagation  of  the  leaky  Lamb  wave  around  the 
shell,  but  neglects  the  possibility  of  a  direct  ray  contribution.  A  direct  contribution  from  a 
leaky  Lamb  wave  is  briefly  discussed  in  Sec.  4  of  Chapter  4  and  is  illustrated  by  ray  jB"A" 
in  Fig.  1  of  Chapter  4  (where  0/  0).  The  results  of  a  Sommerfeld-Watson 

transformation  of  the  form  function  for  backscattering  from  an  elastic  sphere,'  which  are 
extended  to  the  shell,  are  based  on  the  assumption  a/ »  P/.  This  assumption  may  omit 
direct  ray  contributions  associated  with  a  leaky  Lamb  wave  when  P/  ~  a/.  Furthermore, 

CTW  All  OX/&tlV/  TT  JT  TTAUA  V/A  AAVA^W  A  UWAAiAkAV/A  I 

damping.^  It  is  noteworthy  that  the  use  of  the  revised  coupling  coefficient  shown  in  Figs. 

5  and  6  does  not  remove  the  anomaly  near  -  71  in  Fig.  7  of  Chapter  4.  That  is  because 
l/n  I  is  small  since  p/  is  large  (p/  >  2)  near  this  anomaly. 
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Appendix  H 

Leaky  Lamb  wave  parameters 

The  contents  of  this  appendix  are  five  tables  containing  the  parameters  used  in  the 
ray  synthesis  calculations  presented  in  Chapters  2  and  4  and  Appendix  I.  The  method 
employed  in  determining  these  parameters '  *line  in  Appendix  A  of  Chapter  2  or  4.  The 
format  of  each  table  is  four  columns  of  numbers  for  the  dimensionless  size  parameter  ka, 
normalized  phase  velocity  of  the  leaky  Lamb  wave  ci/c,  the  radiation  damping  parameter  P/, 
and  a/  =  /?e(v/)  from  the  Watson  transform  methodology.  These  tables  are  for  a  440c 
stainless  steel  shell  with  inner-to-outer  radii  ratio  bla  =  0.838.  The  material  parameters  for 
440c  stainless  steel  are  given  in  Sec.  5  of  chapter  2. 

The  computer  code  used  in  determining  v/  for  the  Ith  leaky  Lamb  wave  at  a  fixed  x 
is  essentially  the  same  as  the  FORTRAN  program  ZNU  implemented  by  W  illiams^  for 
solid  elastic  spheres.  This  program  is  based  on  the  winding  number  theory  from  complex 
analysis  which  is  discussed  by  Williams.  Two  modifications  were  made  to  ZNU  for  the 
computation  of  V/  for  ^e  evacuated  spherical  shell.  First,  the  external  function  F(v)  is  the 
determinant  D^ix)  in  the  denominator  of  the  partial-wave  coefficient  where  the  integer  index 
n  is  replaced  everywhere  by  complex  v.  For  spheres,  D^ix)  is  a  3  x  3  determinant  and  for 
evacuated  shells  D\/(x)  becomes  a  5  x  5  determinant.  The  elements  of  Df^{x)  for  a  shell  are 
listed  in  Appendix  A  of  Ref.  2.  Second,  the  routine  OLVER  that  is  employed  for  the 
computation  of  Bessel  and  Neumann  functions  of  real  argument  and  com^'lex  order  has 
been  slightly  rewritten.  Previously,  this  routine  returned  single  precision  values,  but  since 
all  other  computations  in  ZNU  (and  in  OLVER)  are  double  precision  OLVER  was  modified 
to  return  double  precision  values.  Ore  other  modification  to  OLVER  consisted  of  the 
rewriting  of  portions  of  the  code  tiiat  are  in  FORTRAN  IV  in  standard  FORTRAN  77. 


Once  a  root  v/  at  a  fixed  x  for  the  Uh  leaky  Lamb  wave  was  determined,  it  was 
inserted  into  a  secon  ^  computer  code  that  tracked  the  locus  v/(jc).  A  discussion  of  the 
necessary  mathematical  tlieory  is  given  in  Appendix  B  of  Ref.  2.  The  algorithm  assumed 
that  V/fx)  is  a  well-behaved  function  of  x  (see  Fig.  A1  of  Chapter  2)  and  that  v/(x)  was 
approximately  the  value  of  V/  at  x+Ax  where  Ax  «  1.  A  square  contour  was  centered 
about  V/(x),  X  was  increased  (or  decreased)  by  Ax  and  then  the  winding  number  was 
determined  from  numerical  integration  of  Eq.  (B 1)  in  Ref.  2.  If  the  winding  number  is 
equal  to  one  (an  isolated  zero,  v  =  V/,  of  D^(x)  is  enclosed  within  the  contour),  then  Eq, 
(B2)  from  Ref.  2  was  integrated  about  die  same  contour.  This  recursive  procedure  was 
implemented  for  the  determination  of  the  parameters  given  in  Tables  n.l  - 11.5. 
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Table  II.  1  The  ao  Ifaky  Lamb  wave  parameters 

ka  ci/c  a/  P/  ka  ci/c  a/  P/ 


7.00 

1.084072 

5.957134 

0.383446 

9.00 

1.199730 

7.001686 

0.319891 

11,0 

1.293156 

8.006322 

0.276207 

13,0 

1.372558 

8.971365 

0.248181 

15.0 

1.441409 

9.906483 

0.229505 

17.0 

1.501791 

10.81982 

0.2)6644 

19.0 

1.555182 

11.71722 

0.207613 

21.0 

1.602686 

12.60301 

0.201267 

23.0 

1.645162 

13,48039 

0.196897 

25.0 

1.683309 

14.35170 

0.194052 

27.0 

1.717691 

15.21878 

0.192434 

29.0 

1.748780 

16.08299 

0.191845 

31.0 

1.776967 

16.94545 

0.192153 

33.0 

1.802577 

17.80712 

0.193257 

35.0 

1,825891 

18.66873 

0.195096 

37.0 

1.847151 

19.53084 

0.197632 

39.0 

1.866561 

20.39404 

0.200840 

41.0 

1.884300 

21.25874 

0.204714 

43.0 

1.900525 

22.12533 

0.209249 

45.0 

1.9i5362 

22.99425 

0.214458 

47.0 

1.928935 

23.86577 

0.220342 

49.0 

1.941352 

24.74014 

0.226932 

51.0 

1.952696 

25.61774 

0.234247 

8.00 

1.145666 

6.482839 

0.349988 

10.0 

1.248506 

7  509575 

0.295577 

12.0 

1.334342 

8.493196 

0.260708 

14.0 

1.408150 

9.442120 

0.237946 

16.0 

1.472554 

10.365476 

0.222498 

18.0 

1.529290 

11.27017 

0.211732 

20.0 

1.579610 

12.16135 

0.204157 

22.0 

1.624505 

13.04259 

0.198866 

24.0 

1.664737 

13,91669 

0.195304 

26.0 

1.700941 

14.78566 

0.193103 

28.0 

1.733624 

15.65114 

0,192022 

30.0 

1.763213 

16.51439 

0.191896 

32.0 

1.790076 

17.37634 

0.192609 

34.0 

1.814507 

18.23787 

0.194087 

36.0 

1.836764 

19.09969 

0.196279 

38.0 

1.857077 

19.96226 

0.199153 

40.0 

1.875630 

20.82616 

0.202694 

42.0 

1.892595 

21.69176 

0.206898 

44.0 

1.908108 

22.55949 

0.211769 

46.0 

1.922300 

23.42967 

0.217318 

48.0 

1.935281 

24.30260 

0.223548 

50.0 

1.947150 

25.17855 

0.230503 

52.0 

1.957989 

26,05786 

0.238197 
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53.0 

1.963046 

26,49886 

0.242336 

54.0 

1.967873 

26.94079 

0.246672 

55.0 

1,972478 

27.38371 

0.251207 

56.0 

1.976872 

27.82759 

0.255943 

57.0 

1.981057 

28.27252 

0.260883 

58.0 

1.985041 

28.71854 

0.266031 

59.0 

1,988834 

29.16563 

0.271386 

60.0 

1.992438 

29.61386 

0.276951 

61.0 

1.995861 

30.06326 

0.282730 

62.0 

1.999109 

30.51382 

0.288721 

63.0 

2.002189 

30.96557 

0.294924 

64.0 

2.005101 

31.41859 

0.301338 

65.0 

2.007853 

31.87290 

0.307964 

66.0 

2.010449 

32.32848 

0.314796 

67.0 

2.012896 

32.78538 

0.321832 

68.0 

2.015195 

33.24363 

0.329068 

69.0 

2.017353 

33.70323 

0.336501 

70.0 

2.019376 

34.16418 

0.344117 

71.0 

2.021263 

34.62655 

0.351917 

72.0 

2.023021 

35.09033 

0.359887 

73.0 

2.024654 

35.55553 

0.368019 

74.0 

2.026168 

36.02214 

0.376265 

75.0 

2.027564 

36.49020 

0.384720 

76.0 

2.028849 

36.95966 

0.393266 

77.0 

2.030026 

37.43055 

0.401925 

78.0 

2.031099 

37.90285 

0.410681 

79.0 

2.032071 

38.37660 

0.419523 

80.0 

2,032949 

38.85170 

0.428427 

81.0 

2.033733 

39.32823 

0.437391 

82.0 

2.034432 

39.80610 

0.446391 

83.0 

2.035046 

40.28531 

0.455415 

84.0 

2.035585 

40.76577 

0.464444 

C\  e’  /\ 

o:j.u 

/JO 

0,4/34/  i 

n  y'  r\ 

oo.u 

r%  y'  1^  y' 

.<i.U.3W.30 

14  r^y^rxy'  1 

/\  1 

v;,4oz4/  / 

87.0 

2.036762 

42.21485 

0.491445 

88.0 

2.037024 

42.70028 

0.500374 

89.0 

2.037226 

43.18685 

0.509245 

90.0 

2.037374 

43.67450 

0.518053 

91.0 

2.037473 

44.16318 

0,526782 

92.0 

2.037521 

44.65290 

0,535432 

93.0 

2.037526 

45.14359 

0.543993 

94.0 

2.037491 

45.63517 

0.552456 

95.0 

2.037417 

46.12766 

0.560822 

96.0 

46.62099 

0.569089 

97.0 

2.037169 

47.11510 

0.577247 

98.0 

2.037001 

47.60994 

0.585298 

99.0 

2.036806 

48.10551 

0.593243 

100.0 

2.036586 

i  48.60179 

0.601084 

Table  II.2  The  Lamb  wave  parameters 


ka 

ci/c 

a/ 

P' 

ka 

ci/c 

a/ 

Pi 

7.0 

5.672502 

0.734023 

0.257460 

8.0 

4.768751 

1.177.588 

9.36004e-02 

9.0 

4.411582 

1.540085 

5.16725e-02 

10.0 

4.235214 

1.861156 

3.24039e-02 

11.0 

4.131961 

2.162174 

2.14516e-02 

12.0 

4.065046 

2.451996 

1.45435e-02 

13.0 

4.018514 

2.735027 

9.91413e-03 

14.0 

3.984383 

3.013719 

6.6948  le-03 

L5.0 

3.958239 

3.289564 

4.4{)952e-03 

16.0 

3.937482 

3.563510 

2.7795  le-03 

17.0 

3.920461 

3.836225 

1.62705e-03 

18.0 

3.906104 

4.108173 

8.40220e-04 

19.0 

3.893684 

4.379697 

3.42250e.-04 

20.0 

3.882687 

4.651071 

7.97300e-05 

21.0 

3.872742 

4.922514 

1.59100e-05 

22.0 

3.863574 

5.194210 

1.24970e-04 

23.0 

3.854973 

5.466319 

3.89240e-04 

24.0 

3.846779 

5.738986 

7.96640C-04 

25.0 

3.838863 

6.012345 

1.33805e-03 

26.0 

3.831120 

6.286527 

2.01034e-73 

27.0 

3.823465 

6.561658 

2.8il40e-03 

28.0 

3.815824 

6.837864 

3.742 18e-03 

29.0 

3.808135 

7.115277 

4.80525e-03 

30.0 

3.800339 

7.394032 

6,00586e-03 

31.0 

3.792385 

7.674275 

7.35021e-03 

32.0 

3.784225 

7.956158 

8.84729e-03 

33.0 

3.775811 

8.239845 

1.05080e-02 

34.0 

3.767098 

8.525516 

1.23452e-02 

35.0 

3.758041 

8.813362 

1.437456-02 

36.0 

3.748588 

9.103615 

1.66i37e-02 

37.0 

3.m693 

9.396507 

1.908346-02 

38.0 

3.728302 

9.692308 

2.18088e-02 

39.0 

3.717362 

9.991310 

2.48179e-02 

40.0 

3.705806 

10.293874 

2.81443e-02 

41.0 

3.693569 

10.60038 

3.18271e-02 

42.0 

3.680577 

10.91125 

3.59l08e-02 

43.0 

3.666747 

11.22702 

4.04499e-02 

44.0 

3.651984 

1 )  .54825 

4.55070e-02 

45.0 

3.636183 

11.87561 

5.11561e-02 

46.0 

3.619230 

12.20989 

5.74824e-G2 

47.0 

3.600981 

12.55200 

6.459 16e- 02 

48.0 

3.581286 

12.90301 

7.26029e-02 

49.0 

3.559972 

13.26416 

8.16586e-02 

50.0 

3.536840 

13.63691 

9.l9230e-02 

51,0 

3.511670 

14.02301 

0.103584 

52.0 

3.484216 

14.42445 

0.116846 
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53.0 

3.454221 

14.84355 

0.1  U  937 

54.0 

3.421392 

15.28305 

0,149076 

55.0 

3.385458 

15.74595 

0.168431 

56.0 

3.346179 

16.23551 

0.190073 

57.0 

3.303382 

16.75504 

0.213381 

58.0 

3.257042 

17.30757 

0.239423 

59.0 

3.207344 

17.89528 

0.265850 

60.0 

3.154775 

18.51879 

0.291861 

61.0 

3.100104 

19.17676 

0.315840 

62.0 

3.044410 

19.86520 

0.336147 

63.0 

2.988909 

20.57793 

0.351540 

64.0 

2.934767 

21.30753 

0.361481 

65.0 

2.882950 

22.04630 

0.366147 

66.0 

2.834155 

22.78737 

0.366245 

67.0 

2.788750 

23.52510 

0.362730 

68.0 

2.746857 

24.25557 

0.356537 

69.0 

2.708440 

24.97592 

0.348488 

70.0 

2.673345 

25.68442 

0.339227 

71.0 

2.64x345 

26.38025 

0.329233 

72.0 

2.612192 

27.06305 

0.318843 

73.0 

2.585635 

27.73292 

0.308290 

74.0 

2.561434 

28.39007 

0.297732 

75.0 

2.539353 

29.03508 

0.287265 

76.0 

,2.519193 

29.66839 

0.276955 

77.0 

2.500764 

30.29058 

0.266844 

78.0 

2.483897 

30.90227 

0.256955 

79.0 

2.468440 

31.50402 

0.247305 

80.0 

2.454265 

32.09632 

0.237902 

81.0 

2.441245 

32.67980 

0.228745 

82.0 

2.429277 

33.25490 

0.219838 

83.0 

2.418263 

33.82216 

0.211181 

84.0 

2.408122 

34.38195 

0.202775 

85.0 

2.398775 

34.93475 

0.194616 

86.0 

2.390154 

35.4K094 

0.186707 

87.0 

2.382201 

36.02085 

0.179049 

88.0 

2.374851 

36.55495 

0.171636 

89.0 

2.368062 

37.08347 

0.164468 

90.0 

2.361786 

37.60675 

0.157544 

91.0 

2.355984 

38.12505 

0.150864 

92.0 

2.350614 

38.63871 

0.144424 

93.0 

2.345643 

39.14797 

0.138224 

94.0 

2.341044 

39.65302 

0.132259 

95.0 

2.336785 

40.15415 

0.126526 

96.0 

2.332840 

40.65155 

0.121020 

97.0 

2.329188 

41.14542 

0.115739 

98.0 

2,325805 

41.63594 

0.110679 

99.0 

2.322670 

42.12336 

0.105831 

100.0 

1  2.319769 

42.60775 

0.101101 
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Table; 

n.3  llieai 

[  leaky  Lamb  wave  parameters 

ka 

cilc 

a/ 

P/  ka 

ci/c 

ai 

P/ 

43.0  16.73859  2.068915  1.28650e-02 


42.0 

28.71326 

0.962739 

6.85368e-03 

44.0 

13.12078 

2.853460 

1.66242e-02 

46.0 

10.014726 

4.093236 

2.17523e-02 

48.0 

8.525793 

5.129975 

2.52l79e-02 

50.0 

7.621105 

6.060729 

2.76363e-02 

52.0 

7.002441 

6.925982 

2.92020e-02 

54.0 

6.547668 

7.7^7211 

3.00173e-02 

56.0 

6.196069 

8.537988 

3.00869e-02 

58.0 

5.913572 

9.307947 

2.93676e-02 

60.0 

5.679032 

10.065181 

2.777 12e-02 

62.0 

5.478043 

10.81791 

2.51732e-02 

64.0 

5.299794 

11.57594 

2.14448e-02 

66.0 

5.135036 

12.35288 

1.65774e-02 

68.0 

4.975240 

13.16768 

1.10206e-02 

70.0 

4.814225 

14.04024 

6.16345e-03 

72.0 

4.654351 

14.96940 

3.66768e-03 

74.C 

4.507755 

15.91615 

3.75392e-03 

76.0 

4.383449 

16.83795 

6.06924e-03 

78.0 

4.281067 

17.71976 

1.05131e-02 

80.0 

4.19^^21 

18. .56568 

1.68916e-02 

82.0 

4.123632 

19.38538 

2.49176C-02 

84.0 

4.060276 

20.18825 

3.43295e-02 

86.0 

4.003290 

20.98233 

4.49350e-02 

45.0 

11.22238 

3.509846 

1.94656e-02 

47.0 

9.164049 

4.628737 

2.36409e-02 

49.0 

8.025656 

5.605420 

2.65350e-02 

51.0 

7.285776 

6.499941 

2.85153e-02 

53.0 

6.759251 

7.341106 

2.970 18e-02 

55.0 

6.361508 

8.145749 

3.01466e-02 

57.0 

6.047717 

8.925045 

2.9830  le-C2 

59.0 

5.791311 

9.687675 

2.86860e-02 

61.0 

5.575073 

10.441560 

2.66050e-02 

63.0 

5.386673 

11.19553 

2.34563e-02 

65.0 

5.216276 

11.96100 

1.91423e-02 

67.0 

5.055015 

12.75416 

1.38237e-02 

69.0 

4.895021 

13.59595 

8.38266e-03 

7L0 

4J33506 

14,49945 

4.56789e-03 

73.0 

4.578581 

15.44380 

3.4l507e-03 

75.0 

4.442661 

16.38177 

4.639 19e-03 

77.0 

4.329773 

17.28384 

8.03265e-03 

79.0 

4.236684 

18.14666 

1.34778e-02 

81.0 

4.158494 

18.97820 

2.07165e-02 

83.0 

4.091011 

19.78839 

2.94636e-02 

85.0 

4.031121 

20.58595 

3.94923e-02 

87.0 

3,976569 

21.37816 

5.06404e-02 
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88.0 

3.950755 

21.77423 

5.65988e-02 

90.0 

3.9012.31 

22.56964 

6.923 14e-02 

92.0 

3.8536.50 

23.37347 

8.27627C-02 

94.G 

3.807183 

24.19017 

9.71307C-02 

96.0 

3.76  U  64 

25.02401 

0.112258 

98.0 

3.715079 

25.87898 

0.128037 

100.0 

3.608520 

26.75895 

0.144322 

89.0 

3.925690 

22.17118 

6.27988e-02 

91.0 

3.877255 

22.97021 

7.5S885e-02 

93.0 

3.830321 

23.77995 

8.98488e-02 

95.0 

3.784153 

24.60469 

0.104604 

97.0 

3.738158 

25.44861 

0.120075 

99.0 

3.691880 

26.31561 

0.136129 
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Table 

a.4  The  51 

leaky  Lamb  wave  parameters 

ka 

ciic 

at 

P/  ka 

ci/c 

at 

P/ 

65.0 

5.709630 

10.88427 

7.561313 

66.0 

5.907279 

10.67266 

6.900472 

67.0 

6.116464 

10.454041 

6.191272 

68.0 

6.327767 

10.246286 

5.421073 

69.0 

6.536299 

10.056^132 

4.554775 

70.0 

6.723228 

9.911665 

3.525174 

71.0 

6.797864 

9.944458 

2.161037 

72.0 

6.254679 

11.01138 

0.693661 

73.0 

5.758624 

12.17664 

0.355892 

74.0 

5.498012 

12.95941 

0.248617 

75.0 

5.331644 

13.56695 

0.194698 

76.0 

5.213098 

14.07866 

0  162056 

77.0 

5.122361 

14.53213 

0.140130 

78.0 

5.049326 

14.94761 

0.124354 

79.0 

4.988286 

15.33710 

0.112390 

80.0 

4.935832 

15.70801 

0.102943 

81.0 

A  QQmAn 

t  *T  t 

1  C 

n  1  Aa^  M 

^  A.  *  -T  y  W- 

o 

A 

U“l  U  U  r' 

:  /1 1 ^ 

Jt  \l,“|  *  4  1 

Q  M 

83.0 

4.811332 

16.75094 

8.30929e-02 

84.0 

4.777253 

17.08333 

7.81386C-02 

85.0 

4.745802 

17.41056 

7.36lj7(k--02 

80.0 

4.716574 

17.73357 

6  96259e-02 

87.0 

4.689251 

18.05307 

6.58708C-02 

88.0 

4.663565 

18.36968 

6.23563C-02 

89.0 

4.639316 

18.68386 

5.90374e-02 

90.0 

4  616330 

18,99601 

5.5878Oe-02 

91.0 

4.594466 

19..30643 

5.2849Ge-02 

92.0 

4.573600 

19.61545 

4.993 15e-02 

93.0 

4.553627 

19.92328 

4.71072e-02 

94.0 

4.534459 

20.23015 

4.43636e-02 

95.0 

4.516011 

20.53626 

4.16917e-02 

96.0 

4.498212 

20,84181 

3.90837e-02 

97.0 

4.481091 

21.14695 

3.65346e-02 

98.0 

4.464319 

21.45184 

3.40403e-02 

99.0 

4.448117 

21.75661 

3.15994e-02 

100.0 

4,432345 

22.06142 

2.92110e-02 
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Table  II. 5  Tlic5'2 

leaky  Lamb  wave  parameters 

ka  ci/c 

ai 

P/  ka 

cjc 

ai 

P/ 

83.0  87.25604  0.451224  4.12350e-02 
85.0  31.51721  2.196939  0.114933 
87.0  21.84771  3.482111  0.137002 
89.0  17.32399  4.637385  0.147201 
91.0  14.62990  5.720140  0.152495 
93.0  12.82195  6.753189  0.155372 
95.0  11.51684  7.748792  0.156940 
97.0  10.526631  8.714724  0.157755 
99.0  9.747499  9.656452  0.158121 


84.0  43.13697  1.447286  9.23560e-02 
8C.0  25.57552  2.862591  0.128295 
88.0  19.25171  4.071021  0.142974 
90.0  15.82827  5.186030  0.150257 
92.0  13.64587  6.241966  0.154148 
94.0  12.12103  7.255115  0.156276 
96.0  10.99016  8.235086  0.157418 
98.0  10.115261  9.188332  0.157981 
100.0  9.416545  10.119606  0.157979 
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